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SYNTHESIS OF OSCILLATING 
ADAPTIVE FEEDBACK SYSTEMS 

Abstract — It has been known for many years that insertion 
of certain types of nonlinear elements in the forward 
loop of a feedback system will produce a limit cycle and 
give the system the laudable property of zero sensitivity 
to plant gain changes. Some such systems have been built 
and tested, but no genuine synthesis theory has appeared. 
In this paper a synthesis theory is developed which 
allows system design to proceed from practical specifi- 
cations on system command and/or disturbance response to 
a design which is very nearly optimal in terms of feedback 
sensor noise effects. The approach taken is to replace 
the nonlinear element by a mean square error minimizing 
approximation (dual-input describing function ) , and then 
use linear frequency domain synthesis techniques subject 
to additional constraints imposed by the limit cycle and 
the approximator. Synthesis techniques are also devel- 
oped for a similar system using an externally excited 
oscillating signal with the above approach. 

The results, to a large extent, remove the design of 
the systems considered from the realm of simulation and 
experimentation, permitting true synthesis and the opti- 
mization that accompanies it. 
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PRELIMINARY BACKGROtJND FOR OSCILLATING 
SYSTEM DESIGN 

1.1 In troduction 

The limit cycling or self-oscillating adaptive sys- 
tem (SOAS) has appeared in the literature several times 

over the preceding decade and before . Indeed one 

. f7l 

version of an oscillating system (externally excited 

and somewhat different from those in this paper, but 
with the same adaptive philosophy) was produced by 
Minneapolis-Honeywell and successfully flight tested in 
the F-94C aircraft in the role of pitch rate controller, 
and separately in the role of pitch attitude controller. 

An oscillating system was also used on the X-15 
aircraft '■ . The adaptive properties of oscillating sys- 

tems have thus been proven to have useful practical appli- 
cation. In fact, out of many adaptive schemes that have 
been investigated, the oscillating system is one of few 
that have been shown to have practical application to 
systems where there is large plant parameter uncertainty 
or variation. The primary attractiveness of oscillating 
systems arises from their inherent zero sensitivity to 
plant gain factor changes. Even so, as pointed out in 
[6] , a careful study must be made to determine whether. 
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in a given problem, a proposed oscillating system is 
indeed superior to a purely linear system. 

Noticably absent from the literature, however, are 
specific synthesis techniques. Rather, the existing 
theory has been presented on a largely qualitative basis 
and actual hardware design entails much experimentation 
and simulation. Synthesis techniques are needed that 
permit the designer to proceed from a realistic set of 
specifications through a step by step design procedure 
which insures that these specifications are met, while 
optimizing in some sense over the free parameters. In 
addition, the designer should be able to see any trade- 
offs involved in the synthesis steps, and to succinctly 
determine what advantages his design has over other 
possible designs, e»g., a linear design. 

In this paper a systematic synthesis procedure is 
presented. The criterion for optim.ization adopted is 
to minimize the effect of output sensor noise at the 
plant input (or the system output) . Rather broad system 
performance specifications are assumed, and from these 
a design that is quite close to optimum is developed. 

The systematic procedure for approaching this nearly 
optimum design is given and at each step the designer has 
a good "feel" for what his chosen degree of optimization 
is costing in terms of system complexity. 

In Chapter I we first consider the class of non- 
linear elements that are suitable for oscillating systems. 
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how these elements will be mathematically described in 
synthesis, and the constraints imposed by this descrip- 
tion. Next, a detailed statement of the synthesis 
problem and the type of specifications to be satisfied 
is given. Finally, the criterion for optimizing the 
oscillating system design is discussed. 

1.2 Description of the Nonlinear Element in an 

Oscillating System 

In the following, some basic properties of a class 
of nonlinearities will be enumerated and the synthesis 
theory of this paper will then apply to the class of non- 
linear elements possessing these properties. The multi- 
ple-input describing function is used for analysis. The 
multiple- input describing function has been treated ra- 
ther thoroughly and elegantly by Gelb and Vander Velde 
and the abbreviated derivation below follows that work 
closely. 

The nonlinearity in Fig. 1.1 has input x(t) com- 
posed of three known signal types. We shall use in 
place of the nonlinear element N a quasi- linear approxi- 
mator for each of the signal components and require that 
the approximators be such that the total mean squared 


error 



E[e^(t)] = E[(y(t) ~ y(t))^] 
be minimized. By quasi-linear it is meant that for a 



Figure 1.1 Representation of nonlinear element 
for derivation of multiple-input describing 
function. 


fixed input signal each of the h^(t) in Fig. 1.1 per- 
form as a linear system to the particular signal type for 
which it is derived, however, an may depend on some 

parameters of its input signal type, e.g., magnitude, 
and therefore is not truly linear. We shall consider 
only time- invariant N and stationary inputs x(t) so 
the resultant h^^ are thus time-invariant. Forming the 
error, minimizing via the calculus of variations, and 
requiring that the input signals be statistically inde- 
pendent, the following decoupled set of equations in h. 
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is obtained. 


o X 


(t"X)di 


yx. 


(t) 


T > 0 , 


( 1 , 1 ) 


In this equation R (x) is the autocorrelation for 

^i 

input X. (t) and R (x) is the cross-correlation be- 
2 - Y-^j_ 

tween nonlinearity output y(t) and input component 

x^(t) . Eq. 1.1 has identically the form of the Wiener- 

Hopf integral of linear mean square filtering theory. 

Letting = A sin(w^t+0) where 0 is the 

random variable distributed loniformly over [0,2Tr] radians 

2 

gives R (x) = hA cos w x. Then 
^ o 

o 


2 po 

h^(t)cosi;co^(t-T) ]dt = 


|cosw^'d h^ (t) cosw^tdt + sinw^rj h^ (t) sinw^tdt j . (1.2) 


R „ (t) = E[y(t)A sin[w (t-x)+0]] = AE [y (0) sin [ 0-w x] ] = 

y2i O O 

O 

A cosw^x E[y(O)sin0] - Asinw^x E[y(O)cos0]. Equating the 
last expression to (1.2) for all x > 0 yields 



(t)coso) tdt = — E[y(0)sinGJ 

O X* 



(t)sinio^tdt = 


- ^ E[y (0) COS0] . 


(1.3a) 

(1.3b) 


These equations are satisfied by 
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h (t) = |[E[y(0)sine]6(t) + ~ E [y (0) cos0] 6 (t)j , (1,4a) 

— o 

or upon taking the Fourier transform 

" |[E[y(O)sin0] + jE [y (0) cos0]J . (1,4b) 

is the describing function for the sinusoidal compo- 
nent of input at frequency If the total input 

x(t) is taken to be a single sinusoid in computing 
(1.4b) yields the more conventional describing function 
obtained using harmonic analysis of y(t). The expres- 
sion here is valid for the sinusoidal component at 
of any stationary, but otherwise arbitrary, input pro- 
vided the sinusoidal component and the remainder of the 
input are independent. If there are, for example, two 
sinusoidal components at different frequencies and 

^1 describing function for each may be obtained 

from (1.4 d) by two separate calculations. The two gains 
thus obtained are known as the two sinusoid input de- 
scribing function (TSIDF) . has quite simple form, 

being a complex gain, but calculation of the expectations 
in (1.4 d) is not in general simple. It is shown in [9] 
that is real for static and single-valued non- 

linearities . 

Taking x_(t) = b, a constant, gives R (t) = b^ 

^f 

- bE[y(t)I = bE[y(0)]. Thus using (1.1) 

again 
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hj(t)dt = 


(1.5) 


which j,s satisfied by 


hf (t) 


E[y(0)3 

b 


6(t) 


(1,6a) 


and this represents a pure gain in the frequency domain, 
viz. , 


N. = 


E[y(0) 3 

b 


(1.6b) 


The nonlinearity must have no constant output component 
when there is no constant input (static nonlinearities 
must be odd) in order that be finite as b ^ 0. 

This condition is always assumed in subsequent discussion. 
If the input to N is restricted to be the sum of a 
sinusoid plus a constant signal, the two gains 
are known as the dual-input describing function (DIDF) . 

For use in an oscillating system we shall require that 

Nq == M^A (1,7a) 

Nf - M^/A (1.7b) 

for A/|b[ > 1 by a sufficient margin, where 
do not depend on x(t) . Actually, we require even more. 
Instead of just a constant we shall allow x^(t) to be 
some transient function of time, therefore not stationary, 
and require it to be slow and small with respect to the 



sinusoidal input in the following sense. 


^t^|Xf(t)[ < I^j/a (1.8a) 

% 1 ^ (1,8b) 

where is the bandwidth of X^(jco) and is the 

minimum value that A can assume in the given applica- 
tion. If for a given nonlinearity there exists a, 6 
such that when (1.8) holds, then (1.7) gives a suffi- 
ciently good approximation to the DIDF with a,B,M , and 

o 

dependent only on the physical parameters of the non- 
linearity, i.e., independent of input provided the input 
has the assumed form, then such a nonlinearity is a 
candidate for use in an oscillating system. It will be 
seen that (1.7) is necessary to obtain the zero gain 
sensitivity property in oscillating systems. It is 
stated in [10] that a > 3 and 3 > 3 give less than 5% 
error in the approximations of (1,7) for an ideal relay 
v/hen x^(t) is a sinusoid and this reference indicates 
these values hold closely for a v;ide range of common 
nonlinearities , The error for any given, nonlinearity 
can be obtained by expanding the TSIDF, obtained from 
(1.4b), and DIDF in Taylor series in b/A with b as 
the bias magnitude or the magnitude of the lower fre- 
quency sinusoid. 

It may appear that by now we have restricted our- 
selves to a nearly empty class of nonlinearities, but 



this is not the case. In order that (l„7a) hold it is 
seen that N must have a saturating property, or for 
large A the output sinusoidal component must approach 
Im^I in magnitude, independent of A. In [11] it is 
shown that for all odd, static, and single-valued non- 
linearities, if the input is the sum of tx«70 sinusoids 
x(t) = x^(t) + x^(t) = bcosWj^t + Acosw^t with 
irrational then 

Urn N (A,b) = N (A,0) + ~ -- 5^7 . (1.9) 

b^O aa 

The limit holds also for ojj^ = 0, i.e., when x^(t) = b. 

Of course dN /dA must exist. The above limit has been 
o 

called the incremental input describing function (IIDF) . 
Applying this limit to the right side of (1.7a) we get 
that =2 ~ M^/2A when is a sufficiently 

small constant. Thus we can conclude that every odd, 
saturating, static, and single-valued nonlinearity is in 
the class of candidates for oscillating systems. In a 
practical design situation one must of course determine 
that (1.7) is a satisfactory approximation with accept- 
able values of a and S. We remark also that there 
may be some lower bound set by the nonlinearity on A^ 
in (1.8a), e.g., for a saturation nonlinearity A^^ must 
be large enough to cause saturation. The above in no 
way limits the possible nonlinearities that may prove 
satisfactory for an oscillating system. Smisek^^®^ 


has 
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studied a dynamic nonlinearity for which is complex 

^vith phase lead, and used this element in an SOAS simu- 
lation. The above conditions, viz., nonlinearities that 
are odd, saturating, static, single-valued, and such 
that dN^/dA is finite are theoretically sufficient for 
(1.7) to hold for some a and 3, but they are not 
necessary. 

The IIDF has another useful application. Since it 

holds for any not rationally related to w it can be 

o 

used to predict stability of the forced signal loop 
transmission for small perturbations. For large signal 
stability the TSIDF can be used. 

Next a random process input x^(t) to the nonlinear 
element is considered. For this input R (t) = 

E[y(t)x^(t-T) ] and Eg. 1,1 becomes 

^OO 

I oh^(t)Rj^r (t); t > 0. (1.10) 

For the general nonlinearity R (x) is cruite difficult 

to obtain, and even if it were known the solution of 
(1.10) is not apparent except in special cases, e.g,, if 
the Fourier transform of R (x) is knov/n and rational 
the Wiener-Hopf type integral can be solved by spectral 
factorization. Hov;ever, as shown in [9], for the special 
case of gaussian zero mean 2 !:^(t) (if it is not zero 
mean the constant component may be included in x^) , and 



for static single-valued nonlinearities only^ R (x) - 

*' r 

2 

R (x) E [y (0) X (0) ]/a o o is the standard deviation of 

X r 

r 

the process amplitude distribution. In this case (1.10) 
is seen to be satisfied by 

h^(t) = E[y(0)x^(0)]6(t)/c^ (1.11a) 

and 

- E[y(0)x (0)]/a^ . (1.11b) 

r ” r 

Thus the random input describing function (RIDF) com- 
ponent is a pure gain under the stated conditions. 

Finally it is shown in [9] that with input x(t) = 
bcosoijj^t + Acosoj^t + x^(t), an odd single-valued non- 
linearity, and x^ gaussian, that 

lim N (A,b,o) = N (A,0,a) . (1.12) 

b->0 ^ 

Eq. 1.12 holds when = 0 as well. This suggests that 
for small inputs x^(t) the mean square error minimizing 
quasi- linear approximator is a gain independent of the 
waveshape of x^(t) . It is found also that if in addi- 
tion to (1.8) 

O < Aj^/a^ (1.13) 


then 
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More importantly, we want (1.7) to hold for forced 
and oscillating signals when there is also a noise input 
to N, i.e., the forced and oscillating signal components 
of the three input describing function, N^(A,b,a) and 
must approach the asymptotic values of (1.7). 
This, of course, occurs when o is vanishingly small. A 
check for two common nonlinearities, ideal relay and 
saturation, indicates that (1.7) holds for = a > 3 . 

The nonlinearity N will be used in a feedback 
system (see Fig. 1.2) , The output of N will always be 
non-gaussian (since only saturating nonlinearities are 
used) and in this structure the output is fed back so that 
the input to N is not in general gauss ian. A linear 
system with gaussian input produces a gaussian output, 
Also, it has been empirically observed that when a 
non-gaussian signal is passed through a linear bandpass 
filter the output tends to a gaussian distribution (this 
can be proved for non-gaussian white noise or for vanish- 
ing bandwidth) . Further, the RIDF has exhibited good 
agreement with experimental measurements obtained from 
a feedback structure using two nonlinearities of the type 
considered here, viz., an ideal relay and a relay with 
dead zone. Therefore we shall assume that the noise input 
to N is reasonably well approximated by a gaussian dis- 
tribution . 

In the development of subsequent synthesis theory the 
nonlinear element will be represented by the quasi- linear 



approximations on the right of Eq„ 1,7 and (1.8) will be 
taken as a specification to insure tlie validity of this 
approximation. Eq. 1.13 will not be taken as a specifi- 
cation because it is not possible to simultaneously 
achieve this and the feedback specifications that are 
assumed in the following section. Rather, for any com- 
pleted design the maximum tolerable noise input to the 
nonlinearity can be calculated and the describing func- 
tion holds as long as this limit is .tespected, 

1 . 3 Statement of the Design Problem 

The system structure of Fig. 1.2 will be assvimed 
throughout. P is the plant transfer function which may 
contain varying or uncertain parameters. The designer 
has access to its input and output only. The remaining 
transfer functions are compensations to be synthesized, 
with the exception of N, which is a nonlinear element of 
the type discussed in the preceding section. A sensor 
with unity transfer function which introduces noise with 
power spectrum n is assumed. D(s) is a disturbance 
input and is the applied excitation, present only in 

the case of an externally excited oscillating system 
(EEAS) . 
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Figure 1.2 Oscillating system two degree-of-freedom 
feedback structure. 

The feedback system in Fig. 1,2 is a two degree-of- 
freedom structure, so called because there are two 
signals that can be independently processed and used to 
control the plant output, viz., the command input R(s) 
and the system output C(s). For linear systems, all two 
degree-of-freedom structures are equivalent with respect 
to the attainable benefits of feedback (sensitivity to 
plant parameters and disturbance attenuation) , and with 
respect to the effect of sensor noise n on the 

[Ig] 

plant. Since we replace N by a quasi- linear 

approximator, all other two degree-of-freedom structures 
are equivalent to Fig. 1.2 in the above respect provided 
two constraints imposed by the nonlinear element are 
observed. First, the nonlinear element must be in the 
forward path in series with the plant to provide the zero 
gain sensitivity property. Second, the nonlinear element 
must be embedded between two compensations to allow for 
simultaneously satisfying the quasi-linearity constraints 




and a limit on output oscillation magnitude » 

The oscillating signal component at the input of H 
in Fig, 1,2 will be taken as sinusoidal so the usual 
filter hypothesis of describing function theory must be 
invoked. In particular the output from N will be a 
periodic signal with fundamental frequency and 

Gj^G 2 P m.ust filter out higher harmonics well enough that 
they may be neglected at the input to N . 

The synthesis problem is to design the compensa- 
tion functions in Fig. 1.2 so that the system output 
always satisfies given specifications in spite of bounded 
ignorance or variation of plant parameters and distur- 
bance inputs. Some a priori knowledge of command and 
disturbance inputs must of course be assumed. 

Below we give a statement of the specifications, in 
their most general form, that the subsequent synthesis 
theory is capable of dealing with. In any particular 
design example one or more of these may be absent. 

Specifications 
Nonlinear Rl^nien^ 

The nonlinear element N is taken to be adequately 
characterized by the DIDF 

Wq = (1,15a) 

Nf = M^/A (1.15b) 

when quasi-linearity constraints 
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Uf ( t) 1 < A/oi ( 1 , 1 6a) 

% - (1.16b) 

are satisfied. Some specific nonlinear element must of 
course be selected before a numerical design can be 
executed, so that M and ft ar-f^ lrnr>wn r>-r r-an Vio 

' -- Q, j., , — 

determined by the designer, is forced signal 

component of input to N arising from command or dis- 
turbance inputs. Wj^ is the bandwidth of X^(jw). A is 
the magnitude of the oscillating signal input to N. 

Eq. 1,16 must hold for all possible and A as the 

plant and input varyf so extremes will be treated in syn- 
thesis. As noted previously there may be a lower bound 
on A, e.g., a saturation nonlinearity must be driven 
well into saturation before (1.15) can hold. Also, it 
was noted in Sect. 1.2 from [10] that a > 3, 3 >3 gives 
5% accuracy for the approximation to obtain (1.15) for 
some nonlinearities when x^ is a sinusoid with fre- 
quency in the closed interval Since x^ con- 

sidered here will be a transient signal, 3 may depend 
to some extent upon its waveshape. Ultimately the choice 
of a and 3 is left to the designer, who must consider 
a particular nonlinearity and the nature of the forced 
signals in his system. It will be seen in design examples 
presented later that in practical systems where the SOAS 
or EEAS is most appropriate, that 


is forced by 



other design constraints to be rather large. 


Plant. 

In the subsequent synthesis theory the plant is 
assumed to be a rational transform with known structure, 
written as 

{ w E W 

( 1 . 17 ) 

K(w) e [K^,K2] • 

w is a vector of all plant parameters bounded by the 
set W and K(w) is the high frequency gain factor, 
i.e., P(s,w) ^ K(w)/s® as |s| ^ oo. K(w) 

and K2 = K{\}) . The plant parameters may have any 

value within the bounding set, or may vary slowly within 
this set. By "slowJ.y” it is meant that the change in any 
parameter over the time period required to process an 
input command or disturbance is negligible. The feed- 
back may be quite effective for fast parameter variations 
[191 

also, but since time-invariant synthesis theory is 

used herein the results are valid for slow variations, 
and fast variations are not considered. 

Extreme ^nput 

The extreme command or disturbance input that the 
system is required to accept must be specified. By 
extreme, we mean the forcing signal input that most nearly 
causes the system to violate the quasi-linearity 
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Constraints, This input denoted by or Dp^(s) 

will be the largest magnitude and/or the fastest input 
to be applied to the system. For example, if the command 
inputs are steps, then Rg(s) = q^/s where repre- 

sents the largest step. 

Transfer Function and Parameter Sensitivity 

The transfer function for commands in Fig. 1.2 is 
denoted by T(s) = C^(s)/R(s) where C^(s) is the forced 
signal component of output. We assume that a nominal 
transfer function is specified along with bounds on 

permissible magnitude variation Aln|T(jw) | as a function 
of w. The nominal response and sensitivity may ini- 
tially be given in the time domain as and an 

envelope of permissible responses, e.g., g^^Ct) £ c^(t) £ 
g 2 ^t), where Cjp(t) is the response due to a specific 
command input, say a step. Such specifications may be 
translated approximately to the frequency domain as 
discussed in [22] . 

Although the oscillating systems treated herein 
may be used for nonminimum phase designs, there is 
not available at present a complete synthesis theory to 
satisfy sensitivity specifications of the above nature, 
i.e,, regarding time response. For this reason v/e 
assume T is minimum phase and since T{j(jj) in this 
case is completely determined by |T(jw) | we treat 

only the magnitude in the design procedure for satisfying 



sensitivity specifications. However, if one is presented 
with specifications on both magnitude and angle varia- 
tion for T(jw) , these may be handled with the synthesis 
method presented, 

.Disturbance Response 

In this paper the disturbance D(s) of Fig. 1.2 is 
treated as a deterministic input (a random input with 
sufficient restrictions could be treated using the RIDF) . 
Specifications are assumed given as an upper bound on the 
disturbance transmission magnitude |Tj^(jw)| = 

I Cjj ( joj) /D ( jo3) [. For minimum phase systems such a bound 
can be obtained approximately from a time domain speci- 
fication of the form |c^(t)| £ provided the form 

f22l 

of the disturbance inputs is known. 

When command transmission sensitivity dominates in 
a design, but nevertheless disturbance inputs exist, it 
is usually desirable to place a damping constraint on the 
disturbance response. Even thouah may not be 

specified over an entire frequency band, we require that 
it have no significantly underdamped poles. Since 
Tj^ - 1/(1+L^) where = N^Gj^G 2 P is the forced signal 

loop transmission in Fig. 1.2, we require 

jw) 

1+L^( jw) 

This places an approximate constraint on the complex pole 
pair nearest to the jw axis in 1/(1+L^) , and thus 




0 ), 


(1.18) 
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3-OOsely 15-m.its sett3J.ng time and overshoot of c.(t). 

a 

Also, one may verify on the Nichols chart that (1,18) 
establishes minimum gain and phase margins. 

Outpu t_0 s ci 1 la tion__Cons t r aint 

For a useful oscillating control system there must 
be some constraint on the magnitude of the oscillating 
signal component of signal at some point in the system. 
In Fig. 1.2 we have a sinusoid of magnitude A at the 
input to N, For this sinusoid N is represented by 
gain (Eq. 1.15a), so the output from N has 

magnitude l^^i' Thus specifying that the system output 
have magnitude less than m, the constraint is written 

where I K (jm ) I = )!. 

‘ m hm -* o ' weW' o ' 

Osciy.at^on Fr^uency Variabipn, 

In the SOAS the frecmencv of oscillation w 

o 

occurs at the frequency of 180° phase lag in the loop 

transmission. This will vary some due to inevitable 

inaccuracies in the loop compensation elements. In the 

EEAS we must have a generator that cannot be perfect. 

Thus, let w' be the center or averaae value of w . 

o ■ o 

The variation due to the above sources is parameterized 
Yj, where 




( 1 . 20 ) 
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is assumed known as a specification in the synthesis 
procedure. Note that vary even more due to 

plant parameter variations, but this component is dealt 
with separately in the design, 

.Limit Cycle Qpe^ching 9.^ii^_Margin 

Let from Fig. 1.2 be the open-loop 

transmission for the oscillating signal component. In 
the EEAS it will be necessary to quench the natural limit 
cycle of the system. As will be shown later, this 
introduces a gain margin requirement on L^, written 

I 1 P ' 

where satisfies - -180®. p will depend 

on the nonlinear element that is used. This specifica- 
tion does not apply to the SOAS. 

Aside from the optimizing criterion, which is 
treated in the next section, the above completes a state- 
ment of design specifications and constraints. Admittedly 
some of the above are rather broad in nature and some 
approximations are required. However, they represent a 
reasonable compromise between the desire to be as analy- 
tical as possible on the one hand, and on the other, to 
get results that have practical application. 



1 , 4 Optimization C.riterion 
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The specifications of the preceding section will be 
taken as inviolate, but upon quasi-linearization of the 
system in Fig. 1.2 these can be satisfied by an unlimited 
number of designs provided the plant is minimum phase. 
Therefore it is desirable to use this reservoir of de- 
sign freedom to optimize in some fashion. In particular, 
the rms noise power at the plant input (or output) due to 
sensor noise input n will be minimized, as shown below. 

As explained in Sect. 1.2 the RIDF, = M^/A, for 

odd, static, and single-valued nonlinearities with 
gaussian inputs is a pure gain. Assuming N in this 
class, the sensor noise r\ in Fig. 1.2 sees an equiva- 
lent linear loop transmission This is 

related to oscillating signal loop 

transmission, by the gain factor N^/N^. (from Sect. 1.2 
No/^r == ^ some common nonlinearities) so we 

can equivalently discuss optimization of which is 

more closely related to the specifications. Letting e 
denote the excess of poles over zeros in L^, a number 
which of course remains fixed while comparing the quality 
of any two designs, then at high frequencies 

l^o^3w)I as w (1.22) 

to 

The optimum L^ will be taken as the one with minimum 

^ 00 * This criterion has been used previously for linear 
[22,23] 


systems . 
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In Fig, 1.2 the noise transmission to the plant 
input is 

"L 


T - 

n 


P(1+LJ 


(1.23) 


Using = QP we see that 


|T^(jw) ! = 


Ift^I ; > 1 

|Q(jw) I ; |L (jo)) I < 1 . 


(1.24) 


V 


n is the power density spectrum of the sensor noise 
input, so writing n(o)) = ( jw) n" (- jw) and = 

(-ju) , the mean square noise at the plant input is 

2 1 _ 

^i ^ Tn l ^ ^ * (1.25) 

with the aid of Parseval's formula^^^^, this may also be 
written 


a? = 
1 



-1 r + 


[Tl\fdt 


(1.26) 


Let us tem^porarily take t) as white gaussian noise. 

In a nontrivial oscillating system design problem it 
will be seen subsequently that the oscillating frequency 
“o constrained to fall above some lower limit. For 
the SOAS the magnitude of is constrained by 

(L^(j(*^o^ = -1) the choice of A similar constraint 

is later shown for the EEAS arising from |L^(jto^) | < p. 



Approximately over the frequency range (0,o)^) the 
designer has little control over noise transmission T 

n 

in (1.24) . First, because at low frequencies the sensi- 
tivity and disturbance attenuation specifications 
generally dictate some loop transmission magnitude 
ib^{jw) I > 1 and the designer of course has no control 
over P. Above this frequency range but still below 

cannot be reduced significantly because of 
the magnitude constraint at is only beyond w 

that can be reduced with a steep slope and the 

accompanying large phase lag. In Fig. 1.3a we show a 
hypothetical frequency response for P. With the 
oscillating frequency at co^ as shown, the resultant 
oscillating loop transmission is shown along with 

^o^a'^^r" '^hese are in decibles vs. log u. The noise 
power integral (1.25) is on an arithmetic scale so the 
mean square noise power at the plant input is propor- 
tional to the area under the |q | curve on the arith- 
me tic scale of Fig. 1.3b. Choosing a smaller oscillating 
frequency, say w^/2, the new and Lj^ are shown 

and noise power is now proportional to the area under the 
dashed |Qj^| curve of Fig. 1.3b. In the frequency region 
where il<j^(ja))| < 1, but |L^( jm) | > IP(jw)[, the noise 
is amplified by the loop compensation Q (Eg. 1.24), 
and from Fig. 1.3a |Q(ju))( becomes larger over a larger 

frequency span as is increased. Clearly, we should 

design for the smallest possible, and because noise 
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power is computed on an arithmetic scales, even a frac- 
tion of an octave improvement in may give a signi- 

ficant improvement. 

Now suppose the noise is not white, but has finite 
bandwidth. If the noise power is concentrated in the 
frequency range above then the design is still 

optimized by obtaining minimum 03^, and hence in 

Eq, 1.22, This is often the case in control systems 
since they generally have small bandwidths and it is 
relatively easy to obtain sensors that are free of very 
lov7 frequency noise. It will be seen in Chapter III 
that is minimzed by making |L^{j(o) | as small as 

specifications permit over the low frequency range 
(w < (jJq) . Thus, even if the input noise spectrum is con- 
centrated in a neighborhood of there is probably 

little improvement possible over the above optimization 
scheme . 

The above optimization criterion was based on the 
nonlinearity having a pure gain RIDF, however, it is 
reasonable to expect that it may also produce an optimum 
design, at least for the case of white noise, for a 
broader class of nonlinearities. 

In summary, as each design specification is con- 
sidered the synthesis theory will be developed with the 
goal of minimizing oscillating frequency for 

the EEAS) which in turn yields minimum and rms noise 

power at the plant input. 





CK:?\PTSR II 

SOAS SYNTHESIS FOR SIMULTANEOUSLY SATISFYING 
QUASI -LINEARITY AND OUTPUT LIMIT CYCLE CONSTRAINTS 

2 . 1 Synthesis of an Elementary SPAS . 

In this section the synthesis of a simple SOAS is 
described „ The SOAS adaptive property of zero sensitiv- 
ity to plant gain changes is shown as well as the basic 
cost in terms of loop transmission bandwidth required to 
obtain this property. It is found that three basic fac- 
tors are the determinants of the loop transmission band- 
width, and that they all influence the bandwidth with 
equal weight. These three factors are the spread of 
plant gain variation, the maximum input that the system 
is required to process, and the maximum limit cycle 
magnitude permitted at the system output. 

The structure shov/n in Fig. 2.1 is used with N a 
nonlinear element characterized by the DIDF components 

No = M^/A (2.1a) 

Nf = M^/A (2.1b) 

when quasi-linearity constraints 

"*^^|x^(t)[ 1 Aj^/a (2,2a) 


% 1 “o/^ 


(2.2b) 
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are satisfied, M , a, and 3 are parameters of 

O JL 

N, is the minimum limit cycle magnitude input to N, 

and x_(t) is the forced signal input at N having 

IT 

bandwidth 


R 


F G^/R 

Q =2> O- 



X 



-1 
< — 


G2H 

— o 


p 


c 



Figure 2.1 SOAS structure. 

In Fig. 2,1 P is the given plant which the 
designer cannot alter, and containing uncertainty in its 
gain factor only, say P = KPj^ = K/[s(s+a)] with the 
gain factor uncertainty bounded by K e [K^,K 2 ] and a 
known. The remaining elements in the structure are com- 
pensations to be obtained in the synthesis procedure. 

It is assumed that no disturbances are acting on the 
sys tern . 

The transfer function is specified, and for illus- 

. 2 2 
tratxon suppose it is T(s) = co /(s+ui ) . The limit 

a. a 

cycle component appearing in C, the system output, is 
to have magnitude less than m, and let the extreme 
input be a step R (s) = q /s. 

G T 




It is nenf^ssary that the describing function 
be valid for the significant frequency components that 
must appear in the system output ^ hence the oscillating 
frequency is constrained to be above the transfer func- 
tion bandwidth by the ratio P of Eq. 2.2b. |T(jWj^) | = 

1//2 gives ~ o)^(/2“l)^ so 
m a 

bi > 3w = (/2~1) , (2.3) 

o — m a 

The loop transmissions for oscillating and forced 
signal components are respectively 


M 

L = N G,G.P = ~ G^G^KP, 
O ol2 A 12h 


(2.4a) 


and 

= N^GjG^P = ^ G^G2 KPj^ . (2.4b) 

It has a valid linear transmission interpretation at 
o 

the oscillating frequency only and is interpreted 

as a linear transmission over (0,w^/S) . To sustain a 
limit cycle in the system we must have L^(jw^) - -1, 

which yields from (2.4a) 

1= |M^Gj(ja.^)G2(j«.^)P^(jai^)| . (2.5) 

The right side of (2.5) contains no uncertainty or vari- 
ation so A/K must be constant and therefore can be 
written A/K = A^^/K^. Putting this in (2.4b), 
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M^K 

< 2 - 6 ) 

*1 

which shows that the forced signal loop transmission 
does not change with plant gain, i.e., for any change 
in plant gain K an exactly compensating change occurs 
in 1/A of N^. Since x = FL^/d+L^) there is no 
uncertainty related to the transfer function in this 
design with only plant gain variation. 

Similarly, the only uncertainty in the forcing 
input to N, arises from changes in the command 

input R. Using the extreme input and the specified 
transfer function, we write with the aid of Fig. 2.1, 

X = ,2 7, 

® N^G^P MjKGjP^ 


We are assuming H - 1 throughout this section. The 

design is carried out to insure that X satisfies 

e 

quasi-linearity constraints (2.2) using R , then for 

e 

all other R, also satisfies these constraints. 

Introducing the maximum magnitude limitation m on 
output limit cycle gives, from Fig. 2.1, 


l’^o^2^^%^^2^h^^"o^ I I"' 


(2.8) 


Solving (2,7) for G 2 and substituting in (2.8), 


X^XjcOo) 


M K 

Vl"' 


(2.9) 
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The last inequality is entirely equivalent to (2,8), 
but written in this form it provides some helpful in- 
sight for choosing X^. The right side of (2,9) is com- 
pletely given by specifications. The left side is 
determined entirely by the choice of X , but this con- 

G 

straint on X^ is partly in the time domain and partly 

in the frequency domain. This is seen by using equality 

in (2.2a) to write |x (juj ) I /A, = 

I Xg ( jto^) I / (t) I ] , is constrained by 

^o ^ ^^m above the transfer function bandv/idth so 

[ ( j (i)) T ( jw) I will be a decreasing function of w in 

the range of permissible Since the objective is to 

minimize to^ , we want to maximize |x (jco ) I /A, at 
o e " o ' 1 

fixed ( 1 )^, subject to the quasi-linearity constraints 
(2.2) on X . Observe that using R = q /s, the 
right side of (2.9) contains the factor g^K 2 /K^m, thus 
for any fixed choice of X^ the oscillating frequency 
increases or decreases with this factor. These 
three specifications all influence the required loop 
transmission bandwidth in the same fashion. This will 
continue to be the case even in the most complex SOAS 
designs . 

The selection of x^ is discussed in some detail 
in Appendix A, and it is seen there that if |x^(jw) j is 
required to be a relatively smooth function of w, then 
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X = 


•*x b 


s+w. 


b 


( 2 . 10 ) 


is a. good choic© that also has th© advantage of simplic- 
ity. Of course will have some constraints imposed 

by the particular plant structure and input form. 

Eq. 2 o 10 is for Type 1 plants (one pole at the origin) 
and step inputs. Adopting (2,10) for gives 

t hence The deter- 
mination of and cOj^ will be discussed presently. 

Using X^ from (2.10) and as above in Eq. 2.7, 

the compensation G 2 is now 

-^l^gT a(s+03j^) R^T 

^2 " M.K.X P, " — * (2.11) 

f 1 e h f 1 h 


To obtain the constraint on due to (2.8), G 2 

from (2.11) is substituted giving 

l”o®2<3“o'Ph'3“o> I 

Since we must have ojj^ £ center quantity in 

(2.12) is approximately independent of therefore 

any < to^/8 will not influence the o)^ satisfying 

(2.12) . It is reasonable to give x^(t) a response time 
comparable to that of the system output. Thus take 

b 

or 0 )^. If we take it at the pole 


near w 
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of cance.l.R one po3.e of T in (Eq„ 2,13,) 

and thus simpl.ifies . 

Using the given specific fom for in (2.12), 




2 

.aq. a) 
'■r a 

Vi 


3“o<3V“a' - "'2 


The quantities of (2.13) are sketched in Fig. 2.2 and 
the minimum due to this constraint is labeled. The 

larger 



Figure 2.2 Constraint on due to 3-imitation 

on output limit cycle. 

of the minimum that satisfies (2.13) 

must be taken as the design value of ^o* can further 

approximate in (2.13) using to < m to get 

a o ^ 
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which again shows the basic importance of the quantity 
q^K 2 /Kim. An approximation equivalent to (2.14) is given 
in [6]. One may observe in Fig. 2.2 that addition of 
more poles to T, which then appear as poles of 
(see Eq. 2,11),: have a very beneficial effect on re- 
ducing , 

Since no disturbances have been assumed for the 
present design and there is no uncertainty in , 

there is no gain— bandwidth requirement for loop trans- 
mission L-. Thus the only requirement for G, is to 
^ X 

shape to establish the limit cycle at the proper 

frequency. Finally, after G^^ is selected, the pre- 
filter F may be obtained directly from the relation 
F = T(l+L^)/L^. The design procedure is next demon- 
strated by a numerical example. 

Example 

Speci fications 

(2.15a) 
(2.15b) 

3 (2.16a) 

3 (2.16b) 


(2.17) 


Nonlxnearity: 


Ideal Relay 
with output 
level M. 


No = M^/A = 4M/^A 
= M^/A = 2M/irA 


Quasi-linearity| max|x^(t)| £ A^/a; a> 
Constraints J g > 


Transfer function: 
T(s) = - 


( 6 )^( 15 )^ 


(s+1) (s+6) (s+15) 
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Plant; 


P(s) = KPj^(s) = 


s (s+1) 


•; K e = 1 : 1 , 100 ] (2 = 18) 


Maximum inputs 


R^(s) = q^/s = 10/s (2.19) 

Maximum limit cycle output: 


m = 0.1 


( 2 , 20 ) 


Design 

The constraint on due to the required band- 

width of T is “q i ~ 3(t/2'-l). Tcilcing Wj^ = 1, 
at one of the poles of T, and using in the form 

of Eq, 2.10 yields 


X 


e 


<7 r, q 

X b _ ^x 

S+(0, stl * 
b 


( 2 . 21 ) 


Then 


MfG2 


3g giving from (2.7) 

A.R T 
1 e 


30 (6) ^ (15) ^ 
(s+6) ^(s+15) ^ 


M G^/2 
o 2' 


( 2 . 22 ) 


Applying the output limit cycle constraint. 
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60(6)^{15)^ 


s(s+l) (s+6) ^(s+15) ^ 


< oOOl = rpHL. ( 2 . 23 ) 

^2 




The quantities in (2.23) are plotted in Fig. 2.3 and 
the minimum oscillating frequency is found to be 

^ 26.5 rps. To see the beneficial effect of extra 
poles in T, note that if the two poles at -15 were 
deleted the minimum oscillating frequency vjould be 
37 rps as shown by the dashed segment on Fig. 2.3. 

In the absence of disturbance inputs and parameter 
uncertainty in there is no gain-bandv?idth require- 

ment for L^. The only specific requirement is that 
LfCju^) ~ ~ ^ ~ choose the minimum 

” 26.5 rps and design accordingly, making |L^(jtjo)| 
quite small for w < as shown in Fig. 2.3. The 

resultant rational loop transmission is 


L 


f 


^o/2 

7.585x10^^ (s+1) 

s [ (s^+2 ( .6) 35s+35^) (s^+2 { . 4) 60s+60^) ] ^ ’ 


(2,24) 


Thus 


G 


1 


Vf 


NfG^P 


K,M.G-P, 
1 r z n 


Aj^3.12xl0^.[ (stl) (s+6) (s+15) 


[ (s^+2 ( .6) 35s+35^) (s^+2.6 (.4) 60s+60^) 



(2.25) 




Figure 2.3 


Frequency response of functions 
applicable to design example. 


’.O 
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Only the prefilter remains to complete the design. 
Factoring the numerator of 1 + the data in Table 2.1 

is obtained. The prefilter is then obtained 


Poles 

1 

ft " 

Zeros 

P 


0) 

z 


0) 


,164 

26.8 

1 




. 813 

48. 





.255 

60 





.468 

67 




K 

, o 

= 1 

K 

<X 

= 7.585x10^^ 


Table 2.1* 

Data for L^/(l+L^) for of Eq. 2.24 


from 


F = T(l+L^)/L^ . (2.26) 

Inspection of the data in Table 2.1 for L^/(l+L^) and 
in Eq. 2.17 for T shows that F has 8 zeros and 7 
poles which is not realizable. Additional far poles 

must be added to T and therefore to X in Eq. 2.7. 
5F 

For any rational function F(s), K and K are 

O 

defined by F(s) ^ as [s| ^ 0 and F(s) ^ 

as I s [ *> 00 * 
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From Pig. 2,1 we can write (H =1) ^ 

so must always have the same excess of poles as 

^ for step R) . If such poles are added to 
X^ at frequencies well above they will have negli- 

gible effect on ^^^|x^(t) | and they will not alter the 
calculation for minimum oscillating frequency in Eq. 2.23. 
Thus we add complex pole pairs to T and X^ at 60 and 
67 rps respectively which cancel the zeros of 1 + at 

those locations. The resultant prefilter has parameters 
as given in Table 2.2. Since H is 


Poles 

1 Zeros 

P 

1 

(0 

z 

c: 

60 

1 

1 

1 

6 

6 

15 

15 




.164 

,813 

26.8 

48, 

K = 1 K = 4.9x10"^ 

o 


Table 2.2. Parameters of 
prefilter F. 


being taken as unity this completes the design of all the 
compensation blocks in Fig. 2.1 except for choosing the 
two constants A^, and the relay output level M. With 
an ideal relay the design theory leaves these constants 
completely arbitrary. For some other nonlinearities, 
e.g., saturation, there will be a lower bound for A^^ 
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given in terms of the parameters of These constants 

JL 

will in practice be chosen from consideration of signal 
levels within a given system, q was not given a value 
either, however, we have used = a = 3 and taken 

= 1, so is determined by the selection of A^. 


9 9 




_s_ ^ ^ ^ ^ X 1 


IjOW0x j-iiyj 


the SPAS Oscillation Frequency . 

In the previous section the smallest possible 
oscillation frequency was obtained under the assumption 
that X was to be in some sense a smooth function of (o. 
It was found that the most significant factor requiring 
to be large was the combination of input magnitude 
q^, range of plant gain variation K 2 /K^, and the demand 
for a small output limit cycle component m, so that the 
basic determinant of was the magnitude of the 

quantity q^K 2 /K^m. Recall the following equations 
(Eqs, 2.7-9) . 


X = 
e 


KiMfG2Ph 


(2,27) 




(2.28) 


A^ M^K^m e o o 


(2.29) 


Eqs.- 2.27 and 28 are obtained from Fig. 2.1 with 


H = 1, 
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a,nd (2.29) is obtained froin substitution of Gp from 

(2,27) into (2,28). Eq. 2,29 shows that if one could 

place large magnitude peaking in X (jw) at to , and 

o o 

still maintain control of the time domain quasi-linearity 
constraints, a significantly smaller to^ may satisfy this 
inequality. Recall the quasi-linearity constraints are 

~ (2.30a) 

% ~ “o/^ • (2.30b) 

To provide for peaking in (x^(jw)l a pair of under- 
damped poles are introduced at in such a way that 

’''^t^lxe(t) I can be simultaneously controlled ^ . x (t) 
is written 

-to, t -C (0 t / ^ 

“ ‘^jjtOj^ (x-l)e ® cos [w^t/l-?^] . (2.31) 

In the frequency domain 

= X^^(s) + x^j(s) 

2 ~~ 2 
s+w. s +2Cto s+w 

b o o 


H X^^(s)/H(s) . 


(2,32) 
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SolvincT for H gives 


H = 


1 


D+« 
o o 


x\ 2 [ (^+1) [w^+ (X-1) ] 


s + 


. (2.33) 


With the above choice of X , (t)| = Xq to. can 

be taken as close as desired to q^to^ by taking X 

sufficiently close to unity, while r, is independently 

adjusted to make |l/H(jw^)| and thus Ix {jto )| as 

large as desired. This means that with the ideal X in 

e 

(2.32), (2.29) can be satisfied for any to . Therefore. 

the bound on output limit cycle component (2.28) does not 

constrain oo . 

o 

We prefer to treat the peaking compensation function 
H as a separate function so in Eq, 2.27 is replaced 

by G 2 H giving 


X 


X = 
e 


el 


A^R T 
1 e 


H 


K^M^G-KP, 
1 f 2 h 


(2.34) 


and 


X 


el 


KiMfG2Ph 


(2.35) 


The function in (2^33) is the compensation H shown 
at the output of N in Figo 2*1 so the constraint on 
output limit cycle from (2*27) is now replaced by 


(2.36) 
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It is undesirable to have the highly underdamped zeros 
of H appearing in the loop transmissions and 

so 1/H is inserted at the nonlinearity input as shox-jn 
by Fig, 2.1, The following example demonstrates a 
design using the ideal limit cycle filter H. 

Example 

Specifications 

The specifications are the same as those in the 
example of Sect, 2.1, Eqs. 2.15 through 20. 

Design 

The constraint on oscillating frequency due to 
required transfer function bandwidth is, as in the 
example of Sect. 2.1, — ^^m. ~ 3(t/T-l) 

the low frequency dominant part of 

is taken as 


X 


el 


_ ^x 
s+w, s+1 


(2.37) 


with (jj, = 1, This constrains w > = 3 due to 

D o — b 

quasi-linearity constraint (2.30b). With the narrow limit 

cycle filter R it is possible to design for any w >3. 

o — 

Let us use t>) = 10, (w =3 could be used, but it is 

o o 

seen later that even 10 is not practical) and choose 
X = 3/2. Then from (2.30a and 31) A, = Xaq u, = 4.5q , 

X X O X 

since a = 3 from the specifications. Using X^^^ 

(2.37) in (2.35) gives 


from 
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Mr - 1 e _ 45(6) (15) 

-^f '■'2 2 2 

^ ^ (s+6)'^(s+15)^ 

= M^G^/2 , (2„38) 

The output limit cycle magnitude constraint is 

MolG 2 (ju)o)H(ja)^)Ph(jWo^ I 

90 (6)^(15)^ 
s(s+l) (s+6)^(s+15)^ 

We find at = 10, | G 2 ( jlO)Pj^ ( jlO) | = -15.7db, so 

that |H(jlO)| = “60 + 15.7 = “44,3db. H, given by 
Eq. 2.33 is completely determined except for Solving 

for ? that provides the required |H(jlO) | gives 
C = 0.0015 and 

H(s) = I 

The quantities on both sides of Eq. 2.39 are plotted on 
Fig. 2.4. 

Shaping G^^ to establish the oscillating 
frequency at 10 rps, a satisfactory compensation is 

A. 1 . 26 x 10 '* t (s+l) (s+6) (s+15) ] ^ 

G = 5 5 — , (2.41) 

[ (s +2(.6)13s+13 ) (s‘^+2(.4)25s+25 )] 


s^+2( .0015) 10s +10^ 


s^+2(.0219) 8.17s+8.17^ 


(2.40) 




and the corresponding loop transmission is 



2-20 


Fig. 2 


)0 


isa 


1-4 Q 


4 Frequency response of functions 


applicable to design example . 


iSi 

U1 


Phase Lag (deg) 
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^ 4.58x10^(3+1) . (2.42) 

s [ (s^+2 ( . 6) 13s+13^) (s^+2 ( . 4) 25s+25^] ^ 

There is no magnitude or bandwidth requirement on 
other than that L^(jo)^) = L^(jlO) = - 1/2, so 
jL^(juj) i iias been chosen rather small, A plot of 
is shown on Fig. 2.4, Factoring 1 + gives the data 

on Table 2.3. 


Poles 

i Zeros 

P 


(0 

z 


w 


.172 

10.1 

1 




.79 7 

18,2 





,278 

24.4 





.455 

27.3 


i 1 

1 

K = 
o 

= 1 


K = 

00 

4.58x10^ 


Table 2.3. Parameters of L^/(l+L^) 
using from Eq„ 2,42. 


The prefilter is found from 

F = T(l+L^)/L^ , (2.43) 

and comparing data for T (Eq. 2.17) and (1+L^)/L^ 
from Table 2.3, it is found as in the previous example 
that F has an excess of zeros over poles, and is 
therefore not realizable. With the justification given 
in the design example of Sect. 2.1, we add additional 
far pole pairs to T and at 24.4 and 27.3 rps with 
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proper damping to cancel the zeros of 1 + at 

these locations. These are well above and so have 



Poles 


i_- 

Zeros 

1 

p 

C 

( 1 ) 

r 

1 

1 

! Z 



1 

1 

1 

6 

6 

15 

15 


i 

1 

1 

j 

\ 

1 

.172 

.497 

1 

1 

1 1 

10,1 

18,2 

1 ^0 = 1 


« 

0 

(1 

8 

24 


Table 2.4 Parameters of 
prefilter F. 


negligible effect on the design. The design is now 
complete, except for a choice of q and the relay 
output level M, neither of which are important to 
what is demonstrated in this example. Note that H in 
Eq. 2.40 is unrealizable because it has an equal number 
of poles and zeros. In the structure of Fig. 2.1, G^^/H 

and G 2 H may be realized as single compensation func- 
tions so it is not necessary to make H realizable by 
itself . 

Consider again the idealized 

Xg(t) = x^^(t) + 

M r 4 ./ 2 ^ 

= cos[a)^t/l-i;^] . 


(2.44) 




The quasi- linearity constraint w > Bw, implies tb.at 
should have no significant frequency components 
above but is at frequency 

In fact will be a very lightly damped sinusoid 

at this frequency. In Sect, 1.2 it was observed that 
for odd, static, and single-valued nonlinearities the 
incremental input describing function (IIDF) , Eq. 1.9, 
holds for small sinusoidal components of any frequency 
not rationally related to oj^. Also it was noted that 
the IIDF was the describing function being used for 

the total signal x (t) . Hence, if X in (2.44) is taken 
near unity, one should expect to satisfactorily 

describe the nonlinearity transmission for both x (t) 

0X 

and x^ 2 ^^) • the following section we will not be 

able to make arbitrarily small, so a simulation 

will be necessary to verify the validity of for the 

total Xg (t) • 

We have shown in the above example that the SOAS 
oscillating frequency can be reduced to the limit 
required by quasi-linearity, i.e., = max[BWj^, Boj^] , by 

using narrow and high peaking in = X^^/H. Reducing 

“o general reduce the sensor noise power that 

is transmitted to the plant input. As discussed in 
Sect. 1.4 this noise power is approximately proportional 
to the area under the | ( jm) /P ( jco) | ^ = 

Nf 1 ^1 ( j(u) G 2 ( joi) j ^ curve on an arithmetic scale. The 
quantity { jco) G 2 ( jto) | in db has been plotted 
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in Fig, 2,3 and 2,4 to illustrate its significant re- 
duction as a result of the decreased oscillating fre- 
quency , However, the ideal design scheme in this 
section is not practical. Consider, in Fig. 2.4, the 
effect of a small variation of say from 10 to 

9 rps . The system output limit cycle magnitude is given 
by K 2 M^| G 2 ( jw^) H ( j(o^) ( jw^) I , and at 9 rps this 

quantity has increased by 60 db or a factor of 10^. Thus, 
even small variations in which will always be 

present due to imperfections of the compensating elements, 
can cause the limit cycle component of system output 
to far exceed its specified bound. 

In Sect. 2.1 a rather simple synthesis scheme gave 
a large value of basically determined by the ratio 

present section an idealization in the 
design of has essentially eliminated this constraint 

and allows (o^ to be determined by other, and frequently 
less severe, constraints. The idealization has very 
limited, if any, practical use, but it indicates that 
one can possibly do considerably better than in Sect. 2.1. 
In the following section we shall explore the middle 
ground between these two extremes, and in the process 
develop a synthesis scheme that can be applied to real 
systems , 
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2 c 3 Satisfying Quasi-Linearity and Output. Limit Cycle 
Constraints with Minimum Oscillating Frequency 
and Practical Specifications . 

In Sect. 2.1 it was shown that the combination of 
SOAS output limit cycle magnitude limit, range of plant 
gain variation, and command input magnitude places a 
severe limitation on the minimxjm achievable loop trans- 
mission bandwidth (or oscillation frequency) . This 
phenomenon is accentuated when the plant gain variation 
is large, which is just the situation where the SOAS is 
most useful. In Sect. 2.2, by using a very special kind 
of limit cycle filter the above constraint on bandwidth 
of loop transmission was removed completely, however, it 
was pointed out that this limit cycle filter was too 
much idealized to work in practice. Nevertheless, this 
ideal filter has provided the insight and motivation for 
the design procedure to be developed in this section. In 
the following development the anticipated range of 
oscillation frequency variation is taken as a specifica- 
tion and a notch limit cycle filter is used to reduce 
the oscillation frequency as much as is practical in a 
given design situation. It is found that the cost of 
reducing the oscillation frequency is the necessity to 

increasingly complex filter and its inverse 
in the system. 

The same feedback structure is used, and it is shown 
again in Fig. 2.5 for easy reference. The constrained 
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Fig. 2.5 SOAS feedback structure. 

part of the system, the plant P, is now assumed to 
have uncertainty associated with its dynamics as well as 
the gain factor, written 

{ w e W 

(2.45) 

K(w) e [Kj^,K 2 ] . 

V7 is the vector of all plant parameters and K (w) is 
the high frequency gain factor (P(s,w) K{w)/s® as 

K^(w) = K(w) and In 

subsequent notation the functional dependence on w will 
be omitted unless specifically needed for clarity. The 
oscillating frequency to^ is assumed to vary about a 
center value with the variation quantified by the 

given parameter y, i.e., 

“o ^ t^^o/YrW^Y] = ^ (2.46) 

is not known at this point, rather it emerges from 
the synthesis . Y = Y^Yj^ is made up of a component of 
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specif 5-cation to compejasation uncertainty from 

Eq. 1.20, and a component from plant parameter uncer- 
tainty treated in Chapter III. 

In Fig. 2.5 let be the forced component of 

plant output signal due to either a command or distur- 
bance input and is the corresponding input to N. 

Then = Z^/[N^G 2 HKPj^] = AZ^/ [MjG 2 HKPj^] . There is 
some more restrictive set of plant parameters C W 
and input that produce the extreme nonlinearity input 
X^, i.e,, the fastest and/or largest X^ which comes 
closest to violating quasi- linearity constraints 


max 

t 



(t) I £ A^/a 


(2.46a) 


% - ' 


(2.46b) 


All the extreme parameters and functions are so denoted 
with a sub-e and we have 




A Z 
0 G 

^e “ MlGlHiTp: 


*f'"2““e"he 


(2.47) 


For the remainder of this section the discussion is re- 
stricted to command inputs so is made up of the 

extreme input R and an extreme transfer function T 

^ © 

which is the transfer function only when is at the 

value Pj^^, Thus, 


A R T 

V — e e e 
e M^G..HK P, 


JL j-vj j: m 

f 2 e he 


(2.48) 
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There is a distinct pairing of T with P, and 

Appendix B discusses how is obtained and paired 

with T^. Here we assume the pairing has been made 

and both functions are known. In general may be 

independent of some elements of the parameter vector w 

so that when it is fixed at P. , K can still have 

ne e 

uncertainty e tKel'^'e2^ C [K^,K 2 ] . However, 

in (2.48) is constant by the same arguments used earlier 
when P^ had no uncertainty at all. Thus any value in 
the set ^^el'^e2^ used in the following syn- 

thesis calculations. We shall simply denote this choice 


by K and when A is obtained it is the A that 
e 0 

exists when P = To clarify the above consider 

the following plant. P(s) = K^a/[s(s+a)] where 
w = with and a varying independently in 

and r giving K e [Kj^,K23 = 

* Then ^he ” ( s ( s t a^ ) } whi le 

Kg e ^^el'^e2^ ~ ^^'tl^e'^Jl2®e^ ' might choose 


K 

e 




to use in the synthesis calculations. 


It is again convenient to separate into the 

smooth'* part and the part due to the complex limit cycle 


filter H as 



with 


(2.50) 
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^el chosen in accord 

with Appendix A and assumes step inputs and one pole at 
the origin in the plant. In Sects. 2.1 and 2,2 it was 
found that cannot be realized with a single pole, 

so allows for the addition of far poles. To 

emphasize that the absolute value of A^ or q are 
not determined by the design procedure, the notation 

A^ H (2.51) 

is introduced . 

Summarizing the above we have 


X . fel ^ 

* H M.K G,HP^ 
r e 2 he 


(2.52) 


where 


X 


el 




M.K G»P. 
f e 2 he 


s+to, X ' 

D 


(2.53) 


and X^ must satisfy quasi-linearity constraints 


max 

t 


Ug(t) 


< A /a = Xq 0 ), 
— e ^x b 


(2.54a) 


% 1 • 


Simultaneously the output limit cycle constraint. 


(2.55) 
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rmist be satisfied:, where Ik P, (jw )| = )!. 

m hm o ' weW* o ' 

Solving (2.53) for Gp and substituting in (2.55) gives, 
after some rearrangement. 




Fe<J“o> 


K 


m 


-Km 

e 




M. 




IP 


hm 


<3“o> 


he 




(2.56) 


The task now is to find a suitable combination of 
center oscillating frequency f "^he filter H, and 

~ satisfy (2.56) . Observe that x depends 

on H because each H will produce a different 
^^^[x^(t)j for Eq. 2.54a, and therefore demand a 
different x. Given values for oj' and X, the corre- 
spending |H(jw)| required is obtained from (2.55) as 


H( jw) 


/ 



K m 




e 

f 



^ Xaq w, K 
D m 

M 

O 



X^^(jo,) 

R^n^TT^Tjc^ 




1 



( jeo) 

TW ' 


(li e ^ 


; W ^ fi, 
(2.57) 


The ideal limit cycle filter of (2.57) is shown in 
Fig. 2.6a using a straight line approximation over Ji. 

An obvious modification modification of (2.57) will give 
H with finite slopes shown dashed in Fig. 2.6a. Sub- 
stituting (2.57) into the first relation of (2.52) 
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loaio 


(b) 


Fig. 2.6 Ideal limit cycle filter H 
and corresponding X^. 


yields 



Using from Eq. 2.53 combined with jH(jw)| shown 

in Fig. 2. 6a, the plot of Fig. 2.6b for |x^(jt)i))| results. 
Observe that (2.57) is independent of q , 


while 



57 


q is just a scale factor in (2,58). Similarly, since 
|H(jw) [ is evaluated at ^ in (2.57) where 

coj^/w, the value is closely independent of 
0 )^, i.e., the choice of has an insignificant influ- 

ence on the determination of minimum w , Therefore it 

o 

is reasonable to take near the maximum transfer 

function bandwidth, a: o) . 

D in 

is minimum phase, so using the Hilbert 

r?7l 

transform as applied to this problem by Bode •' the 

phase can be obtained from the magnitude characteristic. 

With both phase and magnitude the time response may then 

be calculated as x (t) =X”^[X (s)]. The following 

scheme is used to obtain satisfactory combinations of 

and H(jw). First choose a value for w' . Then 
o o 

iterate on X to find a solution of 




I = Ag^io^ (2.59) 

For each X pair H(jto) is given by (2.57). In this 

manner one can obtain a plot of X vs. u)' as shown in 

o 

Fig. 2.7. A good estimate for a starting value of is 
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obtained by letting H = 1, Then from Eq„ 2.53 
(assuming negligible contribution to the time response 
X (t) from the far poles in $ ) x = and 

M^X^^(jw)Pj^^(jw) I as shown in Fig. 2.8a. Where this 
quantity equals m/K^ (Eq. 2.55) gives a suitable value 
for w'/y* This is the solution found in Section 2,1 



Fig. 2.8 Behavior of functions of interest 
as 0 )^ is reduced. 
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generalized to permit variation in w . From here we 

o 

work to J.ower frecmencies with w". Each reduction of 

o 

^o H to become a deeper and more complex 

filter and it causes I increase, thereby 

increasing the gain of = ^“qx%^e'^e^^ Ve^elS^he^ 

at all frequencies. The behavior of H, and X 

2 G 

with decreasing is illustrated in Fig. 2.8. 

Assume now that Fig. 2.7 has been constructed in the 
manner explained above. How does the designer use this 
curve? First there is some lower bound dictated for the 
oscillating frequency by the bandwidth of the transfer 
function, ^ Also, a minimum oscillating fre- 

quency is dictated by disturbance response and/or 
parameter sensitivity considerations which are treated 
in the next chapter. Let the largest of these lower 
bounds be shown in Fig, 2.7, such that some freedom 

remains to reduce to^ by taking H / 1. Any 
“o then be chosen as the center value of 

oscillation frequency. The designer must choose a value 
by making an acceptable tradeoff between loop transmission 
bandwidth and the coinplexity of the limit cycle filter H , 
which must be realized twice in the loop of Fig. 2.5, 
Having made such a tradeoff, and A are available 

from Fig, 2.7. The magnitude characteristic for H{j(o) 
is given by (2.57), or its finite slope replacement, from 
which a rational approximation is obtained. The compensa- 
tion function G^ is obtained from (2.52) as 


D e 


60 



M^K X tP. 
f e el he 


( 2 , 60 ) 


The compensation is then shaped to provide the level 

of loop transmission required for disturbance attenuation 


.and/or parameter sensitivity over the low and inter- 

mediate frequency ranges, and to provide 
M 

^oe ~ h~ ^e^l^2^he 180“ of phase lag at the chosen 

e 

oscillating frequency w'. Finally, 

“f"e 

^fe A_ ®l°2'’he ' and the prefilter of Fig. 2,5 can 
e 

be designed from. 


P = T^a+I-£e)/Lfg . (2.61) 

We mention that for a significant class of design 
problems some simplifications are appropriate. Recall 
that was defined (Eq. 2,45) so that it has vanishing 

uncertainty at high frequencies, i.e., Pj^(s) 1/s® for 
large |s[. It folloxvs that when the specifications are 
such that a high oscillating frequency is required, 

*'‘ml^hm^^“o^ I \ ^ ^2' result can be 

used to simplify Eqs. 2.56, 57 and 58. Also, some far 
poles above may be assigned coincident in and 

giving a reduction of order of G 2 seen in 
Eq. 2.60. Further, if these poles are a reasonable 
distance above neglecting them will have no signifi- 

cant effect on the calculation of [H(jw)| or [x^(j(jij)| 
from Eqs. 2.57 and 58 respectively, nor will they effect 



61 


x^(t) significantly. This means their location need not 
be known during the preparation of the curve of Fig, 2,7; 
all we must know is that they lie well above the oscilla- 
ting frequency, say above the resulting when H = 1 

in Fig. 2.7 and 2.8a. Therefore they need not be assigned 
at any specific frequency until w' is determined and 


^4= 

fe 

designed. 

Since Lj- (jw') = 

fe'-* o 

180°, there are 

in- 

evitable poles 

of above 

fe o 

which result in 

far 

off 

zeros of 

1 + L_ . Once 

fe fe 

is known these 

zeros 

of 

1 + can be assigned as far poles in T^ 

and 


Xg which by (2,60) reduces the order of G 2 C and by 
(2.61) reduces the order of F. This is precisely the 
technique used in the examples of Sects. 2.1 and 2,2, 

The synthesis technique of this section is next 
demonstrated with a numerical example. 

2.4 Design Example 

In this section a design example with significant 
plant ignorance is worked out. The objective here is to 
assume that the minimum oscillating frequency is deter- 
mined by the combined constraints of output limit cycle 
magnitude and quasi-linearity, eind then select the 
oscillating frequency and design and H using the 

methods of Sect. 2.3. This portion of the design is later 
incorporated into a complete SOAS design and simulation 
in Chapter V. 
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Specif loat.^ons 

Nonlinearity; 

Ideal relay) N = M /A = 4 m/ttA (2„62a) 

with output) 

level M j = M^/A = 2M/7TA (2.62b) 


Quasi-linearity'^ 

/ 


"*t^ I ( t) I Ag/a= ; a> 3 

% - P i. 3 


Tran.sfer function; 


(2.63a) 

(2.63b) 


T^(s) 


(6)^(25)^(50)^ 

(s+6)^(s+25)^(s+50)^ 


(2,64) 


Plants 


P(s) = KPj^(s) 


K 

s (s+a) 


a e [aj ,a2l 
K e 


11,103 

(2.65) 

[1,103 


The plant parameters a and K are independent. 
Extreme command input; 


Kg(s) = qj./s = 1/s (2.66) 

Maximum limit cycle output; 


m = 0.1 (2.67) 

Limit cycle frequency variation; 

The limit cycle may vary ±12% about the nominal 
, i . 6 . , 

“o ^ t‘*’i/T,wyY3 = Qj Y = 1.12. (2.68) 
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The first step is to determine what plant parameter 
set constitutes the extreme plant This is con- 

sidered in Appendix B and in accord with the conclusions 
there we select 


he 


s (3+32) 


1 

s (s+10) 


(2.69) 


We also note that for the plant given in Eg. 2,65, 

obtained from 


K P, (s) = 
m hm 


^2 

s (s+a^^) 


10 

s(stl) 


(2.70) 


without regard to the particular that emerges from 

the design. 

The low frequency dominant part of is chosen as 


X 


el 


^x% 

s+w^ 


6q 

^x 

s+6 * 


(2.71) 


and the far poles that will eventually appear in X^^ 

are neglected for the present. If no limit cycle filter 

is used, i.e., H = 1, then x (t) = x^, (t) and 

Eqs . 2.63a and 71 give "*^^|x^(t)| = 6q^ = A^/a, while 

= 1 is satisfactory. Therefore 


MfG2 


. aq to, R T 
- X b e e 

K X ,Pi, 
e el he 

= 3(6) ^ ( 25) ^ (SO ) ^ (s+TO ) 
(s+6) (s+25)^(s+50)^ 


= M^G2/2 . 


(2.72) 
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This is the design value for only if H = 1, which 

permits X = 1 in Eq. 2,63a. Otherwj.se Gp in (2.72) 
will be mu.ltiplied by X > 1. Using Gp from above and 
^hm (2.70) 


M 




6(6)^(25)^(50)^(s+10 ) 






< .01 = m/K, 


m 


(2.73) 


The functi-ons on each side of inequality (2.73) are 

plotted on Fiq. 2.9, from which the minimum w* is read 

o 

as = 48.6 y = 54.5 rps . A smaller oscillating fre- 
quency is desired if possible so the limit cycle filter 
design steps are carried out. From Eq. 2.57 


/ 


K^m 

Mf 




Phe<J“> 

( V“Jx“b 

M 

o 


jw)Te(jw) 




w e 


V 


1 


; w ^ Q 


4 


2.77x10 10s (s+1) (s+6) (s+25) (s+50) 

X 6(25)^(50)^(s+10) 




w e 

(u fi. 


(2.74) 


Then 



Magnitude (db/100) 



Fig. 2.9 Frequency responses for functions 
determining lower bound for 


c?> 
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'x^(jw) 


Cf 


X 


K Aao), . 
m b 

M 

o 

jw)T^( jw) [ 


; w e 



'■he '3“)' 

|X3i(jw)|/qj^ 


? w ^ 


/ 


3.6AX10 


= < 


(6) ^(25) ^ (50) ^ (s+10) 


10s (s+1) (s+6)^(s+25)^(s+5Q)^ 


DOJ 


s+6 




; w e 

1 

; w ^ 

(2.75) 


The next task is to use Eq. 2,75 to calculate w' \ 

o 

pairs to produce the X vSo plot* This is implemented 
by getting the phase of using the Bode method^^^^ , 

followed by inversion to time response using the Discrete 
East Fourier Transform (DFFT) . Several points 

are plotted on Fig. 2.10 and connected by a continuous 
curve . 

From consideration of the rapid rise in A below 
30 rps in Fig. 2.10 and the complexity of the required 
filter H, 0)^ = 2 8.2 rps is chosen as a center value of 
oscillating frequency that appears achievable. The ideal 
filter with 500 db/dec slopes centered at = 28.2 rps 
is shown in Fig. 2.11. A rational approximation to the 
ideal filter is also shown in Fig. 2.11 and has the 
P^J^3.meters tabulated on Table 2.5. Also the plot of 
^o^^2 ^ been added to Fig. 2.9. 






















Fig. 2,11 Ideal and rational 

limit cycle filters. 


Shown in Fig, 2.12 are time responses of the 
original (t) from Eq. 2.71 and the resultant x^(t) 

after modification by H. Of course some far poles 
must eventually be added to X^(s) when designing the 
final loop transmission and these will set x^(0) =0. 
With \ = 1.22 we get A^/q^ = = 22 which is very 

close to the peak value of x^(t) shown in the figure. 

The design value of ^^^2 now obtained by 

multiplying the quantity on the right of (2.72) by 
X = 1.22, completing the portion of the design to be 
demonstrated in this section. Use of the limit cycle 
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Fig. 2.12. Time response ax^j^(t)/q^ 

and ax (t)/q using rational 
limit cycle filter from Table 2.5. 

filter has reduced the oscillating frequency by the 
factor 54.5/28.2 = 1.93, or nearly an octave. 

The filter in Fig. 2.11 is not quite as deep as it 
should be at some frequencies. The maximum error in the 
interval is about 2.5 db^which will manifest itself 

in the completed system design by producing an output 
limit cycle in excess of the specification by the factor 
1.34. Of course this is only a violation when the plant 
has maximum gain = 10.. Note the clear visibility of 
the tradeoffs among various specifications. The com- 
promises available to avoid redesigning H are a) 
allow the larger limit cycle, i.e., set m = .134, 
b) demand smaller reduce lG 2 (jw)| (Eq. 2.73) 
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to bring the limit cycle in tolerance,, c) reduce 
|G 2 {jw) 1 and require smaller inputs, j< 1/1.3 = .75, 
d) reduce [G 2 (jw) | and attempt to get by with a small 
violation of quasi- linearity constraint {2.63b) when the 
maximum step input is present. We shall consider the re- 
mainder of the system design for this example in 
Chapter V. 



CHAPTRH III 


SATISFYING SOAS SENSITIVITY SPECIFICATIONS 
WITH UNCERTAINTY IN PLANT DYNAMICS 

3 . 1 The Approach 

The constrained part of the system, the plant, in 
Fig. 3.1 is described by transfer function 

j w e W 

P(s,w) = K(w)P (s,w) ? / (3.1) 

[k{w) e [K^,K2], 

(P(s,w) ->■ K(w)/s^ as |s| ^ <») with the parameter 
bounding sets W, [Kj^,K 2 l , and the plant structure known. 
The SOAS property gives zero sensitivity to variations in 


R-& 


F 


G^/H 

■o- . * o- 



G,H 
— =5^ 




Fig. 3.1 SOAS structure. 

the high frequency gain factor K(w) . In this chapter a 
synthesis method is given for satisfying fixed specifi- 
cations on sensitivity to plant dynamics with a minimum 
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bandwidth loop transmission, i,e«, a loop transmission 
that is optimal in the sense of Sect, 1,4, 

The forced signal loop transmission is 

M^K 

Lf = = ^ 1 ^ 2 ’^h ' ( 2 . 2 ) 

A 

and the transmission for command inputs is 

T = FL^/(1+L^) , (3.3) 

It is assumed that specifications are given as 
bounds of acceptable variation Aln[T(jw) ] as a function 
of 0 ). A method for satisfying such specifications in 
linear systems has been given by Horowitz and Sidi^^^^, 
and this method will be adapted to the present problem. 

The prefilter F is a fixed transmission so Aln[T(jw)] = 
Aln[F(ju))L^(jo))/(l+L^{jw))] = Aln|L^(jw)/(l+L^( jco) ) I + 
j Ayiif ( joj) / (1+L^ ( jco) ) . Similarly, if we assume A/K does 
not change with then Am [L^ ( jw) ] = Ain [Pj^ ( jw) ] 

since the remaining fixed parts of cancel in the 

difference. Thus the specifications and plant uncertainty 
are both in terms of L^(jw), viz, Aln[L^( jo))/l+L^(ja)) ] 
and Aln[L^(j(o)] respectively. The two quantities are 
conveniently related on the Nichols chart which has 
coordinates of ln[L^(jw)| vs. ( jw) . For a given 

range of plant uncertainty Aln[Pj^(jw)] = Aln[L^ ( jw) ] , 
the specifications Ain[T(j(o)] = Aln[L^{ jto) / (1+L^ ( jw) ) ] 
may be translated into bounds on acceptable L^(jw). The 
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details of this translation are given in Appendix C. A 
loop transmission is then synthesized which has 

minimum gain-bandwidth while satisfying the sensitivity 
bounds o 

3 . 2 The Equivalent Plant for Forced Signals 

To apply the method of [22] to the SOAS it is 
necessary to write the loop transmission as = GP^, 
where contains all the parameter ignorance of the 

loop. P^ is the equivalent linear plant seen by forced 
signals in the SOAS. Prom the nonlinearity describing 
functions 


= M^/A (3.4a) 

= M^/A, (3.4b) 

the forced and oscillating loop transmission are 


and 




(3.5) 


L, = 


NfGiGjKPh 




N M 

fr '■o = »T ^ 

o o 


(3.6) 


To sustain the limit cycle L^(joj^) = -l, and using this 
in (3.6) gives 


1= I ♦ 


(3.7) 
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If only the plant gain factor changes, A/K in (3.7) 
remains constant. However, when the parameters of 
change then [Pj^(ju))[ will be changed as well as 
y^h ^ phase change in in turn causes the 

oscillating frequency to shift and the net result 

is that the right side of (3.7) changes in a rather 
complicated way with the varying parameters of Pj^. To 
obtain an expression for the equivalent plant, A/K from 
(3.7) is substituted back in of (3.6) yielding 


= 


^f^l^2^h 


M 




— G 


^f^h 




= GPf ^ 

where we have defined 


(3.8) 


and 


G — 


^ 1^2 


(3.9a) 


Pf = K^Pj^ = M^Pj^/|M^G(ja)^)Pj^(j03^) I . (3.9b) 

The quantity = Mj^/[m^G( jm^)Pj^( jw^) [ is a varying 
gain factor that is part of the equivalent plant P^. For 
example, suppose = 1/2, G = l/(s+10), and 

Pjj = l/[s(s+a)] with a e [1,4]. Then = 

M^/A[s(s+a) (s+10) ] . For a = 1, the oscillating 
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frequency is determined by = -180*^ which 

gives uij = 3.16. Solving for |lj^(j3.16)| = 1/2 
gives = 55, Thus = 55/[s (s+1) (s+10) ] and 

this function is plotted in Fig. 3,2. Now letting a = 4, 



Fig. 3.2 Illustration of gain factor 

variation with parameters of Pj^. 

we get ~ (s+1) / (s+4) which is shown dashed in 

Fig. 3.2. ( j = -180 gives the new oscillating 

frequency at CO 2 = 6.32, but at this frequency 
I^f2^^^2^ I is -20 db so A must change value so that 
Lf 2 (jw 2 ^ = -1/2. Solving, M^/A^ = 280^ and L ^2 = 

280/ [s (s+4) (s+10) ] . As a varies from 1 to 4, A varies 
by the ratio A^/A^ = 280/55 = 5.1. This gain variation 
must be expressed as part of the uncertainty of the loop 




transmission and accounted for in the design to satisfy 
transfer function sensitivity specifications. The 
magnitude of the gain uncertainty depends on the compen- 
sation as well as the parameter uncertainty of Pj^. 

In the eljove example ^fl 


ratio of these two loop transmissions gives 


7 \ / 1 \ 

“ 1 ^ 


I G ( jw^) /G ( jw 2 )Pj ^2 ^^^ 2 ^ ^ write the 

equivalent plant P^ as a function of G as given by 
Eq. 3,9b, In summary, we have the loop transmission 


Lf = GP^, (3.10a) 

V7here 

^f = ^f^h “ Vh/I (3.10b) 

which is to be synthesized to satisfy given transfer 
function sensitivity specifications with minimum band- 
width, Due to the constraint L_(i(o ) = -M^/M , 
minimizing loop transmission bandwidth is equivalent to 
obtaining the smallest possible subject to the 

sensitivity specifications. Neither co nor G are 

o 

known at the outset, while the equivalent plant 
depends on both of these quantities. Thus an iteration 
scheme will be implemented whereby we converge to the 
minimum oscillating frequency and the required compensa- 
tion G simultaneously. 
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3.3 Synthesis Df SPAS Loop Transmission 

to Satisfy Sensitivity Specifications 

In this section a method for determining the minimum 
bandwidth (minimum w^) L^(jw) to satisfy given sensitiv- 
ity specifications is developed, Assiame that bounds on 
the acceptable magnitude variation, ln|T(jw)|, are 
given as shown in Fig, 3.3a and b. Also shown in 



Fig. 3.3 Permissible variation of ln|T(jw)| 
vs. variation of ln|Pj^ ( jw) | . 


Fig. 3.3b is a hypothetical set of bounds for plant 
uncertainty Aln| P^^ ( ju) | . In general the loop transmis- 
sion is = GP^ = GK^Pj^ with P^ dependent on G and 
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in the manner given by Eq. 3.10b. However, since 
neither of these are known at this point, we first neg- 
lect the gain factor of and use Pj^(jto) to 

obtain a set of bounds for L^(j(o) using the specifica- 
tions on Aln|T(jw) | and the method detailed in Appendix 
C. Since varies with Pj^, the bounds apply to some 

reference value of plant parameters denoted by L 

f 1 

GPhi. The bounds are conveniently displayed on Nichols 



chart coordinates {ln|L^(jo)) [ vs. /L^(jw) ) as shown in 
Fig. 3.4 labeled ( 1 ) 2 , . . L^^(ju^) must lie on or 

above the bound at in order to satisfy the 
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sensitivit-y specification Aln|T(jw^)[. It is noted in 
Appendix C that these bounds can also include specifica- 
tions on Ay^T(jco^) if desired, however, for minimum 
phase design the magnitude specification is sufficient. 
There are in theory an infinite number of bounds, but in 
practice the designer chooses enough frequencies to 
provide high confidence that L^(jw) ivill satisfy the 
intermediate bounds that have not been computed. At some 
frequency near, or possibly as much as an octave above, 
the Aln|T(jw)|, Aln|P^(jw)l crossing in Fig. 3.3b the 
bounds will cease to exist or become insignificant 
because the permitted variation in ln|T(j(o)| is larger 
than that exhibited by the plant, i.e., no feedback is 
required in this frequency range. Several other bounds 
are shown in Fig. 3.4 whose construction is investigated 
next. 

In the SOAS we permit only one frequency to satisfy 
( j oi) = -180°, which is the oscillating frequency 
Therefore, there is some maximum phase lag boundary, 
in Fig. 3.4, given by ( jto) ^ 6^ > -180° which 

applies over the intermediate frequency range. In the 
next frequency range the bound comes into effect. 

This bound is derived from the disturbance response 
damping constraint 


jw) 


1+L^ ( jm) 




(3.11) 
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The Significance of (3.11) has been discussed in 
Sect. 1.3. At the oscillating frequency is con- 
strained to pass through the point = ~M^/M 

Of course Si in (3.11) must be large enough so that 

B. excludes this point. Above w , where the final 

o 

poles of are introduced, some limit on maximum 

peaking is imposed, |L.(joa) | < -A (db) for u > w . 

i — 0 

This gives the horizontal segment of of Fig. 3,4. 

Finally, must have some specified excess of poles 

over zeros e, so that eventually -90e“, 

giving the vertical segment of B 2 . 

The boundaries at are applicable to 

^fl but they are constructed to guarantee satis- 

^sction of the transfer function sensitivity specifica- 
tions over all possible plant parameters provided the 
Specific value satisfies its bounds (and assuming 

for the present the gain factor variation in is 

negligible) . This allows the designer to work only with 
bfi in full confidence that all other will be 

satisfactory. 

The above is not true for the heavy boundaries Bj^, 
and B^ of Fig. 3.4 which apply to all L^. It is 
possible to design to satisfy these bounds as 

indicated by Fig. 3.5, and then find for another plant, 
S3.y that L ^2 violates the bounds as shown. It 

can be seen from Fig. 3.5 that the situation would not be 
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improved by using as the reference plant to design 

with, A reasonable approach to this problem is to choose 
^hl parameter set that gives maKimum phase lag 



Pig. 3.5 Violation of boundary 

which applies to all L^. 


over the low frequency band. Then must be designed 

to avoid the higher frequency portions of B 2 by 

some margin. The margin can only be found by checking 
the remaining L^(jo)). 

Finally r one additional constraint is placed on the 
oscillating frequency Recall that the describing 

function has a valid frequency interpretation over 

( 0 , 01 ^/ 8 ), where 3 is the quasi-linearity parameter 
from Eq. 1.16b. Therefore, to insure the validity of 
over the frequency range where sensitivity reduction is 
required we require that 


OJ > 

o — 


3w, 


n 


(3.12) 
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for all 

A summary of properties of L^(jw) which will be 
quite close to optj.mal in the gain-bandwidth sense of 
Sect. 1.4 is given belov/. 


Property 1 . 


3. 


5. 


Lj;{jw ) = for all P, . 

f o f' o h 

o) > p(D for all P. . 
o — n h 

L^(jm) does not cross the bounds 
and B 2 for any Pj^. 

, i = l,2,...n, is on or above 
the respective bound at If 

is above the bound at then it has the 

maximum phase lag permitted by Property 3. 
(For example, if L^^(j« 2 ) maximum 

phase lag of all possible L^(ju) 2 ), then 
if above the bound of 

Fig. 3.4, it is on B^^.) 

For w > (0 , L-- (iw) is as close to 

boimds B^ and B 2 as permitted by 
Properties 1. and 3. 


These properties follow closely those given in 
[22,23] for linear design, but are modified to allow for 
the special constraints of the SOAS loop transmission. 
They are intended to serve as a guide to the designer 
for shaping a real rational loop transmission. A 
practical L^^(jw) might appear as shown by Fig. 3.6. 




Fig. 3.6 Hypothetical loop transmission 
satisfying all boundaries. 

Next we consider the gain factor of that has 

been thus far neglected. We have 


(3.13) 

with 

The sensitivity bounds at are first 

obtained by neglecting the gain factor uncertainty in 
and using Pj^ in the manner of Appendix C. Then 
Lf, (jw) = GP , { jw) is designed such that satisfies 
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the sensitivity bounds and all satisfy the bounds 

®d/ ®2“ the design for the value 

of for plant may be obtained, and the com- 

pensation is G = Lfi/P^i* ^^th any other plant value 
the loop transmission is = GK^Pj^ where 

= ^f/ I ]^qG ( jw^) ( j W q) [ is the gain factor required 

to satisfy Lj^(jw^) = For a significant class 

of design problems this gain factor will exhibit very 
small variation, and can therefore be neglected com- 
pletely. The reason that this is true is that w tends 

o 

to be well above the varying plant dynamics and Pj^ is 
defined with vanishing uncertainty at high frequencies 
(see Fig. 3.3b). If it is ascertained that the gain 
uncertainty is negligible, then the design of for 

sensitivity is complete. 

In any particular design problem the gain factor 
variation may not be insignificant, in which case we 
must insure that the sensitivity specifications are met 
with the smallest possible bandwidth loop transmission, 
taking the uncertainty into account. From the 

design of above, a function g is obtained. This G 

can now be used to calculate P^ for any plant parameter 
and with p^ we can obtain a new set of Nichols 
chart bounds for sensitivity at to , < jD 2 f . . . . Then we 

reshape as required to satisfy the new bounds, 

obtaining a new G and repeating. The above procedure 
is summarized in the following design algorithm. 
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Algorithms 

Step 1. Choose a reference plant parameter set de- 
noted by ^hl’ Take and 

approximate P^ ~ p^^ for remaining plant 
parameter values. Use as the plant for 

obtaining bounds on L^^(jo)^}, i = l,2,...n 
in accord with Appendix C. 

2. Synthesize ^fl = GPj^j^ approximating the 
optimal lity properties given earlier in this 
section. 

3. Use G = to calculate P^ in 

Eg. 3.14, and use to obtain new bounds 

for L^{jo)^), i = l,2,...n. Repeat Steps 2 

^rid 3 until the realized and the bounds 

rl 

computed using G from this are 

arbitrarily close. 

If the gain factor M^/|m^G( jw^) Pj^{ | of P^ 

has negligible uncertainty this will be evident the first 
time Step 3 of the algorithm is executed and this portion 
of the design is complete at that point. 

Each time an ^fl is formed it is constrained by 
Property 1 to pass through the value and the 

frequency at this point is the oscillating frequency 
It will be minimized in accord with how well the remaining 
properties are satisfied. 

When a satisfactory design has been completed by the 



methods of this section one may check the extremes of 
oscillating frequency variation and obtain 


and 


0 ) ' 


/CO (0 

V °max 


o 


itlin 


(3cl5) 


V. 


•n 


iii\ / (j) 

i/”o ' o . 

V max min 


( 1 


. 6 ) 


parameterizes the maximum variation in oscillating 
frequency due to parameter uncertainty in This 

number is used in the design procedure of Sect. 2.3 to 
determine the required width of the notch limit cycle 
filter. ( 0 ^ is the minimum center value for oscillating 
frequency permitted by the present design consideration, 
i.e,, transfer function sensitivity to parameter uncer- 
tainty Of . However, other design considerations , 

e.g., those of Chapter II, may dictate a higher oscilla- 
ting frequency. 

3 . 4 Design Example 

A numerical design example is carried out below to 
demonstrate the design techniques of the preceding 
section. 
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£pe c i f i cat ions 


Nonlinearity: 


Ideal relay | 
with output \ 
level M j 


N 


N, 


Quasi-linearity 

Constraint 


} 


= M^/A = 4M/7TA 
= M^/A = 2M/itA 


3 


(3.17a) 

(3.17b) 

(3,18) 


In Eq. 3.18 is the highest frequency at which 

there is a nontrivial sensitivity specification. 

Plant; 


P (s) = i — y a e [a,,a 2 l = U,10] (3.19) 

s (s+a) 

This plant has been deliberately chosen with the 
pole variation not confined to the low frequency range so 
that it is likely to produce some variation in oscillating 
frequency. 

Disturbance Response Damping: 

L^(jw) 

1+L^( jo)) 


5 = 2i\f 


( 1 ) 


( 3 . 20 ) 


Transfer Function Sensitivity: 

The upper and lower bounds for acceptable transfer 
functions are shown in Fig. 3.7a below. The lower bound 
is = [(6) (25) (50)/[(s+6) (s+25) (s+50)]]^. In 

Fig. 3.7b the permissible variation in the transfer 
function is compared with the existing plant variation or 
uncertainty. 
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ro 

o 


max Pj^(joi) [j \ 


: miniTCjw) 


10 w-rps 


Fig. 3.7b Comparison of plant variation 
with permissible variation of 
transfer function, 

De£i£n 

The first step in the design algorithm in Sect, 
is to obtain bounds for L^j^(jw) using the method of 
Appendix C. The forced signal loop transmission is 


GP^ » 


with 




Vh 


To obtain the initial bounds, the uncertainty of 
M^:/ 1 M^G { j (ij )P, (jti) ) I is neglected, and we use P^ - 
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As a reference plant 


P 


hi ~ 


1 

s (s+a^) 


1 

s (s+1) 


(3.23) 


is chosen and the resultant bounds for 
shown as the solid bounds at oa = .4,. 
. 3.3. Also shown xs the boundairy 



GP 


hi 


8/. ..^4. 8 in 

B, due to 
d 


are 


Eq. 3.20. 


A rational function is found as plotted in 

Fig. 3.8. The plotted points correspond to the boundary 
frequencies in sequence, a convention that we shall 
follow throughout. The parameters of the compensation G 
are listed in Table 3.1a and further data on is 


Poles 

Zeros 1 

P 

C 

I CO 

z 

c 

to 


.6 

15 

3 




.4 

60 

18 




.4 

100 





.5 

120 

. 

1 


K 

= 21. 

23 

K = 

4.59x10^^ 1 

o 



00 


1 


Table 3.1a Parameters of G for L, 

f 

shown in Fig. 3.8. 

given in Table 3.1b. has been shown dashed on 

Fig. 3.8 when = l/[s(s+10)]. Of course each 

must pass through the point = -1/2 (-6 db,-180“) . 




Fig. 3.8 Bounds for L^^(j(jo) and plots of 
rational using G of 

Table 3.1a. 
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o 


K 

O 

K 

00 

a 

16 . 1 

1.0 

21 „ 23 

4.59x10^^ 

1 

17.3 

1.15 

8.11 

5.26 

3 

18.4 

1.31 

5.56 

6.01 

5 

19.4 

1.49 

4.51 

6. 82 

7 

20.8 

1.78 

1 3.77 

8.14 

1 

<0^ = 18.4 


Yh = 1.14 



Table 3.1b Parameters of L_ = GP^ 

f f 

with G given in Table 3.1a. 

Obtaining a rational function for to give a 

reasonable fit of the boundaries is a difficult task/ 
which is largely accomplished by a process of cut and try. 
Some help may be obtained by using the average relation- 
ships between magnitude slope and phase for minimum phase 
functions as explained in [22] . The problem here is 
considerably more difficult than in linear system design. 
For example, in Fig. 3.8 the last sensitivity boundary is 
at 4.8 rps v/hile L^^(jl6.1) = -1/2. If this were a 
linear system loop transmission one could add a complex 
pole pair, at say 30 rps, having virtually no effect at 
4.8 rps and below where the sensitivity bounds are 
located. However, for the SOAS loop transmission the same 
pair of poles at 30 rps will have a significant phase 
effect at 16.1 rps, hence shifting the oscillating fre- 
quency. The gain factor is then adjusted so that at the 
new oscillating frequency, say cd^^, L^^(jw^) = -1/2. 
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The resultant effect is that is shifted away from 

all the sensitivity bounds . 

Next we use the gain factor obtained from the 

first compensation G to calculate bounds for the 
plant is listed in Table 3,1b where a 

spread of 1.78 (5 db) is noted. The new bounds are 
shown dashed on Fig. 3.8 for comparison with those 
obtained previously. Note they are considerably lower 
allowing a smaller bandwidth than the first bounds 

yielded. The new bounds are shown solid on Fig. 3.9 
along with the second iteration for Data for this 

is listed on Tables 3.2a and 3.2b. Comparing the 
data of Tables 3.1b with 3.2b, the oscillating frequency 


Poles 

Zeros 

P 


W 

z 


(0 

1.5 

.6 
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.4 

100 
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Table 3.2a Parameters of G used 
for of Fig. 3.9. 




Fig. 3.9 Bounds for L^^{jw) after 
obtaining first G, and 
using second G from Table 3.2a. 






Table 3.3a Parcumeters of G for 

of Fig. 3.10. 
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1 O 


K ^ 

o 

K 

CO 


11.1 

1.00 

10.8 

1,43 10^^ 

-- 

1 

12.2 

1.21 

4.36 

1,73 

3 

13.1 

1.46 

3.15 

2.08 

5 

13.9 

1,74 

1 2.68 

i 

2.48 

7 

15,0 

1 j 

2.21 

1 . 

2.39 

3.16 

10 

I = 

13.4 

^h 

1.17 



Table 3.3b Paraineters of using 

G from Table 3o3a 


with the previous set, and are shovm solid on Fig» 3.10. 
Data for the next is given in Tables 3.3a and 3,3b 

while is sketched on Fig. 3.10. 

The next set of bounds are calculated using 
from Table 3.3b, These are shown dashed on Fig. 3.10, 
where it is observed that they are quite close to the 
previous set, i.e., the algorithm has conveirged to an 
acceptable G. 

Comparing data from Tables 3.1b and 3.3b, the itera- 
tion has reduced oscillating frequency by the ratio 
1®»4/13.4 — 1.37 while lov/ering the high frequency 
asymptote of by 201og[4. 59/1.43] = 10.1 db. This 

xs a significant improvement that should not be ignored. 

It is notable that if the entire plant were P = p 

h 

from Eq, 3,19, i.e., no high frequency gain variation, 
then the loop transmission required in a linear design 
is specified by the bounds in Fig. 3.8, For the SOAS the 




Fig. 3.10 Bounds for L^j^(jo)) after obtaining 
second G, and L^j^(jo)) using 
third G from Table 3.3a. 
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lower bounds of Fiq» 3„10 are sufficient. Hence^ pro- 
vided the remaining SOAG design considerations allow the 
loop transmission to be obtained, the SOAS is superior to 
a linear design in this example even without any plant 
gain factor uncertainty. 

3-5 Additional Considerations in the Parameter 
Sensitivity Design 

The method of synthesis in Sect. 3.3 which has been 
demonstrated in Sect. 3.4 provides the designer with the 
means for finding the minimum mean oscillating frequency 
(0^ permitted by parameter uncertainty, as well as 

the uncertainty which is written 

“o ® 

then be used in determining the width required for 

the limit cycle filter in Sect. 2.3. There are other 

lower bounds for the bound established by the 

spread of gain uncertainty in Chapter II. Thus it may not 

be possible to achieve the dictated by Sect. 3.3. 

When this is the case some freedom is permitted to reduce 

the uncertainty parameter y. , The variation of ui is 

n o 

caused by the phase variation of P, in the interval n, , 

n h 

Recalling that has vanishing uncertainty at high 

frequencies (see Eq* 3«,1 and Fig. 3.3) ^ any increase in 
can be expected to decrease the variation. On 

Tables 3,1b and 3,3b in the numerical example of the pre- 
ceding section it is seen that a change of from 13.4 
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to 18,4 reduces y.^ from 1.17 to 1,14. In addition, for 
a change in plant phase A0 in the neighborhood of 
the corresponding change in is proportional to the 

slope of the loop transmission phase characteristic. 

This phenomenon is illustrated in Fig. 3.11, where we see 
that it is desirable to make the /L^(jw) characteristic 
as steep as possible in this neighborhood. When oo' must 




Fig. 3.11 Illustration of oscillation frequency 
variation proportional to phase slope. 

be significantly above the value found from Sect. 3.3, 
the frequency range between and the required can 

be used to shape ( jco) for to e with a very steep 

slope. This is done by removing some phase lag above 
and then introducing underdamped poles near to^. Fig, 3,12 
qualitatively illustrates the appearance of on 

Nichols chart coordinates when the above is implemented. 
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Fig. 3.12 Shaping of to reduce 

w variation, 
o 


The validity of the disturbance damping boundary 
and boundary 82 ^ is open to question at frequencies above 
Wq/P (3 is the quasi-linearity parameter of Eq. 1.16b) 
where the forced signal describing function becomes 

ambiguous. The conseirvative approach is to avoid these 
boundaries, although it might be shown experimentally that 
this is not entirely necessary. 

Clearly the synthesis procedure of Sect. 3.3 entails 
a considerable amount of design effort, while we have 
noted that other design considerations such as adapting 
to large gain uncertainty with small limit cycle in 
Chapter II may preclude the use of the minimum found 

in Sect. 3,3. Therefore, the designer should ascertain 
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approximately the achievable to handle gain varia- 

tion and output limit cycle requirements before carrying 
out the synthesis of Sect. 3.3 in detail. In many situa- 
tions the tedious design of that section may be avoided 
completely, because the gain uncertainty and output limit 
cycle magnitude considerations obviously demand a larger 
loop transmission. 



CHAPTER IV 


DESIGN CONSIDERATIONS FOR 
DISTURBANCE INPUTS TO THE SOAS 

4.1 Satisfying Quasi-Linearity and Output Limit 
Cycle Constraints for Disturbance Inputs . 

Disturbance inputs are taken as deterministic signals 
D applied to the plant output as depicted by Fig. 4.1. 
Disturbance signals that act internally, or at the plant 
input, can be translated to equivalent signals at the out- 
put. The nonlinear element N is characterized by DIDF 



(4.1a) 

N^ = M^/A, (4.1b) 

when quasi- linearity constraints 

^t^|xe(t)[ < Ag/a (4.2a) 

“b - “o^^ (4.2b) 

are satisfied. Simultaneously, the limit on maximum mag- 
nitude of output limit cyclo component ^ 
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must be satisfied. In this equation Ik P, (iw ) I = 

P ( I ' where the plant ignorance is characterized 

as in Sect, 1.3, Eq. 1.17, Chapter II presented a design 

technique to satisfy (4.2) and (4.3) when x (t) , the 

e 

extreme forced signal input to N, is due to a command 
input R in Fig. 4.1. Here some modifications of that pro- 
cedure are considered when x^(t) is the forced input to N 
due to disturbance D. 

The plant output signal due to a disturbance in Fig. 
4,1 is 


= -DL^/(1+L^) , (4.4) 

where is the forced signal loop transmission 

hf = N^G^G2 KPj^ . (4.5) 

Let us take an analysis viewpoint temporarily. Suppose 
is the forced input at W due to commands and Xj^ is that 
resulting from disturbances. From Fig. 4.1, 


V = TR 
R N^G.,HKP, 


Aj^TR 


f"2 


KiMfG2HPh 


A^R 


KiMfG2HPh 


1 + L. 


F 


(4.6) 


where we have used A/K = A^/K^, and T 


= FL^/(1+L^) 


For a disturbance input -D, we get similarly. 


X. 


D 


A^D 


KiM^G2HPh 


X, 


R RPl 



1 + 


(4.7) 



Suppose has been designed as in Chapter II (for step 
commands and a Type 1 plant) , so that the dominant part 
of is 


X, 


q 


UJ. 

X b 


s+w. 


b 


(4.8) 


•The frequency response, |Xj^(j(o)[, is sketched in Fig. 4.2 
showing as is often the case because "hhe 

limit cycle constraint. Eg. 4.3, forces to be high. If 
the disturbance inputs <?C“*[D(s)] have approximately -the 



Fig. 4.2 Comparison of forced inputs to N 

for commands R and disturbances D. 

same speed and magnitude as the prefilter output signal in 
Fig. 4.1, ^ [R (s) F (s) ] , due to command inputs, then by 

(4,7) the two forced signals at N will be the same, and 
will both satisfy the quasi-linearity constraints. Since 
J^R(t) is slower than required to satisfy (4.2b) (Uj^ < 
10^3), we could speed up Xj^(t) to x^(t) = 
without violating (4.2b). This gives 
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Xd(s) 


s+w^/3 ' 


(4,9) 


which is also sketched in Fig. 4.2a. 


For the system structure being considered here (step 
R and Type 1 plant) F tends to be a low-pass filter. 

Thus, if D(s) and R(s)F(s) are to be comparable signals 
with respect to their propensity to violate the quasi-lin- 
earity constraints on N, then D(s) must have smaller band- 
width than R(s) , In particular, a step D(s) will probably 
be an unacceptable input to an SOAS that has been designed 
on the basis of step R(s) . We noiv return to the synthesis 
problem. 

The nonlinearity forced input is obtained in terms of 
the plant output signal as 




(4.10) 


with 


= RT = RFL^/ (1+L^) (command input), (4,11) 

and 

= -DL^/(1+L^) (disturbance input). (4.12) 

When there is plant uncertainty there is some extreme set 
of plant parameters which produce the extreme 

due to disturbance input -D^ , so that 

In view of Eqs . 4.10 thru 4. 12, the synthesis procedure 
when and X^ are due to a disturbance input is 

exactly the same as explained in Sect. 2.3 when is 
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due to a conimand Input, except t.b,at In all th.e equations 
of that section Z = R T is replaced bv Z = 

DgLfe/ Cl+h^p,) . Thus we can parallel many of the equat- 
ions from Sect, 2.3 in the following. Beginni.ng with 


= Xaq w, , 
e b ' 


(4.14) 


we get 


X = 
e 


el 

H 


A Z 
e e 




laq W, D L- 
^x b e fe 

M^K G-HP, (1+L- ) ' 

f e 2, he fe 


(4.15) 


and 


el 


Aaq to. D L- 
b e fe 

M^K G..P. “ (1+L^ ) 
f e 2 he fe 


b 


(4.16) 


and tOj^ are parameters of due to a disturbance, and 
the chosen form for X^^ on the right of (4.16) is approp- 
riate only if has one pole at the origin. Recall that 
^fe^^“o^ ^ Using this, solving for from 

(4.16), and substituting in (4,3) yields 





A 

e 


lV<3“o> 


K |m I 

- K m|M -M I 

e f o 


hm 


< 3 %) 


(jw ) 


(4.17) 


he o' 

This is the equation corresponding to (2.56) , but now 
applying to X^ generated by a disturbance input. 

Suppose D = gj/s, and no limit cycle filter is used 


e '^d 
max 


so H = 1, and ^ [x^(t) | - giving A^ = 

Then (4.17) reduces approximately to 


aq CO 


- x'^b • 
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1 

s+oT 


,aK q j 
xt! < 


o 




{4.18) 


or approxxmatrng 


iw +a), 

■ o b ' 


we get 


„ e ' f o' 

- “ dK [M 1 
m o ' 


(4.19) 


For step D and H = 1 there will be a value of lo 
e Q 

satisfying (4.17) only if the step magnitude satisfies 
(4.19). In other words, without a limit cycle filter the 
system will only tolerate very small step disturbances 
and continue to operate in a quasi-linear fashion. For 


example, using numbers from the design of Sect. 2.4, 

i.e., = m = .1, a = 3, and |(M^-M^)/M^| = 1, 

yields q^ < .0033. This is smaller than the smallest 

value of limit cycle magnitude, which has peak value 

m/K = .01. 
m 


Continuing with the design equations (which are not 
just valid for step D^, but for any with one pole at 
the origin) , we assume that can vary in 

w e (4.20) 

Solving for |H(jto) | from (4.17), 


H( jw) 




Phe<3“) 


Dg(jto) 



W E 


(4.21) 


^ ; w ^ U . 

Again taking = q^s, we get |H(jo)) | =; K^m | -M^ | / 

[K^Xaqd[Mo|] for w e fl, which is entirely given by 

specifications and is independent of the choice of co . 

o 

least one additional finite pole, then the 
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right sa.de of (4,21) will increase with w and the higher 
the choice of the less depth required in H(jw). In 

any event, the design problem is to find a suitable com- 
bination of X, and HCjw), and this is accomplished 

by iteration on and A in precisely the manner exp- 
lained in Sect . 2.3 when the extreme X is due to a 

e 

command. The only difference is that a new defining 
function for H( jco) , (4.21), is used. If no narrow limit 

cycle filter is used, then the minimum permissible 
oscillating is the smallest frequency satisfying 

(4.17) with H = 1. 


Suppose now that X, and H have been found. How 

is the design completed? Solving Eq. 4.15 for Gp gives 


Xaq «, D 

^x b e fe 


(4.22) 


We cannot determine Gp uniquely until has been 

designed. can be designed in the manner detailed in 

Sects. 3.3 and 4.2 to satisfy both parameter sensitivity 

and/or disturbance attenuation specifications over the 

low and intermediate frequency bands, while constraining 

it to satisfy L^^(jai^) = at the selected 

above. Vilhen is thus obtained, G^ is completely 

determined and G^^ is solved from 

Xao). 

^I'^'^x ^ M_K G„P, * (4.23) 

f e 2 he 

Lastly, the prefilter F is obtained in the same manner as 


in Sect. 2.3 as 
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(4„24) 


This completes the synthesis steps for the case when 
minimum is determined by the need to satisfy quasi” 
linearity constraints (4.2) and output limit cycle con- 
straint (4.3) with a disturbance input. 

We have seen that the constraints of this section 
are difficult to satisfy (if they can be satisfied at 
all) when the disturbance inputs are steps. However, 
this will not always preclude the use of an SOAS when 
step disturbance inputs are present. Usually one is 
more interested in simply attenuating the effect of the 
disturbance at the system output, rather than shaping 
the disturbance response function in any particular 
manner. Therefore, it may, in some cases, be acceptable 
to design the system so that it remains quasi-linear for 
all commands, but operates nonlinear ly for some distur- 
bance inputs . 

4 . 2 Design to Satisfy Disturbance 
Attenuation Specifications 

As suggested in Sect. 1.4, specifications on the 
maximum disturbance transmission magnitude as a function 
of frequency may be given. The synthesis scheme for 
satisfying such specifications is considered below, and 
shown to be a simple extension of the technique used in 
Chapter III for sensitivity specifications. 

We assume the bounded plant ignorance 
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p7 e W 

P(s,w) = K(w)P, (s,V7) ; < (4.25) 

[^K(w) £ [K-,,Kpl, 

with P(s,w) -» K(w)/s© as [sj -> ■», and the bounding 
sets W, [K^,K 2 l» and the plant structure known. As 
derived in Sect. 3.2, the equivalent plant for forced 
signals (command or disturbance) is 

Pf = , (4.26) 

where 

Kf = M^/lM^G(ja)^)Pj^(jw^) | , (4.27) 

when the forced signal loop transmission is written 

From the system structure of Fig. 4.1 the disturbance 
transmission to the system output is 


T 


D 


1 

I+L^ 


(4.29) 


Specifications are taken as an upper bound on ln|Tjj(ju) | 
vs. 10 , which is to be satisfied over the range of par- 
ameter ignorance for the equivalent plant P^. The 
specifications on ln|Tp(jto) | = In 1 1/ (1+L^ ( ju) | is 
easily translated to a set of boundaries on L^^(j(o) = 
GCjfo) P^ 2 ^ ( jfo) when bounds on the ignorance of P^ are 
known. The method for doing so is explained in Appendix 
C. As in the sensitivity design, the ignorance in is 
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not known at the outset so it is neglected in finding the 

first set of bounds for which = P, is used. The net 

f h 

result is a set of Nichols chart bounds for . Once 
such bounds are obtained the design procedure for finding 
an optimal satisfying them is exactly the same as 

given in Chapter III to satisfy sensitivity bounds. In 
fact, if both sensitivity and disturbance attenuation 
specifications are given, the design should be carried 
out to satisfy both simultaneously. At any frequency, 
say w., a boundary, B in Fig. 4.3, for L^^(jw.) is 
obtained by the technique of Appendix C to satisfy the 
sensitivity specification. Similarly the boundary of 



Fig. 4.3 Composite boundary, for 

due to sensitivity and ^ 

disturbance specifications. 

Fig. 4.3 is obtained from the disturbance response 
specification for the same . Then the boundary 

with greater must be respected to satisfy both 



sensitivity and disturbance response specifications. 
Henpe, the oniy alteration of the design technique in 
Chapter III is in the calculation of the bounds for 



CHAPTER V 


EXAMPLE OF SOAS DESIGN AND SIMULATION 

5 . 1 Design Specifications 

Listed below are the necessary specifications for 
a complete SOAS design in the context of this paper. It 
will be evident that the specifications are selected to 
take maximum advantage of the design efforts previously 
invested in Sects. 2.4 and 3.4. 

Structure 

The design will be implemented with the same feed- 
back structure of previous sections, shown again in 
Fig. 5.1. 



Fig. 5.1 SOAS structure. 


Nonlinear element ^saturation) 

The nonlinear element is taken as a saturation with 
parameters of the characteristic defined by Fig. 5.2. 

The DIDF is 



1.14 


N = M /A = 4M/itA 
o o 

= M^/A = 2M/ttA, 

with quasi-linearity constraints 

^f^|x^(t) I < A/a; a ^ 3 




4J 



Input 


(5.1a) 

(5,1b) 

(5.2a) 

(5.2b) 


Fig. 5.2 Saturation characteristic. 


Pl^t 

The plant is specified by 


P (s) 


KPj^(s) 


K 

s (s+a) 


{ a e [aj^,a 2 ] 
K e [K^,K 2 ] 


where a and K are uncorrelated. 


tl,10] 
El, 10) 
(5.3) 


Disturbance response_dam£ing 

To avoid excessive overshoot and ringing due to 
disturbance inputs the specification 


( jw) 
1+L^ ( jw) 


< a = 2; y to 


(5.4) 
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is imposed. This specification was discussed in Sect. 1.3, 
Eq. 1.18. The disturbance inputs are assumed sufficiently 
small that no specific transmission requirements are im- 
posed, and also that such inputs will not intolerably 
.violate the quasi-linearity specifications. 

Transfer_f unction and sen^sitivit^ 

The transfer function magnitude must lie between 
the boxinds given in Sect. 3.4, Fig. 3.7a, repeated here 
as Fig. 5.3 for convenient reference. These bounds 
establish both sensitivity and the absolute level for 
|T(jw) j. 

Extreme conmand in£ut 

The extreme (largest and fastest) command input 
anticipated is 

Re<s^ = qj-/s; kj.! = 1. (5.5) 

The maximxam output limit cycle magnitude over all 
plant conditions is 

m = 0.1. (5.6) 

Note that this is a 'peak' value. 

Oscillation frequency certainty 

We assume in the absence of plant uncertainty the 
oscillation frequency may vary by + 6% due to parameter 



Fig, 5.3 Specification envelope for |T(jo))|. 
variation in the compensating elements. Thus 

“o ^ “ ‘^C-' ■'c " 

when the plant parameters are fixed. 

5 . 2 System Design 

The first step is to determine the minimum os cilia' 
ting frequency allowed by the transfer function sensi- 
tivity specification. This has been done in Sect. 3.4, 
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and the minimiDn frequency found to' be. = 13.4rps with 
variation parameter due to plant uncertainty of 1.17. 

The data for this minimum bandwidth sensitivity design 
is given on Tables 3.3a & b. 

Next the combined constraints of quasi-linearity 
(5.2) and output limit cycle limit (5.6) are considered. 
This part of the problem with the present specifications 
was treated in Sect. 2.4. There it was found that 
0 )^ = 28.2 rps is an achievable oscillating frequency 
with a moderately complex limit cycle filter H, which 
is significantly better than is possible without the 
filter. This value of course takes precedence over the 
smaller value demanded by sensitivity. 

Before proceeding with the design we need to con- 
sider two additional requirements on the oscillating 
frequency. First, recall that co^ must be high enough so 
that the assumed is appropriate over the transfer 

function bandwidth. Fig. 5.3 shows that this bandwidth 
is permitted to reuige over approximately [4,13] rps. 
However, the actual variation will be considerably less 
because of the inherent SOAS property of removing 
nearly all high frequency gain variation (of course at 
high frequency is ambiguous anyway, but the effect 

extends to the intermediate frequency range also) . The 
crucial question is does the transfer function vary from 
4 rps upward, or from 13 rps downward. The former is 
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true in the present problem because in Sect, 2.4 the 
extremes of plant and transfer function, T and P, , 
used for the extreme were the small gain-bandwidth 

cases. Thus = 28 rps is judged to be sufficiently 

high for transfer function bandwidth requirements. 

Lastly, must be large enough to render the 

sensitivity design meaningful. There is little question 
that 28 is sufficient, since the last sensitivity 
specification was at about 3 to 5 rps (Fig, 3.10) . Thus 
we proceed to design for = 28 rps. 

The design value of determined in 

Sect. 2.4, and is obtained by multiplying Eq. 2.72 by 
the appropriate X = 1.22 foxand in that section. The 
parameters are listed on Table 5.1. 
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Zeros 1 
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10 



K = 220/M. K = 2.06x1oVm. 

0 f 00 ' f 


Table 5.1 Parameters of G 2 . 

Now must be obtained such that the Nichols chart 
sensitivity boundaries are satisfied over the low and 




intermediai-e frequency ranges, and so that- the limit 
cycle is estab.ij.shed about the chosen nomina3. = 28 
rps. For this purpose it is easier to work with the 
entire 3.oop compensation G, defined in Sect. 3.2 as 
Gj^G 2 / and the forced signal equivalent plant i.e., 

Lf = G (5.8) 

M.P 

'■f = TTTGTTypJTJ^ = Vh- 

A function for G which satisfies the sensitivity 
bounds, but sets (jj^ too small, is available from 
Table 3.3a. This function is modified, primarily by 
moving the complex pole pair near up in frequency, 

enough to establish the new center oscillating fre- 
quency near 28 rps . This is again a cut and try process 
and the result is tabulated on Tab3-e 5.2, By increasing 
ti)^ significantly the spread of variation is 
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Table 5.2a Parameters of G for 
shown in Fig. 5.4. 
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8.26 

1.30 
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5.14 

1.35 

5 

29.1 

1.11 

3.82 

1.40 

7 

30,0 

1.19 

2.85 

1.49 

10 

1 Wq = 28.4 


= 1.057 


Table 5.2b Parameters of using G 
from Table 5.2a . 


considerably reduced/ as is the spread of variation; 
both may be checked by comparing data from Table 3.3b 
and 5.2b. Recalling the iterative procedure for satis- 
fying the sensitivity bounds, from Sect. 3.3, the new 
G has to be checked and possibly adjusted somewhat to 
insure that the final satisfies the bounds corre- 
sponding to the final G. Here only the completed design 
is given with the data of Table 5.2 and the plots on 
Fig. 5.4. In this figure with 

= l/[s(s+l)]. Frequency responses of for 

several values of P^^ are given in Fig, 5.5. 


Thus 


M 

To extract G, , = GP, = GK. P, = -i-K G, G-P. . 

1 fe fe fe he A e 1 2 he 

e 


G 

1 - 


K^MfGj 


(5.10) 


Recall that = aXq^Uj^ is the value of A when the 

plant is P = K^P, / where K = 

e he e 


1 was used for this 




Fig. 5.4 Bounds for and plot of 

Lfj^(jo)) using G of Table 5.2. 










(db/lOo) 



10”^ • 10® o)-rps 10 ‘ 10^ 10® 

Fig, 5.5 Frequency responses of L^. 
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l?3 


design in Sect., 2,4. is the from Table 5,2b 

that pairs with = l/[s(s+10)] used in Sect, 2,4, 

i,e,, K = 1.19. Also from Sect, 2,4, X = 1.22 
f s 

and = 6. ^ given by Tables 5.1 and 

5.2a respectively. Collecting the data, has the 

parameters tabulated on Table 5.3. For completeness 
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Table 5.3 Parameters of 

the parameters of H are given in Table 5.4. 

The system prefilter F is obtained from 

(5.11) 


{5.12) 


where 

T = (6)^(25)^(50)^ 

® (s+6)^(s+25)^(s+50)^ 
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Table 5.4 Parameters of H. 

This T^ is the lower bound in Fig, 5,3 and the ex- 
treme transfer fimction used in Sect. 2.4. The para- 
meters of L^^/{l+L^^) are listed in Table 5,5. One 
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Table 5.5 Parameters of L_ /(1+L_ ). 

re re 


finds F in (5.11) with 10 zeros and 8 poles, hence 
unrealizable. We assign additional pole pairs to 

and T^ at 98 and 121 rps with damping factors to 
cancel these zeros of 1 + These poles are 

sufficiently far off that they should have no appre- 
ciable effect on the time responses x^(t) or c(t) . 
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Table 5.6 Parameters of F. 

The resultant prefilter is given by Table 5.6. The 

design is now finished, with all functions in Fig. 5.1 

completely identified except for the selection of q 

or . 
e 

Before proceeding to simulation of the completed 
design we can do some verification by using the quasi- 
linear model, i.e., we calculate the linear response 
functions for Fig. 5.1 when N is replaced by the 
gain factor M^/A. The particular A that goes with a 
given plant is accounted for by using the appropriate 
gain factors from Table 5.2b. First the transfer 
functions T(jw) are calculated and the magnitudes are 
shown by Fig. 5.6. In this figure the lower response 
(a = 10) is identically T^ which is also the lower 
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bound in Fig, 5.3. The upper boioind (a = 1) is very 
close to the upper bound of Fig. 5.3 up to about 5 rps 
and then falls beJ.ow the specification bound. Thus our 
goal of barely achieving the sensitivity specifications 
is met and not exceeded except at higher frequencies 
where it is unavoidable (and where T(jw) is ambiguous). 
Finally we verify that the bandwidth ^ ~3db point / lies 
within [3. 5, 7. 5] rps which seems compatible with 
OD^ 28 rps. 

Shown in Fig. 5.7 are the system outputs c(t), 
which will be useful in evaluating simulation results 
later. For a = 1 the overshoot is rather large 
(=i 35%) , but this is permitted by the transfer function 
specifications , 

Next, shown in Fig. 5.8 are several x^(t) . From 
these we can check the validity of the selections made 
for P. and T when designing the limit cycle filter 
in Sect. 2.4, i.e,, does the a priori selection made 
truly yield the largest and/or fastest x^(t). Figure 
5 . 8d permits easy comparison. At small t they are 
identical since for large w in the frequency domain 
the X^(jw) are identical. For larger t the de- 
signed x^(b) (a = 10) consistently has larger magnitude 
peaks than the remaining two signals. It does not seem 
possible to detect any meaningful difference in speed 
of the Xj(t) by simple inspection of Fig. 5.8. 
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Fig. 5.8d Simultaneous plots of x^(t) . 

The conclusion drawn is that the selections of P. , 

he' 

are quite satisfactory. However, the change in 
x^(t) with in the present example is rather small 

so the selection of Pj^^ is probably not a critical 
factor in this design. Note that the proper pairing of 
with Pj^^ is verified by Fig. 5.6. That is, we 
designed the transfer function T^ on the lower speci- 
fication bound so that for values of P, other than 

n 

^he' |T(j(ii)) I should increase, as it does in Fig. 5.6. 
Finally let us compare the resultant in Fig. 5. 8a 

with the design value in Sect. 2.4, Fig. 2.12. The two 
signals are identical after an initial transient in 





Fig, 5,8a which is the effect of the added far poles. 
This transient does exceed the Sect. 2.4 design value 
of (t)!= 22. by about 35%. This narrow pulse 

should not cause any appreciable effect on the output 
c(t) , because even if the nonlinear element does trans- 
mit it, the smoothing effect of the plant integrator 
will essentially remove it from the system output. 

5 . 3 Digital Simulation of Command Response 

The system design of the previous section was 
simulated using the MIMIC digital simulation language 
as adapted for the CDC 6400 computer. 

It is noteworthy that the simulation was first 
attempted with an ideal relay nonlinearity. In propa- 
gating the limit cycle through the ideal relay the zero- 
crossing times must be detected with rather high accur- 
acy, and errors in the detection produce low frequency 
transients, or ideally subharmonics in the nonlinearity 
output signal. The difficulty is highly accentuated 
in the oscillating system because of the extremely high 
low frequency gain compared to the gain at the limit 
cycle frequency w^(see Fig. 2.9) between the non- 
linearity output and the system output. This causes the 
low frequency error terms to be greatly amplified at 
the system output, and these components swamp out the 
limit cycle signal. The net result was that an 
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extremely small digital simulation integration step size 
was required, making the simulation very expensive in 
terns of computer running time. This difficulty was re- 
moved by changing the nonlinear element to a saturation. 
For an ideal relay with output level M the output 
error is 2M when the simulation fails to detect a zero 
crossing. For saturation, the output error will be pro- 
portional to the time increment by which the z.ero-cross- 
ing is missed, and therefore reduces rapidly with inte- 
gration step size, making it possible to execute the 
simulation at reasonable cost. The lesson of this 
experience is that, although the synthesis is independent 
of the particular nonlinearity used, there will indeed 
be practical considerations that give some nonlineari- 
ties advantages over others in a given application. 

Shown below are simulation results using a satura- 
tion nonlinearity with the parameters listed on Fig. 5. 2, 
and using system input 

is a unit step = Qj-Zs = 1/s in all cases. 

Fig. 5.9 shows the output from prefilter F of Fig. 5.1. 
This signal is of course the same for all plant parameter 
cases. Figure 5.10 shows the nonlinearity input signal 
for two parameter cases. The forced signal component 
x^(t) is essentially indistinguishable in these plots. 
The system output c(t) is shown for several plant 
parameter cases in Fig. 5.11. Comparison of these step 



responses with those of Fig, 5.7 shows th.at the non- 
linear system response closely approximates the response 
of the quasi-linear approximation used in synthesis. 

The specified limit cycle magnitude has been scaled on 
several plots of Fig, 5.11 and shows satisfactory agree- 
ment with the observed limit cycle. 



Fig. 5.9 Digital simulation plot of 
prefilter F output. 
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Pig. 5.11a System response to iinit step 
applied at t = 1. sec. 
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5,4 Digital Sitnulatibn of Distuirbance Response 

This section presents some disturbance response sim- 
ulation results obtained from the system design of Sect, 
5.2. Any evaluation of these results must be rather 
qualitative, since the design does not conform to any 
specifications on the disturbance response. Nevertheless, 
they do provide useful information about system perform- 
ance for disturbance inputs. All simulations were per- 
formed with P(s) = K^Pj^^(s) = l/[s(s+10)], because this 
parameter set is the worst case in terms of magnitude 
quasi- linearity constraint (5.2a). 

It is obvious from the treatment in Chapter IV that 
the system will not be quasi- linear in the presence of a 
unit step disturbance. Regardless, a unit step input, 

D = = -1/s in Fig. 5.1, was simulated and the 

resultant system output is shown in Fig. 5.12. In Fig. 
5.12b the disturbance input is applied at one half the 
limit cycle period (ir/w^) later than in Fig. 5.12a. 

The radical difference in the two responses demonstrates 
the well known phenomenon that the oscillating system 
response is time dependent when the quasi-linearity con- 
straints are violated. Both responses exhibit a predom- 
inant time constant x which is believed to be the time 
required to return to quasi-linear operation. Later plots 
show this time constant to be dependent on the magnitude 
of step input. Both responses of Fig, 5.12 appear quite 
undesirable for most applications. 
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Nextf a smaller step, = “»025, was simulated. 

The magnitude of this input was obtained by using the 
design values of ^he' ^hm' other required 

parameters from this chapter in Eq. 4.21. Evaluating 
this equation at = 30 rps (corresponding to ^he from 

Table 5.2b), the given q^ is obtained. This does not 
gaurantee quasi- linear performance however, because the 
present design values of and H have not been 

selected to produce quasi-linear operation for any step 
disturbance input. In Fig. 5.13 the frequency responses 
of the forced input to N are compared for R = 1/s and 
Dg = -.025/s. It is observed that the signal due to 
is rich in frequency components above This will be 

the case irrespective of the step magnitude, i.e., a step 
disturbance input can be made small enough to satisfy the 
magnitude quasi-linearity constraint, but the frequency 
constraint cannot be satisfied. Fig. 5.14 shows simu- 
lated outputs for the input = -.025/s. These exhibit 
time dependence and overshoots an order of magnitude 
greater than the input, hence they are quite undesirable 
responses. 

We know that if the signal f (t) e ^ [R^ (s) F (s) ] 
is applied as a disturbance input, then will be 

the same as for command inputs. The signal f(t) is 
shown in Fig. 5.9. This signal was time scaled to inr 
crease its speed by factors of /2* and 2, and applied as 
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a disturbance input in the digital simulation. The sys- 
tem outputs obtained are shown in Fig. 5.15, and appear 
quite reasonable. 



Fig. 5.12a System output for unit 
step disturbance. 
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••>0 .49 .14 l.t* 1.72 t.lS 2.81 9.4i 9.44 S.IT 

time-sec 



Fig. 5.14a System output for step 
disturbance input. 
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Fig. 5.15a System response to 

disturbance f(/2t). 
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Fig. 5.15b System response to 
disturbance f{2t). 




lOc(t) 
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Fig. 5.15c System response to 

disturbance f(2t)/2. 






CFAPTES VI 


SYNTHESIS OF AN EXTERNALLY 
EXCITED OSCILLATING ADAPTIVE SYSTEM 

6.1 Introduction 

It is possible, with the proper type of nonlinear 
element, to externally excite the oscillating signal 
which circulates in an oscillating feedback system rather 
than using the limit cycle which occurs naturally in 
the SOAS. Such a system, discussed below, will be termed 
an externally excited oscillating adaptive system (EEAS) . 
The two degree-of -freedom feedback structure to be used 
is depicted in Fig. 6.1. A^ is a zero mean periodic 
excitation signal of frequency which we take to be 

sinusoidal. This input need not be sinusoidal provided 
G, G„P is of sufficient low pass character to maintain 

JL ^ 

the validity of the describing function used for N . 

F, G. , and are linear time-invariant compensations 

to be synthesized, and P is the plant or constrained 



Fig. 6.1 EEAS feedback structure 
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part of the system for which parameter ignorance as 
detailed in Sect. 1,3, Eq. 1,17, is assumed. The general- 
ity and limitations of this structure have been pointed 
out in Sect, 1.3, 

System performance specifications and optimization 
criterion for the EEAS are identical to those taken for 
the SOAS and stated in Sects. 1.3 and 1.4 respectively. 

With the exception of one new constraint, to be de- 
veloped in the following section, it will be observed 
that the EEAS synthesis methods developed below in large 
measure parallel those for the SOAS. 

6 . 2 EEAS Nonlinear Element 

At the outset it is required that the nonlinear 
element N in the EEAS possess all the properties re- 
quired for the SOAS. These properties were treated in 
Sect, 1,2. In particular it is necessary that the DIDF 
components 

^o" ^o^'^ (6.1a) 

M^/A (6.1b) 

represent a satisfactory quasi-linear approximator when 

^^^|x^(t) I < A/a, (6.2a) 

and 
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All quantities in the preceding two equations are as 
previously defined, e.g., in Sect. 1,2, Eq. 1.7 and 1.8. 


If is zero or sufficiently small in Fig. 6,1, 

and the loop transmission has in excess of 180® phase lag, 
then the system will self -oscillate. Indeed, with A^ = 0, 
we have the SOAS treated in the preceding chapters. Thus 
if the DIDF in (6.1) is to be used, the system must be 
synthesized under conditions that prevent the occurrence 
of the natural limit cycle, i.e., there must be only one 
sinusoidal input at N. We shall give sufficient condi- 
tions for limit cycle quenching below. 

Assume there are two sinusoidal inputs to N 
giving x(t) = A sin w t + B sin co t, where w is the 
limit cycle frequency and is the frequency of the 

excitation A^. Using the techniques of Sect. 1.2, the 
two sinusoid input describing function (TSIDF) may be 
calculated by repeated application of Eg. 1.4b, and the 
components labeled N^(A,B) and Ng(A,B) . Letting Lg 
denote the loop transmission seen by the B sin w^t com- 
ponent, the limit cycle will be quenched provided 

|Lg(jio^) I < 1. (6.3) 

Suppose that ™g^|Ng(A,B) [ = ^(A) exists and is finite. 
Then the limit cycle quenching condition may be written 

C(A) |G^(j(o^)G2(j(o^)P(jw^)| < 1|\/ P. (6.4) 

In general this condition may be quite complicated to 
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satisfy. A survey of a few odd# static, single-valued 
nonlinearities seems to indicate the following. For 
nonsaturating nonlinearities g(A) does not exist. 
However, these are of no interest because their DIDF’s 
are not of the form required in (6.1). For three satvir- 
atiiiy nonliuearities (ideal relay, relay witii dead zone, 
and saturation) whose TSIDF's are plotted in [14], g(A) 
has the characteristic shown in Fig, 6.2, and asymptotic 
value 5(A) •> i/j/A where ij; is a constant. The de- 
parture from the 


C 

Fig, 6.2 Character of (A,B) | = g(A). 

asymptote occurs when A is so small that the nonlin- 
earity is rapidly changing character, e.g,, saturation 
when A is so small that the element is linear. For 

the ideal relay with output M, 5(A) = \|;/A = .855M/A for 
r 1 cl 

all A > 0 . We restrict the EEAS synthesis theory 



to nonlinear elements such that i|) exists and 5(A)<i|(/A 
for all A > 0 . Using this restriction 



is sufficient to quench the limit cycle and DIDF theory 
may be applied. Defining s p and recalling the 

oscillating signal loop transmission is = N^Gj^G 2 P = 

M 

^ we get 

iLoCjta^) I < P, (6.6a) 

or 

|Lf(j‘^,j)l < pIMj/M^I . (6,6fo) 

Eq. 6.6 gives the limit cycle quenching constraint clear- 
ly in terms of the loop transmission to be synthesized, 
i.e., it is merely a minimum gain margin. This is the 
gain margin mentioned briefly under specifications in 
Sect. 1.3, Eq. 1.21, and it is taken as a specification 
for the EEAS design. It is not surprising that a gain 
margin emerges as the condition for limit cycle quench- 
ing, however, it may be somewhat surprising that P > 1 
is satisfactory for some non line ari ti es . For the ideal 
relay P = = 4/(.855ir) = 1,49. We have verified 

satisfactory limit cycle quenching by analog simulation 
with |L^(jw^)| = 1.26(+2db). 

In all subsequent discussion of the EEAS, Eq. 6.6 
is taken as a specification and it is therefore assumed 
that the nonliearity input signal is composed of a 

forced signal component x^(t) and a single oscillating 
component due to the excitation input A^. 
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6 , 3 EEAS Sensitivity to Plant Gain Changes and Behavior 
of L ^{joj^ ) with Changing Plant DynaMics 
The most attractive property of the SOAS is its 
inherent zero sensitivity to plant gain changes. In this 
section it is shown that the EEAS retains this property. 

If the oscillating component of X at the input to 
N in Fig. 6.1 is denoted by then 


X = -A GtG„P -XL 
o o 1 2 o o 




thus (6.7) becomes 



<j) Tan [a^/M^ +*^0034) “ An (6.11) 

provided is real. If is complex, 4» is replaced 

<}> + and by |m^ 1 in the Tan"^ term of 

(6.11). In either case <f> is independent of K, and it 

follows that A/K in (6.9) is also independent 

of K. Substituting this in L_ = N^G, G„P = -i— ^G,G-P^ 

r r 1 2 At 1 2 n 
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shows the desired zero sensitivity to pure gain changes. 

Similar analysis will show that the zero gain sensi- 
tivity property is retained with the excitation in 

Fig. 6.1 applied to any point between the nonlinearity 
output and the varying gain factor. The inpiit point shown 
has been chosen because it yields system equations in 
convenient form. If, for example, A^ is applied at the 
input to N, the equation corresponding to (6.9) is A/K = 
|a^/K - e^ ( j G 2 ( Pj^ ( jio^) I , and A/K is not con- 
stant due to the A^/K term on the right. 

Now consider the behavior of L^(jo)^) as plant 
dynamics change. From Eq. 6.11 the plot of (f) vs. 
of Fig. 6.3 is constructed. Taking the excitation freq- 
uency ( 1 )^ constant, ( jto^) will vary with plant par- 
ameters and cause to change along its respective A^/M^ 
contour. However, if A^/M^ < 1 there are either two 
solutions, or no solutions for <j>. Fig, 6.3 shows that a 
unique <j) is obtained for every if A^/M^ > 1. 

For this reason A^/M^ >1 is assumed for the remainder 
of the EEAS synthesis. 


Rearranging (6.10) slightly 

-N. Ae^ 

= 2 = 


o o 


[A^ + M e^^l [A /M + e^'*’] 


o o 


This is the relation that fixes the |L^(jw^)| character- 
istic at a particular level. The magnitude of (6.12) is 
shown on Fig. 6.4, and it is easy to show that if 









IfiO 

^ ^ ^ 1 - 120 ®, + 120 ®]. 

From Fig. 6.3 this translates to the same result for 

^ +180° e [-60°, +60°], or in an even larger range 
for larger A^/M^. This means that if is taken close 

to the loop transmission will be less than unity even 

though the limit cycle quenching gain margin p might per- 
mit a larger value. 

From (6.12) we also get | ( jto^) / (1+L^ ( jw^) ) [ = 

ft 

I -el<]>M^/A^ I = M^/A^, which is constant. Thus, as plant 
parameters vary L^(ju)^) varies along a contour of fixed 
[l^ ( jo)) / (1+L^ ( joj) ) I on the Nichols chart as shown by Fig. 

6.5. 



Fig. 6.5 Variation of L^(jaj^) on Nichols 
chart coordinates. 

The above description applies to real, as it is 

for all odd, static, and single-valued nonlinearities, m 

( 

will be assumed real for the remainder of this chapter. 
Although somewhat more tedious, the extension to complex 
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is straightforward, 

6 , 4 Synthesis of an ElemeBtary EEAS 

In this section we shall develop the necessary 
equations and synthesis steps for a simple EEAS where the 
plant is assumed to have only gain uncertainty. The 
approach of Sect. 2 = 1 is follox»7ed x*7ith the intent of 
stripping away the complexity of the most general design 
and exposing the basic design conflicts. 

The structure in Fig. 6.1 is taken with N having the 
usual DIDF, 

= M^A (6.13a) 

= M^/A , (6.13b) 

under quasi-linearity constraints 

”‘^^ix^(t)| 1 Aj^/a (6.14a) 

“b - • (6.14b) 

In addition, we now have the limit cycle quenching gain 
margin 

I 1 P 7 (6.15) 

where satisfies L^(jaj^) = -180®. 

Let the plant be P = KPj^ = K/[s(s+a)] with the gain 

factor ignorance given by K e [Kj^,K 2 l , and and 

a are known. To complete a set of specifications take 
2 2 

T = w / (s+ 0 ) ) , R = q /s, and assume the maximum mag- 
SL a e r ^ 

nitude of oscillating signal at the system output is 
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bounded above by 

The output oscillation maanitnde Ic (iw ) I is 

' o - o ‘ 

written from Fig. 6,1 as 




U + 


(6,16) 


Taking maxiinum plant gain K^, and imposing the magnitude 
constraint gives 


|1 


< m 


(6,17) 


From (6,12) A^/| 1+L^ ( jw^) | = M^/ | ( jw^) | giving in 
(6.17) 

K 2 M^|G 2 (j< 0 o^^h^^%^ I <m|Lo(jWo)|. (6.18) 
The extreme forced input signal to N is 


R T 

V =r 

e W^G^P 
f 2 


Aj^R T 

Vl«2^h ' 


and solving for 62 ^ 


(6.19) 


A^R T 
r = 1 e 

2 M^K-X P, 
f 1 e h 


(6,20) 


Evaluating at and substituting in (6.18) yields 


Fe<i“o> 

^ 2 0 


Af 

— K^mM^ 



( 6 . 21 ) 


This inequality is the counterpart to Eq. 2.9 for the 
SOAS, and indeed reduces to (2.9) if we recall that in the 
SOAS |L^(jw^)| = 1. Thus, the same arguments used in 
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choosing to itkinimize would seem to apply to the 

EEAS, What we have not seen yet for the EEAS is that 
minimizing will optimize the design in the sense of 

noise transmission discussed in Sect, 1,4, 

The constraints on L^(jw) are summarized as 
follows . For limit cycle quenching 


I 1 P ' (6,22) 

where p is now the design value of quenching gain margin, 
presumably smaller than the theoretical quenching limit. 
Due to the limit on output oscillation magnitude 


K-A,M 

i TTifr 

1 f 


Re(jWo)T(jto^) 




(6.23) 


with A^ satisfying quasi-linearity constraint (6.14a). 
The loop transmission magnitude is fixed by 


Lo(j%)| = 


A /M + e 
o o 




(6,24) 


In contrast to the SOAS, the designer must fix |L^(jw^^)[ 
in the EEAS by selecting a suitable combination of 4> 
and A^. This will make the loop transmission shaping 
task more difficult and some cut and try may be necess- 
ary. 

Temporarily choose in the form previously used 

for the Type 1 plant and step inputs, i.e., 

X = 

e * 


(6.25) 



Then inserting all the given quantities 

in the right side of (6.23) results in 


aq K«M . to (s+tjt), ) 

S cl 


( 6 . 26 ) 


3W 


o 


.The bounds for |L^(jw)| are sketched on Fig. 6.6. 



Fig. 6.6 Boundaries for EEAS oscillating 
signal loop transmission. 


The foremost goal in shaping L^(jiii) within the 
constraints of Fig. 6.6 is to optimize in terms of noise 
transmission as discussed in Sect. 1.4, i.e., to minimize 
the high frequency asymptote of the loop transmission. 

We assume that any sensitivity requirements due to param- 
eter uncertainty other than plant gain will permit 
to be below the lower bound in Fig. 6.6. Otherwise, the 
lower bound is not a constraint (in the present problem 
there are no such sensitivity requirements, but they will 
exist in the general design problem). [L^(jm)| cannot 



assume its final asymptote until it has encountered the 

lower bound and has been established along this 

bound. This reasoning indicates that to should be min- 

o 

imized, i.e., the final asymptote should be assumed at 

the lowest possible frequency. One may postulate oo^ 

lower than to_, say at 6to,_ , in Fia. 6.6. In whirb 

11 ~ w -• ■ ■ 

|L^(jto^)| > p. It is difficult if not impossible to 
satisfy (6.24) at this frequency with any reasonable 
when p > 1/2. However, ignoring this, note that the 
lower bound has an excess of two poles (more in practical 
designs), so if |L^(jto^)| > p, [L^(j«) | must then 
decrease slower than the lower bound to avoid having 180° 
of phase lag before reaching the upper bo\jnd at p. 

The conclusion is that the optimum L^(jw) has close- 
ly the characteristic shown on Fig. 6,6. The shaping is 
now explained in detail. As w approaches the boundary 
intersection of Fig, 6.6, |L^(jw)[ lies between the 

bounds, joining the lower bound at the intersection, and 
following it until the far poles may be inserted. The 
poles must be inserted far enough out so [L^(jw^) | < p 
and to permit |L^(jo)^)j to be established on the bound- 
ary. For example, suppose p > 1 and w is at the 

7T 

boundary intersection. Further suppose that we attempt 
to choose 0 )^ = 03^. From Fig. 6,3 <}> - 0°/ and on Fig. 

1/2. Hence the necessary [L^(jto^)[ 
cannot be achieved? instead must be positioned far 

enough above cou- such that (|) differs from zero by the 
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amount required to obtain some |l (io) ) I > 1/2 » Let L 

O " o o 

shown be a trial value and assume that it establishes w 

TT 

at the point labeled. How do we find to^? First, 

recalling from Sect. 6.3 that 7\ A-1 > 1, choose A ~ 

^ o o 

1,2M^ for example. Using this and j/li^(jw) in 

Eg. 6.11 to get for Eg. 6.24, |L*(j(o)| is plott- 

ed as on Fig. 6.7 over the frequency range where the 
trial |L^(jio) I is coincident with the lower bound in 
Fig. 6.6. |L*(jo))| represents the achievable |L^(joi)^)| 

with this specific trial. The oscillating frequency is 



Fig. 6.7 Location of oscillation frequency. 

then chosen at the intersection. If the two curves do 
not intersect may be adjusted, or the far poles of 

the trial loop transmission adjusted. If p < 1/2 (see 
fig. 6.4) it will be beneficial to take lo = 0 ) because 

O IT 

a smaller A^ results. In this case the design is con- 
siderably simplified in an obvious way, but for p > 1/2 
this is not possible unless the optimality of the design 
is compromised by allowing [L^(jw^)| to be smaller than 
required. 
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It is now clear that should be chosen to minim- 

ize the lower bound given by (6.23), and this is accomp- 
lished by maximizing |x (jw) |/A- at w . This is the 
Scime criterion obtained for in the SOAS design, and 

it is treated in some detail in Appendix A. 

The remainder of the design for plant gain uncert- 
ainty and no disturbances proceeds in essentially the 
same fashion as the SOAS design in Sect. 2.1. The selec- 
tion of X^ determines G 2 completely from Eq. 6.20, 
then is shaped to obtain the desired L^. Lastly, 

the prefilter is obtained using the specified transfer 
function from F = T(l+L^)/L^. 

To compare the EEAS design of this section with the 
SOAS design in Sect. 2.1, note that the preceding argu- 
ments have shown that |L^(jw^)| _< |L^(jw_^)j < p. Using 
this in (6.26), and approximating << gives 



(0 

a 


M 






M^K^mp 


(6.27) 


Here it is seen K 2 q^/Kj^mp assumes the role in the EEAS 
which K 2 qj,/Kj^m has for the SOAS in Eq. 2,14. If p < 1 
the EEAS loop transmission must have larger bandwidth. 

For p > 1 it is conceivable that an -EEAS design could 
have smaller bandwidth, however, p must be considerably 
greater than unity to realize any significant advantage 
from the EEAS design. The ideal relay and saturation non- 
linearities have p £ 1.49, which is not large enough to 
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give more than trivial advantage to the EEAS for the 
design problem considered in this section. 

6,5 EEAS Design to Satisfy Quasi-Linearity and Output 
Oscillation Magnitude Constraints with Practical 
Specifications 

Having looked at a rather idealized EEAS example in 
the preceding section, a more general design scheme is 
now developed. The plant will be characterized by 

{ w e W 

(6.28) 

K(w) e [K^,K 2 l, 

where w is the plant parameter vector and only the 
bounding sets W and [Kj^,K 2 l are known. The oscilla- 
ting signal at frequency is supplied by an external 

generator which has frequency tolerance parameterized by 
Y/ with 

“o ^ (6.29) 

A notch filter H is inserted at oscillating frequency 
as in the SOAS, giving the structure shown by Fig. 
6.8. This filter has width Q and it is presumed that 



Fig. 6.8 EEAS structure with notch filter added. 
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Y in (6.29) is adjusted to account for inaccuracies and 
small parameter variations in H. The nonlinear element 
N is characterized by DIDF (6.13) under quasi-linearity 
constraints 

I - (6.30a) 

% - " (6.30b) 

The sub-e again denotes the extreme nonlinearity input 
signal (fastest and/or largest) and the function and par- 
ameter values that exist in the system simultaneously 
with (t) . 

With uncertainty in plant parameters of p, the 

n 

limit cycle, quenching constraint is 

if- <6-31) 

In the previous section we saw that w tends to be less 

TT 

than ( 1 )^, and Fig. 6.5 shows that when this is true the 
left side of (6.31) contains the plant parameter set pro- 
ducing maximum phase lag in the neighborhood of 
This parameter set will be used to define L 

om 

Constraining the output oscillation to have magni- 
tude bounded above by m. 


max 


weW 

L° + J 

max 

IS'3“o'”<i“o>PhO“o)r 

weW 

° 1 
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Variation in the factor 1 P. ( jw^) /L ( jw ) | with plant 

parameters is dependent upon the shaping of 

near w . For this reason the maximizing parameter set 
o 

for (6.32) cannot be determined until the design is comp- 
lete, Further, we want to use (6.31) and (6.32) as 
bounds for [L^(jw)|, requiring that both inequalities 
refer to the same i.e. , both be determined by the 

same plant parameter set. Let contain the plant 

min 


Vhm ■aeflned as follows. Ah<j“o> 

'h 


and Pi. = P, if and only if w e W^C W, Let = 


max 

weW 


K 


m 


hm 

Then take 


m* = m 


m 


K P, (jO) ) 
m hm o' 

max 

"^0 

3 

o 

L (jw ) 
om ^ o 

weW 

KPh(j(^0>^ 


< m. 


(6.33) 


Then 


M K 
o m 


G2(j‘^o)H(ja)^)Phj„(jw^) 


L — 
om o 


< m* 


(6.34) 


defines the lower bound for L^^(jw^) , while (6,31) 
defines the upper bound for L^^(j(o^) . -The ignorance in 
(6.32) has been lumped in the factor m* . The designer 
must estimate this from study of his particular plant, 
and then verify the estimate or make adjustments when the 
first design is complete. 

Continuing the synthesis steps, there is a pairing 
Pj^^ and T^ as discussed in Appendix B that gives the 
extreme forced input to N, 
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X 


X = 
e 


H ' 


. A E T. 
0 © 0 


\a.q T 
b e e 


(6.35) 


In the preceding section it was shown that the optimal 
fxmctional form for is the same as for the SOAS. 

thus, in accord with Appendix A, assuming step commands 
and a Type 1 plant, tb.e form 


X 


el 


^x b 

S+O), 


X 


(6.36) 


is chosen. 0^ will be used to include the necessary far 
poles in X^, 

Solving (6.35) for G 2 and substituting in (6,34) 
yields, after some rearrangement. 


Pel'>o> 


X (jw ) 
e 0 

AjTTjS^ 


A 

e 


K M 

^ ra 0 

R (jw )T (ju ) 
e 0 e 0 


— K m* 

e r 

L (jw ) 
om 0 

'’he'iV 


(6.37) 


This is the generalization of Eq. 6.21, and the EEAS 
counterpart to SOAS Eq. 2.56. 


Eqs. 6.31 and 6.37 respectively provide the upper 

bound for L ( jw ) and the lower bound for L (jo) ) . 

oin oin o 

If no notch filter is used (H = 1) the resultant w is 

■fT 

Shown on Fig, 6.9a. ^om^^^o^ must be situated on the 

lowest (solid) boundary of Pig. 6.9a above u) 

tr 


and this 
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Fig. 6.9 Behavior of and borinds for L as 

e om 

oscillating frequency (o^ is decreased. 

is accomplished by shaping /L^^< joj) and selection of 
as detailed in Sect. 6.4. must be above the 

bound for all e fi. Study of Fig. *6. 3 and 6.4 shows 

that as 0 )^ moves away from } increases, 

thus one need only work with the worst case value of w 
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which is w®/Yo 

If the notch filter is employed we require its mag- 
nitude characteristic. The basic equation for obtaining 
iH(joj) [ is (6.37), however, the value of |l (tw )| is 
unknown. We can assume that |li^^{jw) | is flat over 

and take |L^„(jw^) | = p, yielding from (6.37) 



^ ra*pK 



H ( ju) I = < 

Aag^w, K 
“X ID in o 




1 




; w e 

; tij Q, 

(6.38) 


where = Xaq^Wj^. This should produce reasonable re- 
sults when p < 1/2 and to^ and co^ can be taken quite 
close together, i.e., 4» « 0 on Fig. 6.4. If p > 1/2 
so that 0 must differ significantly from zero, requir- 
ing that A'oTn ^ ^ ^ differ significantly from -180®, 
the approximation |L^^(ja)^) | =: p will deteriate. The 
alternative is to use smaller p or m* in (6.38), but 
how much smaller must be determined by cut and try. 

Given a satisfactory p for (6.38), the problem of 
determining A, and H is identical to that in the 

SOAS detailed in Sect. 2.3, and therefore is not repeated 
here . 


Once 0 )^ has been selected by the method of Sect, 
2.3 and all remaining constraints, sensitivity to Pj^, 
disturbances, etc. permit the use of this value, the de- 
sign proceeds by the same steps detailed for the SOAS in 
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Sect, 2,3, except for the requirement to shape /l (jw ) 
and select A^, X is known from a X vs, to" plot as 
shown by Fig. 2.7, |H(jto) | is given by (6.38) from 
which a rational function is synthesized, and is ob- 

tained from Eq. 6.33. Then G, is used to shape L (jto) 

i om 

in accord with sensitivity and disturbance attenuation 
requirements (treated in the following sections) , and to 
shape ( jto) at as explained in Sect. 6.4. 

Lastly, the prefilter is obtained from F = T (1+L. )/L. , 

0 Z0 f0 

Where 

e 

The discussion in Sect, 2,3 regarding efficient 

assignment of far poles to the $ part of X and to 

X e 

Tg also applies to the present design. 

To make a brief comparison of SOAS and EEAS, first 
note the position of p in H of Eq. 6.38. Comparing 
this H with that given in Eq. 2.57, the EEAS requires a 
deeper notch when p < 1. When p > 1 there is a poss- 
ibility of obtaining a smaller bandwidth loop transmiss- 
ion with the EEAS. The aspect of the EEAS most likely to 
make it superior to the SOAS in some applications is that 
does not change with plant parameters. In the SOAS a 
significant portion of the oscillating frequency interval 
is contributed by the variation of plant dynamics in 
Pj^, and this effect is accentuated as the design is opti- 
mised by reduction of In the EEAS will exper- 
ience this variation, but will not, because this 

frequency is presumably fixed to much smaller tolerance 
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by the excitation source. Only must remain in the 

notch filter in the EEAS, hence the notch can be narrower, 
having less effect on x^{t), and ultimately can be 
moved to a lower frequency than the SOAS permits. There 
is no difficulty if is not within the notch filter 

as long as [L^(joj.^) | ^ p, and this has been insured by 
using the worst case in (6.31). Of course, the require- 
ment of an excitation source is a disadvantage in the 
EEAS, 

6*6 Synthesis of EEAS Loop Transmission to Satisfy 

Transfer Function Sensitivity Specifications 

The transfer function sensitivity specifications for 
the EEAS are assumed in the form detailed at the beginn- 
ing of Sect. 3.3, i.e., we are given bounds on the max- 
imun transfer function variation Aln|T(j(u)[ (and poss- 
ibly A /t ( joi) ) as a function of frequency. The system 
structure remains as in Fig. 6,1, and the plant ignorance 
is as described in Sect. 1.3, Eq. 1.17. The approach for 
synthesizing the forced signal loop transmission 
will be basically the same as developed earlier for the 
SOAS. 

First, it is necessary to write the forced signal 
loop transmission as the product of a compensation trans- 
mission G, and the equivalent plant which contains 

all the loop transmission ignorance, i.e., = GP^. 

We have previously shown in Sect. 6.3, Eq. 6.12, that 
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M. 




M 


[Aq/M^ 


(6.39) 


.3<J> 


where the oscillating input to N is Ae-''*', Thus 

1 


K/A = 


|Gi{ja)^)G2(jt0o)Ph(j(^oH 1^0 \ 


jT' 


and substituting this in the forced signal loop transmis- 
sion, 

Lj = NjGjGjP = ^ KG^GjPi, 


|Gl(j%)G2(ja,^)P^(ja.^)||A^ + 


= G 


|G(j%) I [Ph^^t^o^ I 1*0 


= GK^P^ = GP^ . 


In the last equation the definitions 


(6,41) 


G = G^G2 


Pf " Kf^h 


M. 


Kf H 


G(jw ) I |P (jco ) I [a + M e^*^ 


(6,42a) 

(6.42b) 

(6.42c) 


have been applied. Recall from Sect. 6.3, Eq. 6.11, that 


4> - Tan 


-1 sin4> 

A^/M + cos 
o 


(6.43) 
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As with the SOAS , it is seen that the ignorance of the 
dynamics of introduces a high frequency gain factor 

in the equivalent plant , and contains some 

parameter uncertainty. In the EEAS is set by the 

excitation source, and therefore does not change with the 
parameters of P. Thus we shall assume constant in 

this section, and it follows that there is no uncertainty 
associated with |G(jw^)|. The quantities IPj^(ju^)| 
and |a^ + both vary with plant parameters and 

contain the ignorance of , albeit small for many 

practical designs. 

Now, using the method in Appendix C, we wish to 
calculate sensitivity bounds on the Nichols chart for 
L^(ju). ); i = 1, 2 , n. Such bounds apply to a par- 

ticular plant parameter set. Specifically, we shall 
find bounds for = G( P^^( jw^) , where the 

parameter set is the set maximizing the plant phase 

lag at and i.e., at high frequencies above the 

varying plant dynamics. In the synthesis procedure of 
this section we will take larger than the maximum 

Then the above selection of P^^ gives 
maximum magnitude over the plant parameters, as illus- 
trated in Fig. 6.5. The limit cycle quenching constraint 
is 

Given a plant parameter vector w, we still require 



“o'" ^ order to evaluate at some given 

frequency. Since and G are not availaJDle in the 

first stage of the design, we shall neglect the uncer- 
tainty in and approximate = P^^ for use in 

finding the first set of Nichols chart bounds. Such a 
set of bounds is shown in Fig. 6.10 labeled 



L, ( 310 ) (deg) 


p M./U 


fi: Q_ 




Locus of L-(jw ) 


Fig. 6.10 Illustration of forced signal loop 
transmission shaping for the EEAS. 
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w^, The remaining bounds, B 2 , and B^, 

on this figure are generated in exactly the manner 
explained for the SOAS in Sect. 3.3. 

In Sect, 3.3 (page 82) five properties of the 
minimum bandwidth loop transmission were given. Pro- 
perties 2. through 5. transfer directly to the EEAS by 
replacing in the Chapter III statement with 1*^^. 

Property 1, is replaced by constraint (6.44). Using 
these properties as a guide, the designer shapes 

~ *^^hm ~ *^^fm' thus obtains the first value for 

and the compensation function G. 

Now a value of must be chosen, and it must be 

chosen such that is realizable as a 

solution of Eq. 6.39. At this point we have a plot of 
^fm' shown in Fig. 6.10. The following method 

for selecting is suggested. From consideration of 

Fig. 6.4 and the specified p, choose a reasonable 
say A^ = 1.2M^. Then using the available 
jto) plot the curve of achievable I over 

the frequency range above oi^. On the same figure plot 
®oth plots are depicted in Fig. 6.11, and 
is chosen at the intersection. The technique is the same 



0 ) 



as we have used previously in Sect, 6.4, With (o^ 
selected, L^(jw^) will vary along a contour of constant 
) ( jw) / + L^(jui)) I as depicted on Fig. 6,10. 

All necessary information is now available to cal- 
culate P^(jw) = K^Pj^(jto) without neglecting K^. There- 
fore Pj may be used to calculate a new set of Nichols 
chart bounds for algorithm given in 

Sect, 3.3 (page 85) can be applied directly to the EEAS 
by replacing s\3b-l with sub-m throughout. Using this 
algorithm the designer iterates G to obtain an optimal 
satisfying the Pj^ parameter sensitivity specifica- 
tions , 

As with the SOAS, the design steps for satisfying 
transfer function sensitivity specifications with mini- 
m\am bandwidth are rather tedious. Hence, the designer 
should have some feel for whether the reduction in band- 
width is precluded by other design constraints before 
undertaking the procedure of this section. 

6 , 7 EEAS Design for Disturbance Inputs 

When the extreme forced input to N in Fig . 6.8 
is due to a disturbance input to the system some change 
to the synthesis steps in Sect. 6.5 are necessajry to 
insure quasi- linear operation. We shall denote the 
extreme disturbance input by With this input 

present, there is a plant parameter set P^^^ which most 
nearly causes the corresponding to violate 
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quasi-linearity constraints (6,30). Paired with this 

plant is loop transmission and the plant output 

signal is written from Fig. 6.8 as D L„ /(1+L_ ). This 

e fe fe 

quantity is evaluated at and substituted in place 

of ^^^o^^e ^^^o^ (the plant output for a command 
input) in Eq, 6,37 to give 


Xei<j»>o> 





A 

e J 


K M 
m o 


- K mM- 
e f 


D^.(ju) )L- 
e'-* o fe 






hm 


(j®o> 


l^he‘^“o* 


(6.45) 

All quantities in this equation are defined as in 
Sect. 6.5. We are again faced with the problem that the 
quantity |l^^/ (L^^( l+L^^) ) | is not known precisely un- 
til the design is complete. Hence this factor, at io) , 

J o' 

is approximated by (M^/M^) | l/(l-pM^/M^) | = , 

and (6.45) becomes 



K M. 
^ mo 

Id (jtt) ) 

e o 




-Km* 

e 

M - M- 
O f 




m has been replaced by m* to absorb the ignorance of 
the unknown factor discussed above. 

For disturbance inputs the design procedure to 
satisfy quasi-linearity and the output oscillation mag- 
nitude constraint is precisely as detailed for command 
inputs in Sect. 6.5, except that we use Eq. 6,46 in 
place of Eq. 6.37. The ignorance in m* must be 
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eliminated by iteration. For example, the first design 
is performaed with m* = m, , then if the oscillation is 
too large m* is reduced appropriately and the design 
carried out again. 

As was the case for the SOAS, the EEAS design 
cannot be made to satisfy the frequency quasi- linearity 
constraint (Eq. 6.30b) for step inputs. The SOAS design 
for step disturbances is discussed at length in Chapter IV 
and the comments there apply to the EEAS as well. 

The above relates to satisfying quasi- linearity and 
the limit on output oscillation for disturbance inputs. 
There may also be specifications on the disturbance 
attenuation as a function of frequency. If so, these are 
incorporated into the Nichols chart bounds on forced 
signal loop transmission used in Sect. 6.6. The detailed 
technique for doing this is covered in Sect. 4.2 and in 
Appendix C. 



CHAPTER VII 


EEAS SIMULATION 

In this chapter some results are presented from 
digital simulation of an EEAS, Due to the strong simi- 
larity of the EEAS and SOAS design schemes an EEAS design 
is not detailed starting from initial specifications. 
Instead, an excitation signal is applied to the SOAS de- 
sign of Chapter V and some theoretical predictions of ■ . 
Chapter VI are verified by simulation. 

For reference the EEAS structure is shown again by 
Fig. 7.1. The nonlinear element N is taken as a sat- 



uration with the parameters given on Fig. 5.2. Study of 
the saturation two-sinusoid input describing function 
indicates the limit cycle will be suppressed for p < 
1.49, i.e., the quenching constraint is 

I 1 P = 1.49 . (7.1) 

The plant P in Fig, 7.1 is specified by 
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P(s> 


KPj^(s) 


K 

s(s+a) 


{ a e 

K e [K^^K^l 


[ 1 , 10 ] 

(7,2) 

[ 1 , 10 ] 


with K and a independent. This is the plant used for 
the SOAS example in Eq. 5.3. The compensations P, 

G 2 , and H are given respectively by Tables 5.6, 5,3, 
5.1, and 5.4. These compensations determine the phase 
characteristic of the loop transmission. This is shown 
on Fig. 7.2 for two extremes of the plant pole over a 

band containing the oscillation frequencies of 

interest. 



Fig« 7.2 Phase characteristic of loop transmission. 


The excitation frequency is chosen at = 30 rps. 

Using the SOAS design this choice is rather limited if w 

is to be kept within the notch filter H (shown in Fig. 

2.11). Prom Fig. 7.2 we see that when the plant pole is 

at -a = -10, 0 )^ = 0 )^ = 30 rps. with the oscillating 

component at the input of N written X (iw ) = Ae^*^, 

0 0 
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w = (0 implies that (j> = 0°{see Fiq. 6.3) , = 1«2K 

Q T[ ^ T ri Q Q 

= 3,0 5 is chosen;, and from 




-ej* 

(A^M^ + e^'i’l 


(7.3) 


•this results in lL^(jtb^)| = |L^(j30)[ = .454 . For 

M 

either EEAS or SOAS ^ Gj^G 2 KPj^. For SOAS operation 

jL^{j30)| = 1, while for the EEAS this decreases to .454 
due to an increase in A and a corresponding increase in 
the system output oscillating component. The input to N 
and the system output for a simulation started at t = 0 
as an SOAS, and with excitation A^ applied at t = 4sec, 
is shown in Fig. 7.3. After t = 4sec. -the oscillating 
signal magnitude should increase by 1/.454 = 2.2 as the 
system assumes the EEAS mode. Fig. 7,3 verifies that this 
does occur. The large transients in Fig. 7.3b are present 
because initial conditions were not set to give an immed- 
iate steady state condition. It is evident, by tracing 
•the effect of variation due to the plant pole 

change in Fig. 7.2 through Figs. 6.3 and 6.4, that the 
variation of |L^(jw^)| will be small. 

Before simulating the step response, consider briefly 
the sensitivity of the present EEAS to parameter changes 
of Recall for both systems we write the forced sig- 
nal loop transmission as = ^f^h* Neglecting the small 

variation in K^, the EE2^ loop transmission is just the 
SOAS loop transmission of Chapter V reduced in magnitude 
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by 201og(.454) = -6.8db. This magnitude reduction 
causes the EEAS loop transmission to violate its respec- 
tive sensitivity boundaries on the Nichols chart. Using 
the the boundaries are computed in accord with 
Sect, 6.6 and Appendix C and shown on Fig. 7,4. Para.-,., 

meters of the EEAS are tabulated on Table 7.1. The 

boundries apply to ^fm and it is seen that they are 



Fig, 7.3a Input signal to nonlinear element 
in SOAS and EEAS modes. 
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Fig. 7.4 Bounds for using 

G of Table 5.2b and A. = 1.2M = 3.05. 


"o 

^f 

<}> (deg) 

K 

o 

n 

1. 

-30. 

12.7 


.99 

-23. 

4.2 

30. 

.99 

-16. 

2.5 

30. 

1. 

-9. 

1.8 

30. 

1.02 

0. 

1.3 


6.66x10' 

6.60 

6.59 

6.63 

6.78 


Table 7.1 Parameters of EEAS 
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violated by up to 6.5db. Since jLjp(jw) ( is smaller for 
all to in the EEAS it should be expected that the step 
response will be slower (it could of course be speeded up 
by changing the prefilter F) . Further, since the sensi- 
tivity bounds are not satisfied, the EEAS step response 
can be expected to change more with 

EEAS step responses for two plant pole extremes are 
shown in Fig. 7.5. These may be compared with similar 
SOAS responses (on a different time scale) in Fig. 5.11a 
and 5.11b. Some time parameters from the two systems are 
collected on Table 7.2, and these figures support the 
qualitative predictions made above regarding the differ- 
ences expected in the two system responses. 


System 

a 

ft 


SOAS 

10 

1.1 


SOAS 

1 

,45 


EEAS 

10 

2.2 


EEAS 

1 

.55 



Table 7.2 Comparison of SOAS and EEAS 

step response characteristics. 














CHAPTER I IX 


SUMMARY AND CONCLUSIONS 

•8.1 Sxarttnary of the Oscillating System 
Synthesis Philosophy 

It has been found that the design of either of the 
oscillating systems considered herein is optimized with 
respect to the effect of sensor noise on the plant by 
designing for the minimum oscillating frequency Lower 
bounds are imposed on by several design requirements. 

In particular, these are: 

a) the frequency guasi-linearity constraint which 
requires that the forced signal describing function 
component for the nonlinearity be valid over the 
bandwidth of the fastest forced signal input to the 
nonlinearity, as well as over the bandwidth of the 
fastest system output. 

b) the combination of magnitude quasi-linearity con- 
straint and the upper bound on oscillating signal 
component at the system output. 

c) achievement of the required sensitivity reduction 
for ignorance or variation of plant parameters, 
excluding the gain factor. 

d) achievement of a specified disturbance attenuation. 
Techniques have been developed in detail for determining 
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the lower bound on oscillating frequency for each of the 
factors above. The designer must first ascertain which 
requirement dominates in his particular design, i.e., 
which yields the greatest lower bound on oscillating 
frequency. Whenever possible, this should be done with a 
rough approximate design to avoid tedious optimization 
with respect to all of the above requirements. Once the 
dominant requirement is determined, the lower bound on 
oscillating frequency is minimized using the techniques 
and considering the tradeoffs detailed in the preceding 
chapters . Portions of the design emerge as a result of 
the oscillation frequency minimization procedure, and the 
remainder is relatively straightforward and explained 
herein. 

8,2 Concluding Comments 

The oscillating system stands out as one of the few 
nonlinear adaptive schemes that have been put to practical 
use. Heretofore, there has been no synthesis theory for 
designing such systems to meet practical quantitative 
specifications. In this paper a synthesis theory is de- 
veloped which proceeds from quantative specifications in 
a transparent fashion to a design that is optimal in the 
sense of sensor noise power at the plant input. All sig- 
nificant constraints on the design are considered, and the 
synthesis theory is demonstrated with numerical examples. 

There is a distinct class of design problems for 
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which the oscillating adaptive scheme is superior to a 
linear adaptive design. The specifications and optimiz- 
ation criterion used in this work are sufficiently close 

(221 

to those used for linear adaptive design •' to permit a 
valid comparison of the best possible design by both 
methods. Thus, it is now possible to quantitatively 
evaluate which of the two design methods yields the 
superior solution to a given adaptive problem. 
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APPKNDIX A 


iH 


Nonlinearity Forced Input Fimction 

In this appendix the nonlinearity forced input 
signal is considered in some detail. The goal is to 
choose this signal, subject to the appropriate con- 
straints, such that the oscillating system loop trans- 
mission bandwidth is minimized. The structure of 
Fig. A.l is assumed and the forced input to N 
denoted by X^. 



Fig. A.l Oscillating system structure. 

X^ m\ist satisfy the quasi linearity constraints 

™^^j x^(t) I ^ Aj/a (A. la) 

% — (A. lb) 

which were introduced in Sect. 1.2. In Eq. A.l, is 

the oscillating frequency, is the bandwidth of 

X^Cjw), Aj is the minimum magnitude of oscillating 
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signal component at N, and a?3 are the quasi-lineayity 
parameters of N discussed in Sect, 1,2. In Sect. 2.1, 
Eq, 2.9, it was shown that must satisfy 

U(u^) = (jtOj,) Xj^(t) [J > [U(jo)^)|/ (A. 2) 

where U is completely determined in terms of given 
system parameters and specifications. In the frequency 
range suitable for to^ , ]u(jw)j is decreasing with u 
so that the optimum X^, which will permit minimiam , 

is the function that maximizes y in (A, 2). In Sect. 

2.2 it is shown that the left side of (A. 2) can, in an 
idealized sense, be made arbitrarily large by choosing 
- X^j^ + where X ^2 given a constrained 

maximum effect on the time response and a complex pole 
pair at , allowing any value for 

hence for |x^(jw^j)|. However, this X ^2 impractical 
because of the acute system sensitivity to small 
parameter variations. The design of a practical X ^2 
is described in Sect. 2.3. In this appendix we want to 
determine the function X^^, which is to be the rela- 
tively 'smooth' low frequency dominant portion of X^, 
that will minimize the remaining burden on X^ 2 » 
allow the minimum if X ^2 = 0- 

For notational convenience we will consider 
X^2 = 0 and = X^. 


In order to deal with specific 
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functional forms for some assumptions must be made 

with regard to the system input and the plant. We 
assume the plant has one pole at the origin (Type 1) and 
R is a step input. Using the Pinal Value Theorem and 
Fig, A.l, c(«) = c(t) = sC(s) = 

s^Q sN^X^ (s) (s)P (s) . For c(«>) finite in the pres- 
sence of step inputs, X^(0) must be finite and non- 
zero, excluding poles and zeros at the origin in G 2 . 
Summarizing, X^(jtjj) is to be a relatively smooth 
function having bandwidth and with Xj^(O) 

finite. The simplest such function is 




(A. 3) 


which yields 


'0/1+1/32 


(A. 4) 


Since X^ = RFGj^/d+L^) with = N^G^GjP in 

Fig. A.l, it must have the same excess of poles as RFG^ 
3), thus X^ in (A. 3) can not be realized in the 
assumed structure. However, sufficient far poles 
(>> cOg) can be added to make X^ realizable with no 
significance to the present discussion. On Table A.l 
below several additional possibilities for X^^ which 
allow closed form calculation of ja are tabulated. 

These do not of course cover all the possible choices 
for X^, but by varying the several parameters in the 




s+1 


^ + 1 


2 

1 ? n > 1 


(s+l)"^ 

(n-1) 

3 

; 1<P<U)« 
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5 
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/u/+l /w.^+p ^ 
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fl < P^' < w, 

/to. ^+P 1 ^ 4.K 
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table a wide variety of functions are considered. 

For the purpose of comparing the functions in 
Table A,1 define as the solution of 

|Xf ( jWb)/Xf (0) I = 1//2 . For 3 = 3, Trial #1 gives 
= .316 with = 3, For Trial #2 with n = 2, 

Vig = .575 with Wg = = 3(.644) = 1.93. Thus Trial 

#2 seems to be better than Trial #1 by the ratio ^ 2 /^ 1 = 
1.82 (= 5.2db) . This is true under the stated condi- 
tions, but it is found that for > /e^-1 = 

2.53. Thus for practical oscillating frequencies Trial 
#1 is superior to Trial #2 with n = 2, i.e., other 

design considerations nearly always force 
For n = 3, for > /(e^-2) /2 = 1.64. It is 

clear that in practical design Trial #1 will provide 
larger y than Trial #2. 

Investigation of the remaining functions of Table 
A»1 is quite detailed. We shall simply summarize by 
stating that no acceptable combination of parameters in 
Trials #3, 4, and 5 produce a y that is larger than 

than of Trial #1 with reasonable to . 

0 

From the above it is concluded that X^(s) from 
Eq. A. 3 (Trial #1) gives the largest y that is readily 
available and this function is used as a standard for 
Type 1 plants and step commands in the body of the paper. 
It may also serve as a standard against which any new 
function not considered herein may be compared. 
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For plants with other than a single integrator 
and/or other than step inputs, an investigation similar 
to that above should provide the designer with sufficient 
insight to choose a reasonable X^. 



APPENDIX B 


Deteriidtiatioii of the Extreme Plant 
and Transfer Function 

Using the oscillating system feedback structure on 
Appendix A Fig. A.l, the forced component of signal input 
to N is 

. AZ. 

where Z^ is the forced component of plant output, A 
is oscillation magnitude at the input of N, and = 
M^/A is taken as the forced signal describing function 
for N. The plant is characterised as 

P(s,w) = K(w)Pj^(s,w) ; w E W, (B.2) 

where w is a parameter vector and only its bounding set 
W and possibly constraints among its elements are known. 
Pj^(s) is taken such that Pjj(s) -► 1/s^ as jsj ^ In 

order that the assumed describing function hold, the 
quasi-linearity constraints, 

I < A/a (B.3a> 

% 5 (B.3b) 

are imposed, oi^ is the approximate bandwidth of X^(jw), 
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the oscillating frequency, and a, 3 are the quasi- 
linearity parameters discussed in Chapter I. (B,3) must 
hold for all system inputs and all w e W. Therefore, we 
require that these inequalities hold for the worst case 
which is written 


X 


e 


■A T7 

e e 


^f^e°2^e 


(B.4) 


The sub-e denotes that and the accompanying function 

and parameter values which most nearly violates (B.3), 
i.e., the fastest and/or largest X^. In the oscillating 
system synthesis theory it is assumed that the permiss- 
ible are given by specifications* From this range 

of Z^, and the possible Pj^, the designer must select 
2^ and the accompanying that do indeed produce the 

extreme X^. This selection must be made before the de- 
sign is carried out, so G 2 and the precise pairings of 
Pj^ and Z^ are not known. 

It is admitted at the outset that we do not know how 
to make the above selections in the case of an arbitrary 
plant and specifications on Z^. What is done below is 
to make some observations that give the designer a rears. ■:!. 
sonable chance to make a good selection. If this select- 
ion is incorrect, it will be revealed by the first comp- 
leted design, and the designer should by then have the 
insight required to make the correct selection. 

For a command input Z = R T , and for disturbances 
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. M. 

it is ^fe = ~\^1^2^he * 

© 

Presumably the expected range of inputs can be surveyed 

to determine R or D_. Calculating the ratio X-/X 
e e j. c 

while the command input is held at its extreme value. 


. A 

' R PL . ■ 
s f 



r 2 n 

L f J 

. . A^., 

■ 




and canceling all the coimtion factors 


^f = 


1 + L 


fe 


1 + L, 


X. 


(B.5) 


The same manipulation for a disturbance input yields 
exactly the same result. Therefore, , is the same 

parameter set regardless of whether is the result of 

a command or a dist\irbance. For command input we 

can also write 


AK 

X - — ^ 
f A_K 


" f 

phe 


P, 

^ e 



X 


(B.6) 


The variation in AK^/A^K with plant parameters is of 
the same order as the variation in ^ ' 

which tends to be small because is well above the 

varying plant dynamics and this plant ratio approaches 
unity for large to . Now suppose temporarily that the 
transfer function sensitivity is to be quite small so 
that f/T^ 1 for all plants. Then 
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(B.7) 


Consider a plant Pj^ whose poles and zeros vary indep- 
endently, e.g., with IPj^(ju) | lying in the envelope of 
Fig, B.la. We have postulated in Fig. B.la that the ex- 
treme plant is P^ with zeros set to minimum and poles 

J. 





logo) 



(b) 


Fig. B.l Illustration of extreme plant selection. 

set to maximum. The concurrent variation in X^- in 

r 

accord with (B.7) is shown in Fig. B.lb. As Pj^ dev- 
iates from |x^(jw) j becomes smaller than |x^(jw) j 

for all w, but the bandwidth of X^(jco) increases. It 
is inferred from this observation that |x (t) | will 

likely have the largest maximum value, but x^(t) will 

0 
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be the slowest (“t) • Thus the plant that gives the ex- 
treme X^(t) with respect to (B,3a) is different from 
the plant giving the extreme with respect to (B,3b). 

This observation is supported by the following example. 
Let = l/[s(s+a)l, a e [1,10], and = 

l/[s(s+a^)]. Let = 150/[ (s+6) (s+25) 1 . In Pig. B.2 
the various x^(t) from Eq. B.7 are shown for three 
choices of a , In Fig. B.2a, with a = 1, x (t) is 
simultaneously the fastest and smallest x^(t). 
Conversely, with a = 10 in Fig. B.2c, x (t) is sim- 
ultaneously the slowest and largest x^(t). 

Recall in the synthesis steps where X is used 

e 

(Sect. 2.3 and 6.5), that the designer is most concerned 
with constraint (B.3a). This constraint is used directly 
in the equations of that synthesis procedure to obtain a 
lower bound for the oscillating frequency. Once the low- 
er bound is obtained, (B.3b) merely dictates whether it 
may be used in the design. Since it is not generally 
possible to select P^^^ to give the extreme x^(t) for 
both (B.3a) and (B.3b) simulteineously , it is clearly 
appropriate to choose the extreme with respect to (B.3a), 
i.e., the magnitude quasi-linearity constraint. The pos- 
tulated Pj^g in Fig. B.la is the correct selection when 
poles and zeros vary independently. For the more general 
plant, a good choice of Pj^^ seems to be the parameter 
set that minimizes |Pj^(jw) | over the maximum portion of 
the system bandwidth, obtained by inspection of |Pj^(jw)| 
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over the permissible plant parameter sets. 

Moving now to the case where the transfer function 
variation is not vanishingly small, it is quite reason- 
able that the above selection of P, remains valid. 

he 

Then it is necessary to determine which of the admissible 
transfer functions T is the that is paired with 

With selected as above (zeros minimized and 

poles maximized) it tends to have smallest magnitude and 
least phase lag over the system bandwidth. A glance at 
the Nichols chart, which displays the relationship be- 
tween T/F s= L^/ (1+L^) and = N^G^G 2 KPj^, immediately 

leads to the conclusion that T is the T for which 

e 

|T(jw) I is minimum over the system bandwidth. 

In summary, the guidelines for choosing Pj^^ and 
are as follows. Pj^^ is the plant with minimum 
[Pj^(jm) I and maximum , while is the lower 

bound of specifications on ]T(ju))|. 



Appendix C 

Derivation of Bounds for L(j<>)) front Transfer 


Function Sensitivity or Disturbance Response 
Specifications 


The technique for determining bounds on acceptable 
loop transmission from specifications on transfer func- 
tion sensitivity and/or disturbance response is detailed 

below. The technique applies to any linear time- invariant 

[181 

two degree-of-f reedom feedback structure'' ^ and is 
taken from [22]. The structure in Fig. C.l is used, 
however, all two degree-of-f reedom structures are equiv- 
alent with respect to the specifications considered 



.Fig. C.l A two degree-of-f reedom feedback structure. 

herein; thus any other structure could be used, or could 
be transformed into the structure shown. For the purpose 
of providing insight into the meaning of the bounds for 
L(jw) and how they are constructed a graphical method is 
first presented, then a digital computer implementation 
of the method is given. 

In Fig. C.l the plant P(s) is the constrained part 
of the system given as a rational transfer function, but 
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wit±i bounded uncertainty of some parameters, i.e., P(s) 
= P(s,w); w e W, where w is a parameter vector, and the 
bounding set W is known along with any constraining 
relations or correlation that may exist among the 
elements of w. The loop transmission is 


— n. ( ci \ “o ( ts T.T^ 


/ ~l \ 

V'-o x; 


Since there is no uncertainty associated with G, the 
plant variation may be written at some specific fre- 
quency as 

Aln[P(jo)^)] =-- Aln[G(j(Oj^)P(ja)^) ] = Aln[L(jWj^)] 

= Aln|L( 1 + jA y4.( ju). ) . (C.2) 

The transfer function is 


C _ a, _ FL 
R 1' + L 


(C.3) 


There is no uncertainty associated with F so that varia- 
tions in T{jw^) are written as 

Aln[T(jto^)] = Ain = =-j (C,4) 


= Ain 



L{jo]^) 


L{ jw^) 


1 + L( jo). ) 

X 


1 + L(jw^) 


Thus both the given plant uncertainty and the specified 
permissible transfer function uncertainty are expressed 
in terms of L, Suppose that the maximum permissible 
transfer function uncertainty is specified by 

Aln|T(jo))| < ln[3(u))]; 3(o)) >1. (C.5) 

We shall show how this specification is translated to a 
boundary of acceptable L(jto^) at the frequency 
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First a number of discrete parameter vectors 
r = 1, 2f m are selected and the plant is evaluated 

at for each, giving the complex numbers P_^(j«i,Wr); 
r = 1, 2, ... m. These points are plotted as shown by 
Fig, C.2 on Nichols chart coordinates (ln|p(jw)| vs. 



/p ( jtri) ) . The points are then enclosed to define a region 
of plant variation as shown, which we shall refer to as 
the plant template. The designer must select m large 
enough and the w^ with sufficient care to insure that the 
entire region of plant variation is included in the temp- 
late with negligible error. Since L(jw^) will vary with 
the plant parameters, we will deal with a particular Wj^ 
which gives Lj ( jco^) = G( jto^)P^ ( . To find a point on 
the boimdary for Lj(jw^) the template may be cut out of 
Fig. C.2 and laid on the Nichols chart for L(jw) shown in 
Fig. C.3, (which has been drawn using the same scales as 
the plant template) keeping the coordinate axes through 
Pj ( parallel to those of the L(jo)) chart at all times. 
The L(jo)) chart has contours of constant 
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ln| L( jw) / {1+L { jto) ) j , Choose any angle /g ( , or equiv- 
alently /L^{jo)^) , and position the plant template on 



Fig. C.3 Method for locating one point on the 
boundary of acceptable Lj(jw^). 

this angle grid. Next the template is moved vertically 
along the chosen angle grid until one finds a point where 
Eq. C.5 is satisfied with equality, i.e., until the temp- 
late just fits between two ln| L ( jw) / ( 1+L ( jw) ) j contours, 
say ln6 and ln63- (3. = 3(o).)) as shown by Fig. C.3. In 

this position the origin of template coordinates deter- 
mines a point on the boundary for L^(jw^). Additional 
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points are found by repeating the process at different 
angles. 

A plant template is prepared and the boundary for 
L^{jio) found at a sufficient number of discrete frequen- 
cies to provide a high degree of confidence that (C.5) 
will be satisfied closely at any co by an Lj satisfying 
the bounds at the discrete frequencies. Of course the 
same plant parameter set is used for Pj at each frequency. 

The Nichols chart also has contours of constant 
/l ( j oj ) / ( 1+L ( j bj) ) , so that given a specification of the 

form 

A ^T( jo)) < ’?{(*)), (C.6) 

the boundaries for ( jw) to satisfy this specification 
can be found in the same manner as above. For any ( jtOj^) t 
the greater of the ln[Lj(jw^)| determined by (C.5) and 
(C.6) determines the bound that must be used to satisfy 
both specifications. 

Next we consider the bounds on L(ju) determined by 
disturbance response specifications. From Fig, C.l the 
disturbance transmission is 

^ = T = -A- = lA = _V_ . (C.7) 

D D 1+L 1+iA “ ,1+V 

We assume the specification on is given as 

ln|Tjj( jw) I < ln[p(w) ] . (C.8) 

Here we have a specification on the absolute value of the 

transmission rather than the change. However, (C.8) must 

be satisfied for all values of the plant. To find the 
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boundary for ( jw^) we first construct a template of 
l/P(j(o^) since AlntV(jo)^)] = Aln[l/P ( ju^) ] . Positioning 
this template on a ln[V(jw)] vs. /v( jca) chart, a boundary 
for is found as demonstrated in Pig. C.4a. This 



Fig. C.4 Location of boundary points for Lj(ja)^) 
to satisfy disturbance transmission 
specification. 


boundary is then translated to its equivalent for Lj(jWj^) 
= l/Vj(jo)^) as depicted by Fig. C.4b. 

Next a digital computer implementation for finding 
the boundary for due to transfer function spec- 

ification (C.5) is investigated. Let 
G(jWj^) = ge^® 

Pj^{jw^) == p^e^'^r- j- = l,2,...m. 


(C.9) 

(C.IO) 
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Using Tg, *3.ue to respectively in (C.5) with 

equality and squaring both sides yields 


Ti (jw.) 

2 

Pi~ 

2 

( gP 2 ) ^+2gp 2 cos ( 0+c|) 2 ) + 1 

tTTd^ 




Jgpx) ^+2gpjCOs(0+(|)j)+l_ 


= f(g) = {C.ID 

It can easily be shown that f(+“) •> 1, f(g) 0 , and f(g) 
is finite for all real g and 0+({>^ Rtt. Also, f'(g) = 0 
has two real roots giving the qualitative sketch in 
Fig. C.5. Given P^, P^, and 9, it is a simple matter to 



Fig. C.5 Graphical representation of Eq. C.ll. 

so.lve for the roots g^, g^ of (C.ll). Consider the poss- 
ible results in determining the boundary for ( jw^) = 
gp^e-'' ^1". 

g^ , g^ complex: For this case there is no real g 

satisfying (C.ll) and therefore at the given ( j<*^|) 
= 9+4>j there is no boimdary point on the Nichols 
chart of Fig. C.3. This case can occur if the per- 
mitted variation in transfer function is larger than 
the plant uncertainty. 

gj, g^ < 0: This gives two boundary points on Fig. 
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C,3 on the angle grid at { jw^) = 0+(j)j^+18O® , Thus 
= 0+<j)^ there is again no boundary for 

{ jw^) . 


^ ^2 ^ negative root gives a boundary 

point displaced 180® from the angle grid of interest 
at 0+(j>^. g^ gives a valid boundary point for Lj ( jw^) . 

9 2 ^ This case produces upper and lower 
boundary points for ( jto^) on the angle grid 

( j^j^) = Li(jw^) must lie below the upper 

boiondary or above the lower boundary. 

The above observations allow us to construct an 
algorithm for obtaining a boundary for Lj{jw^) as follows. 
Algorithm for bounds on ( jw^^) due to Eq. C.5 

1. Select plant parameter sets w^; r = l,2,...m e W 
such that the plant template boundary is satisfac- 
torily covered. 

2. Choose 0 = ^(jtOj. ) and solve Eq, C.ll m^-m times, 
i.e,, for all combinations of |t (joj.)/T (jo).)|; 

r 7^ s, getting roots g^^, g^^; ct = l,2,...m* < m^-m. 
Let 


q* = "'in {a } 

l<a<m*^yia-^ 


_ max f , 
^2 l<a<m*^^2a^* 


^la ^2 a if g* < 0 there are no 

bounds. Otherwise, the lower bound for L (joj.) is at 

If g^ > 0 there is also an upper boujid 
. When there are two bounds the lower 
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bound is higher on the Nichols chart (g* > g*) and 
L^(joj^) must avoid the region between the two 
bounds . 

3. Repeat 2. for as many 8 as desired to obtain the 
complete boundary for . 

Treating specification (C.6) in the same fashion 

= y., (C.12) 

i' 

and after some manipulation 

g^PiPa [sin((}) 2 -<(i j) -tcos ]+ (C.13) 

g [Pj (sin (9+(|)2) -tcos (0+(j)2) ) -p^ (sin(9+<f) j ) +tcos (6+(}> ^) ] -t = 0, 
where 

t = tan('Fj^+4>2"<f>i) • (C.14) 

The equation for g locating bounds for is again a 

quadradic. The bounds can be calculated by using the mag- 
nitude algorithm given above with Eq. C.13 substituted in 
step 2 . 

Lastly, the disturbance response specification of 
Eq. C.8 yields 

(gp) ^ + 2gpcos ( 9+<J>) + 1 - 1/p? = 0. (C.15) 

The bounds for this specification are calculated using 
the same algorithm with (C.15) in step 2., which only 
needs to be solved m times because (C.15) contains only 
one value of the plant. 

Listed below is a FORTRAN subroutine which calculates 
the composite boundaries for Lj(jw^) due to the three 
specifications (C.5, C.6, and C.8). A representative 
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output, is shown using the plant P(s) = K/[s(s+a)l with 
K e [1,100], a s [1,10], and specifications as listed 
on the printed output. Also listed are PLANT subroutines 
suitable for evaluating the forced signal equivalent 
plant Pj used in the SOAS and EEAS iterative designs for 
parameter sensitivity and disturbance attenuation. 



nnoooooonooooonrs 
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SUBRCIJTIKT N!CRNn<W,RrtA, A1 A2pPS) 

J, N. SI'^AY, RrVlSFO NOV, 72* 

RPr. NOROl^lT? ANn SIDI, international JOURNAL OF rONfROl-* 

VOL- 16, N0. 2, 1.972, PP, 287"309- 

RCtiTINE I nCATFS ROUNDS TOR G AND L1-G*P1 AT FReOUENcV U<RPS) 
SATISFY the most RESTRICTIVE OF 

1 . h^AXCABSf T|JW> ) )/MIN(ASS<T< JW) ) )=PE J Ad ) K 

2. r-AK(ARG<T(jWn |-MlN(ARGfT(jW)))aflE1Ar2><DFGl ^ 

3- KAX( A0Sf TD( JW) nrRETA<3) ^ 

ir SETAM) is Nr(^ATIVE THE SPECIF I f A T ION IS IGNORED- ^ 

USER RUST SUPPLv SUBROUTINE PLANT WHIrM EVALUATES A 

Pd.DsPLANT magnitude AND PI2,I)5PLANT ANGLE IN RADIaNs N 

FOR N parameter SETS- ^ 

PS IS SWITCH WHICH aLLOwS USER ROUTINE TO RFaD PLaNt PARAMETER K 
DATA ON first CALL OF PLANT. N 

BOUNDS ARF rOMPi^TED FOR ANGLE OF LI ON (A1,A?> aT INCREMENTS K 

DA IN CEGRFES, ^ 

DIMENSION P(2,50),orTA(3) N 

1 FORMATdhj ,10X,*NI^MOLS CHaRT BOUNDS AT FREOUENCt = # • Fl0 - a, /IN ,10»N 

2,*PETAa)a*,F8-4,# PE T A ( 2 > = • , FB • 4, # PE^ A ( 3 ) =•. F 8 . 4 . HR, N 

310X<«RCUND G-ANGLF»,4X»*G-HAG*,3X,«G*MAG(f)BM.2yF 
4*Li-ANGLF«,3X,#Li'*MAG*,2X»#L1-MAC<DRI*//> 

2 FORMATaH ,t5X,Fin,4,22X,F1 0.4, # NO SOLUTION*) 

3 F0RMAT<1H .10X,*LOwFR«,7F10,4> 

4 FORNAT<tH , 10X . »l)POER« , 7F10 • 4 ) 

5 rORMATdh ,i0X,»STAqILITY LIMIT ENCOUMTEREO* > 
ir<PS-EQ-0.) CALL PI ANTCP.W.NtPS) 

IF(PS.EO.0.) return 
MHITE( 6d) W.BET.A 
BETA2=B£TA(i)*»2 
0CTAlal.-BETA2 
0ETA3 = 1.-1 ./RETA<3)<»«2 
BETA4=3.1415927*0ETA(2>/ie0. 
call PLANT(P.W,N,PS> 

AlaAl-10C .«P( 2.1)/.^. 1415927 

10 Al5Al*3.14J5927/18 

glm=-i. 


FUNCTION magnitude 


n))MPd,i) 


GUM-1 ' ^ 2g 

ir(9ETACl>.|T.0-) GO TO 20 

CALCULATE bounds FOR TRANSFER 
DO 15 1=1, N 
DO 15 J=1,N 
IF(I.EG.M> GO TO 15 

C1=(PC1^^ >«C0S< A1 ♦P(2. I >)-RETA2*P(1, I ) *C0S ( A3 ♦P ( 2 ► 

2#P(l#J)*eETAl) 

C2s(P(1#^)«*2-(BETaci>#P(1, I))**2>/< (P < 1 . j ) *P< i , m **2*0 ETa1) 

C3=Cl**2-C2 

IF<C3.LT.0.) GO TO 15 

Gl5-Cl-SORT{f:3> 

G2»-2.*C1-G1 

GUMsCl 
GLMaG? 


IF(Gl-LT.GLM) 
IF(C2.CT.GLM> 

15 CONTINLE 
20 GLAs-1. 

GUAsl«E2e 

IF(BETA(2),LT.0,) CO TO 
CALCULATr rounds FOR 
DO 25 Isl.N 
00 25 J=l,\ 

IF<I.EC-j> GO TO 24 


30 

transfer 


FUNCTION angle SPcC I F T CATION . 


N 
N 
N 
N 
N 
N 
K 
N 
N 
N 
K 
N 
K 
N 
N 
N 
N 
N 
N 

SPECIFICATION .N 
N 
N 
K 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
N 
K 
K 
N 
N 
N 
N 


ICRP008 

|CRi^0i e 

ICR'/01 5 

irpt70?e 

ICP00^^5 

ICRt^0?F 

irR(>03^ 

1CPV'04P 

JCRV045 

1CRV056! 

lCRf''05S 

1C0E706P 

irRt^06« 

ICRP 07C 

irRJ>07? 

ICByOR^ 

ICR^ 

irR{^09^ 

IC0V0C*: 

ICRi>10^ 

IC0pl0? 

ICRPllC 

icR(?iie 

icR'^iac 

rrRj.i?® 

ICBP13^ 

irp(^i35 

ICRP14F 

irPf7i4« 

lCPr’151? 

Irppl5f 

1C0V<16( 

ICP.^165 

ICPV17C 

110^^17? 

irenpp 

ICPV192 
ir.B 10? 
IC0^^2P|2 
1CR720? 
lC0v2d 

TCRV??« 

irP(?225 

lrPt^23C 

irRv?3? 

ICP‘^24? 
ICBv 259 

Kfl<^25? 

iraF26P 

irflP279 

IC0E274 

tCRrPPF 

ICR92P? 

lCP(^29e 

ICRP29? 
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TeT«ri<BCrA4<.P(3.UoP{;>.J>> 

C0=P<l.n«P(’ .J>»<ciN(P(?,n-P(3,J))-T«C0S(Pr?,|>>P(?,J>)) 
Cl = <Pn.no(41N(«1«P(2. l}>-TeC0S(Al*P<3.n>>°P(l,J>« 
2fSIN(Al«P<?..tn«T«rns( tl<>P(2.jm)/42.«CP) 

C2a-T/C8 
CSsCjee 2 »ff 2 
tF(C3.t.C.(l,t GO TO 25 
Cls-Cf-SC»T(r3> 

C23*2>*C1oG1 

]F(G1 .1.T.GLA) GUAaGl 

|F(G2.GT.GLA> GLAsC? 

25 COftTtMtE 
30 GI.D3-1. 

GUD31.C20 

rF<eETA(3).LT,0.) CO 70 4» 

C CALCULATr BOUND' FOB OlSTUReANCE RF<:PONSF SPFCIFI CAT ION . 

DO 3b lal.N 

Cl«C0S<Al»F«2.in/P»l.l> 

C2sCl •»2-BETA3/P(l ,n*o2 
IFIC 2 .LT. 0 .) GO to 35 
Cls-C1-sCBT(C2) 

C2s- 2.»C1-C1 
IFCGI .LT.GLD) GUOsGl 
lF(G2.GT.GLn> GLOaG? 

3b CONTINUE 

40 GLsANAXl(Ct.N,r,LA,GLD) 

GU«AMlM(GLH,GI;A.GMDt 
4l8Al»J8e,/3. 1415927 
AGP8A1*180,»P<2.1)/3. 1415927 

ISafl 

11.80 

IFfGL.LE.O.) GO TO 5;: 

GLDBa22.«aLOG10(GL) 

CLPsGL*Ptl.l) 

CLPOBs20.«ALOG1O<GI.P> 

MB|TE(6>3t At .GL.GLDB. AGP,GLP>GLPOe 
C CHECK IF present LOWER BOUND IS UNSTAPLE FOR ANY PLANT. 

DO 45 I a] .A 

ALaAl*180.*Pt2.1)/3. 1415927 
GPaGL*P(l,I) 

IF(AL.LT.(-1H0. ) .AND.GP.GT.I. ) |Lal 
' 45 continue 
CO TO 55 
90 lSaIS«l 

55 irtGl.LE.B.) CO TO 6 'T 
GUD0a2<;.<A|.0GlP<GII) 

GUP*CU«Pfl,l) 

GUPD8a23.«ALOGi0(C'*P> 

mriteia. 4) At .GU.cune.aCP.GUP.Gupofi 
GO TO 65 
80 tSatS«l 

IF{IS.EQ.2 ) wRITE(6 ,2 » A,,agP 
65 IF(IL.EQ.l) WR1 Te(a,5> 

IF<IL.EQ.1 > RETURN 
IFCAGP.lt. A?) RETURN 
AlaAi-EA 
GO TO 10 

end 


MCB«*308 
MCfl-'30' 
MG8K31T 
MCB?315 
AtCfiP32e 
MCRF325 
A IC«»33T 
A ICBf)335 

A tnrYrtTfji/? 

• » I L n ^ ‘J •! 1 C, 

MCR(0345 
KICB7352 
MCBr.'35* 
MCfl?3^Cr 
Mf = 
MCR/37(' 
Mr0^37e 
MCR»738e 
MC«5^3e5 
MC0P39C 
MC0039« 
MCRt^4?)P 
MCflt^405 
MCP?41C 
MCP«»415 

MrRC^4?5 
MCP/43P 
MfRl»43l» 
MCRf44IZ 
MCfil^445 
MCPtJ^45C 
^ rCP0*455 
MC9P46P 
MCPfr465 
MCP«^47Cf 
MC0:^475 
MCBlC4eg 
MCPP485 
MCR^"49C 
UfR049“ 
MCR(^50I2 
Krcp«^505 
MCB^510 
MCP«^519 
MCR?52e 
MrP(»525 
Mrfl053C 
MCP«53? 
MCBP940 
MCB054" 
MCR^550 
MCP(.555 
MrBPi5ft0 
Mrpt^56« 
MCPt^57? 
MCPf/575 
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SUBftauTINE PLANHP.H.N.PS> 

C ROUTINE GlTNFRATrS VALUES OE PLANT WITH WARYINO GAIN. POLES. 

C AND ZEPCS rOR CALCULATION OF NICHOLS CHART ROUNPS. 

PlMCNSION PfP.S0).n<tf I.Pecit’T.PmBnZCPI.ZriJ.lOJ .PC(2.10> 

1 FORMATT 13.7x.7E1P.4/(8E10.4/>) 

2 FORHAT<lhl»10X.«PLANT PARAMETER VALUER licEO TO DETERMINE NICMOLS C 
2HART 0CUM)S*/1H0.1 x.*ERROR C0NSTANTS*/<'H ,10y,6E12.4) ) 

3 FORMATTIH ,10X,*PO( F*.6F10.4/tlH , 1 7K . AFLO . 4 / ) ) 

IFIPS.EO.I.) GO TO 3«5 

REAL ANC write plant PARAMETER VALUES. 

add one read ANn ONE WRITE CARD FOR EACH VARYING PARAMETER. 
REaOCS. 1) NG. (Cm . fsl.NG) 

MRITE<6,2){g< 1 >.Isl .NO) 

RCAD(S.I) NP1.(PP(j). J sl.NPl) 

WRITET6.3) {P0|J),Jsl.NPl) 

C SET UP constant parameters of plant. 

Pi(i)=ii. 

NP:2 
NZ=0 
NPC=0 
NZCS0 
RETURN 
3D N30 

C ADD ONE no 60 FOR EACH VARYING PARAMETER. 

DO 60 1=1. NG 

DO 60 Jsl.NPl 

Pl(2)sP0«j) 

«oN»l 

call 0CDE(W.G{I»/P1 < 2 ).NP,Pl.NPC.PC.NZ. 2 .NZC.ZC,PR.PI. 0 e.P( 2 ,N>» 

60 P<l.N)sS0RT<PR»*?*Pl**2) 

RETURN 

end 
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BOf’P*W»(5*NP.P,NPC.PC.NZp?,NZC,Zn,flH,6l,D0, 4> 

subroutine evaluates rational function at jw. 

ZEROSsZU). complex P0LESs7ETA,0MFGArPcn.n,PC(?.I 

I). ERROR CONSTANTsG, 

DIMENSION Pll).2(l),PC(2,10),ZC<2.10) 

PU3.1A15927 

VkR _ A > ^ 4^ 

UU.'*U^ 

A«0. 

IFINP.LE.P) SO TO 3 
DO 2 1*1, NP 

IF(ABS(PI1)),LT»1.E-10) CO TO 1 
DBsDB*P(! )••^/(H••9*P(|)••2) 

A*A»ATAM|M/Pn ) > 

GO TO 2 

1 OBaOe/He«Z 
AoA-Pl/2. 

2 continue 

3 IFlNZ.LE.p) GO TO A 
DO 5 Isl.NZ 

lF<A8S(Zn)),LT.l.r-10) CO TO 4 
D8»DB*(«**2*ZI I l••2)/^(I)••2 
AaA*ATAN(H/Z( I >J 
GO TO 9 

4 0Bs0B*N««2 ^ 

A»A*PI/2. 

5 CONTINUE 

6 lF(NPC.LE.e> GO TO 9 
DO 8 1=1, NPC 

D0rDB*PC(2,n**4/nPCI2.|)**2-M..2>«.2*4.a(Pn(Z.|).pr(i,n,M>..91 

A.A-ATANC?.*PC(i.|).PC(2.n,H/(PC«2.n..?-W..?r» ‘ ^ » 

ir<w,cE.PC(?,nt GO TO a ^ ^ 

ir(PC<lpI),LT,0.) GO TO 7 

A«A*P1 

GO TO 8 

7 A«A*PI 

8 CONTTNUE 

9 ir(NZC4LEp0) GO TO 12 
DO 11 IclpNZC 

DBsDS^C CZC(2p I )«« 2«W«*2)**2^49#(ZC(l8 I)»7C(?* 1 /7nfo. i 

A-A*ATANC 2 .* 7 CClpr)# 2 C( 28 l).N/C 2 Cf 2 .n* 4 -W*» 2)1 

]F(NiLE#ZC(2p n> CO TO 11 
1F(2C(1^ I) ) GO TO 10 
A«A*Pl 
GO TO 11 

10 A»A-P! 

11 contznle 

12 O0sSQRT(A8S(DB) > 

GRsOG«COS{ A ) 

GrsD6*SZN(A) 

DB»209«ALOC10(DB) 

RETURN 

end 


ponE 5*000 
FOrE‘ Bi? 'i 
>PODE''010 

ponF70?0 
ROrE(^0:>5 
porFf 0050 
fiODE''0;^5 
PODFv’04fl 
|3ODE('049 

ponrt^0*;0 
POrFF055 
eODFf'0f 0 
pOnFl '065 
FOtFf ^070 
tOrE{/075 
ROCF{?080 
ROf)E</0P9 
RODEc^a^B 
FOrFi^099 

POpEra0*5 

pore^ijB 

POCEini^ 

POHEPISB 

eonE^i?5 

POOEi>i;^0 

eODFP135 

RODE^UB 

F0rEn49 

eOHEt^lSB 

eOCFr>l 58 

eOfiEraaB 

BOrEt'165 

eorE0i70 

ROCE'?175 

BOrE'^lPR 

BOrEf'lfl9 

ROrE(/190 

G0DEP195 

BorFr'200 

BCDEE209 

eODF?'2l0 

P0PE7218 

FODEe2?0 

GODE0225 

PODE0230 

feOrE(^239 

FOnE^?40 

P0DE^24B 

POn|T 17250 

FODEt^259 
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PLft'-T PA-«'MFTirR vftWu--S JS<^D TO OFTERMINE WICHOUS CHAHT BOUNDS 


FPROR rO.sST 


loOOOOF*On 1 




POLF l.POOO }«nOOO 

SoOOOO 

7.0000 

10.0000 


NirrtOUi C»'A'?T bOi, 

'' S AT F 

REQUEMCVs 

• 1000 





= •loOOOO 

BETAOln - 

l.OOOO 


ROUND fi-AM'5LF 


G-MA6(00) 

Li-AN6LE 

Ll-MAQ 

U1*MA6(0B) 

(.OWE? 

^.UTR 

R.<5<S03 

-lOO.UOOO 

31.3223 

29.9171 

LDWER -^<»,9AR4 

1 J.Tl 5? 

??.7A22 

■>120.0000 

136.4490 

69.6990 

LOWER -44,?qR4 

•!<*.4654 

27.77U 

>140.0000 

243.4401 

47.7278 

LOWfR 


30,0416 

-160.0000 

316.1690 

49.9964 

LOWER --i4,i>qq4 

3b,0f>q.? 

30.89S3 

-160.0000 

349.9412 

50.3S50 

LOWER -llta,?p44 

3‘».0017 

30.6300 

-200.0000 

333.3291 

50.5868 

STA^ILiry LlnTT FNCO.fUTF.R 

"0 




NICHOLS Cma»T HOI' 

•s= >s at f 

'»F0II£NCY= 

.1000 



BFTAd )= -1 ,of 

EtA(?)s 1.0000 

RETA(3)= - 

1.0000 


BO'JiMD G“AM(3i.F 

3-V4AG 

0-MaG«D8) 

L1-AN6LE 

LI -MAO 

Ll-MAO(De> 


i-<^wER 

-4, 

67,1 1. 4 4 

35.1360 

-100.0000 

568.3139 

55.0918 

LOWER 

-94.9R94 

52.3539 

34.3790 

-120.0000 

520.9419 

54.3350 

LOWER 

-44 , 

I , E7B4 

39.5145 

-140.0000 

410.7355. 

59.2712 

LOWER 

-64.9R94 

26.2240 

94.0363 

-160.0000 

950.9883 

47,9931 

LOVJ^R 

-H4.9R94 

R.126S 

15.7443 

-180.0000 

60.9613 

35.7011 

LOWER - 

104,9994 

15.6436 

93,8867 

-200.0000 

155.6595 

♦3.8435 

StaItI TTY LI'^rT 

CrjCJiiNTE RPD 





NICf>OL-S 

1 Chart wou''OS at f-»e3Uencys 

.1000 



9FT4 n ) 

s -l.nnno 

HEtA t?) = 

-1,0000 

RFT4r3)a 

.5000 


ROUND 

G • A M G L, ? 

AG G 

-MA6(08) 

Ll -angle 

Ll-MAG 

Ll-MAG(Oe) 

LOWER 


1.8190 

5.1965 

-100,0000 

18.0995 

25.1533 

LOWER 

-?>a,?A94 

2.9029 

6.8599 

-120.0000 

21.9197 

26.8167 

LOWER 


2,5757 

8.2180 

-140.0000 

25.6294 

28.1748 

LOWER 


.'.',859R 

9.1267 

-160.0000 

28.4561 

29.0835 

LOWER 


2.9941 

9.5254 

-IttO.OOOO 

29.7927 

29.4822 

lower - 


2.9509 

9,3970 

-200.0000 

29.3554 

29.3537 


STAdlLiTT tl'fliT E.MCOtiwTtRFO 
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SUBROUTtME PLANT(Pp,n<N,PS> 

e ROUTINE CENERaTPS NORNALIZEO plant VALUES PF FOR ROAS NlChOLS 

C CHART SOUND ITERATION. 

REAL KF.KS.kINF 

DIMENSION Pp<2,90J.oa2>.P0fl0I.Pl(10».Z(10>.Zr(P,10>,PC(P,ie> 
DIMENSION NF(5S> 

1 rORHATII3.7x.7E10.4/<8El0.</>> 

2 fO«MAT<lHl,10X,*PLANT PARAMETER VALUES USED TO DETERMINE NICHOLS C 
2HART BuUNdS«/) 

3 FORMATIIN ,10X,«POLE*.6F10,4/<1H ,17*.6F10.4/> > 

A F0RMATI412I 

5 FORMATtlM0,i0X,«*0»t H0 ITERATION FAILED TO CONVERGE IHI*> 

6 FORHAT(lM0.i0 X, •compensation PARAMETERS*/) 

7 FORMATI8E10.4) 

8 FORMAT(1H0,10X.«RE4L poles*. /<1H ,10X.F10.4,)> 

9 FORMAT(1N0.10X.«REAL ZER0S«,/(1H ,10X,F10.4. > ) 

10 FORMAT<lM0.10Xi*COMPLEX POLES(ZETA.OMEGA)«./(lH . 10X , 2F10 .4. ) ) 

11 FORMAT(1M0.10X.*COMPLEX ZER0S(ZETA.0MECA)«./<1H .1MX.2F10.4. )) 

12 FQRMAT{lM0,14X.*H0«,ex.*Kr«,8X,eK0*.8X,*KlNF*/) 

13 FORMATUH .10X.3P10.4.E12.4) 

IFIPS.EO.I.) GO TO »e 
GHal. 

Gl»l. 

MO. 09 
N9b1 

c read and mRite plant Parameter values. 

hRITE(6.2) 

READ(9,1> NPl.{P0(j).Jal,NPl) 

HRITC(«.3) (P0(J), jal.NPl) 

Pl(l)e0, 

NPo2 

NZ*0 

NPCO0 

NZCO0 

C READ FIRST GUESS OF OSCILLATING FREOUFNCV Ha. 

REA0(9.7) H0 « 

c read and write compensation data. 

READI9.4) NPL.NZL.NPCL.NZCL 
hR1TE{6.6) 

if(npl.eq.np) go to 20 

I0sNP->l 

RCAD<9.7) (P1(I).|bI0.NPL> 
do 19 1«I0,NPL 
15 6lsCl*Pl(I) 

wAITE(6.8) (P1(I).I>I0.NPL) 

20 tP(NZL.EO.NZ) GO TO 30 

I0BNZ41 

READ(9.7) {z(I>.I«I0.NZL> 

DO 29 1*10. NZL 

29 CI*GI/Z<1) 

HR1TE(6.9) czm,|a|0,NZL) 

30 IFINPCL.EO.NPO GO TO 40 
I0eNPC*l 

REA0<9.7) (PC(l.l).PC(2.l).|el0.NPCL) 

DO 39 1=10. nPCL 

39 CI»G1*PC(2.I)**2 

MRITE(6.10) (PC(l.n.PC(2.II.I«l0,NPCL> 

40 IPINZCL.EO.NZC) GO TO 50 
10sNZC*l 

READI5.7) (ZC(l.n.ZC(2.I). I«I0.NrCL) 
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DO 45 l=j0rN?CL 
45 elaCl/(ZC(2. 

NRIT£«6,11> (ZC<l,n.ZCC2,n.I = 10.NZCL> 

90 N°0 

DO 70 Jal.NPl 

pn(2>speu} 

NoN*l 

CtLt. BODE(h,l,/Pl<?>,NP*Pl»NPC.PC,NZ.Z»N7C.ZC.PR.PI.l^0.PF(^*^) > 
IP(NS.€0.()) 60 TO Tp 

C ON FIPST call solve FOR MB AND ADJUST GAIN SO LF ( JW0)s-l/J . 

A0al000. 

DO 55 K«l,25 

call 80DE(Wffl,GM/Pi(2>.NPL.Pl.NPCL.PC»NZL,Z.NZrL,ZC,HR,Hl ,DH/Ai) 
call B0DE( .9S«M0.eH/Pl{2).NPL.Pl.NPCL.PC.NZL»Z*NZCL . 2C . FR,F I • OP . A2 
2> 

AlaAli-3. 1415927 
A2aA2*a.i4j5927 
M13H0>.05*V0«A1/(A1*a2> 

IP<AeS<Al-A0j.LE..'l.AND.4BS(M0-Wl).LT..05*M0> GO iO 
A0aAl 
55 HgaHi 

MRITe<6,5) 

60 KF(N»8,5ZSCRT<HR*«?*HU«2) 

K«sGPaKF(N)/Pl(2) 

KlNPa6I*K0«Pl<2) 

IF(N.Ne.l) GO TO 45 
GMsGM*KF(1) 

KOaGM 

KlNFa6I*K0*Pl(2> 

KP(l)al. 

HRITE<6«12) 

45 HRITE(fi>13> U0>KF(Nt >K0»KlNF 
70 PF(l.N)asORT(PR**2*PI**2)*KFCN) 

N6a0 

RETURN 

end 
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SUBHOUTthC PI.ANT«Pf,H*N,PS> 

C ROUTINE GEMERaTFS NORMALIZED PLANT VALUES PF FOR EEAS NICHOLS 

c chart bound iteration. 

REAL KF,K0,KINF 

DIMENSION PF{2.50>.O(l3l,P0a0).Pl{i0),Z(i0),2cr?,10>,PO(2,je) 
DIMENSION KF(S0> 

1 FOB«AT<l3,7w,7El0,4/(eci0.4/n 

2 FORMAT<IH1,10X.«PLaNT parameter values used to determine NICHOLS C 
2HART BOUNDS®/) 

3 FORMATdH ,10X.«POl F«.6F10.4/flH , 1 71 ? , /ir la, 4 / ) » 

4 F0RMAT(4I2) 

5 FORMAT(iH0,i0X.«W0s*,ri0.4,® A0/Mns*,Ei»,4) 

6 FOflMAT(lH0,l0X,*COMPENSATION PARAMETERS®/) 

7 FORMAT(8E10.4) 

8 FORMAT<IH0.10X,*REAL POLES®, /(IH ,10X,F10.4, ) ) 

9 FORMATtlH 0 .i 0 X.*REAL ZEROS®. /tlH .l0X,Fj0.4, > ) 

10 FOflMAT<lH0,10X,®rOMPLEX POLES <ZET A .OMEGA ) ®. /< IH . 1 /X . 2F1S . 4 , ) ) 

11 FORMAT<lH0,i0X.orOWPLEX ZEROS ( ZFTA, OMEGA ) o, /< IH . IPX , 2F1 g , 4 , ) ) 

12 FORMAT(lH0,i2X,®AflG<LF)®,3K,®MAG(LF)®,5X,*PHI*,7y,*KF*,flX,*Kff®, 
28X.6KINFO/) 

13 FORMATUM , 10X , 5F17 . 4 , E12.4 ) 

IF(PS.EQ.l.) GO TO 50 

GNsl. 

Cl»l. 

NS^I 

C READ AND WRITE PlaNT PARAMETER VALUES, 

WRITE( 6 . 2 ) 

RCADI5.1) NPl.(P0(J),Jal,MPl) 

WR1TE(6.3) <P0(j) , jai.NPl) 

PKDsB. 

NP»2 

NZs 0 

NPCs0 

NZC®0 

C read oscillating FReQUENCT W0 and A0/M0. 

READI5.7) H,40 
WRITE(6.5) W.A 0 

C READ AND WRITE COMPENSATION DATA. 

RCADI5.4) NPL.NZL.NPCL.NZCL 
WRITE( 6 .S) 

IF(NPL.EO.NP) GO TO 20 
iBaNPol 

HEADC5.7) (Pl(I),laI0,NPL) 

DO 15 1*10, NPL 
15 Gl»Gl®Pim 

WHITE( 6 , 8 ) (Pl(l).|aI0.NPL) 

20 1F(NZL.E0.N7) go TO 30 
I0aNZ«l 

REA0(5.7) (Z( I )• lalP.NZL) 

DO 25 I«I0,N7L 
25 GJsCI/ZU) 

MRITE{6,9) (Zn>.l = I0,NZL) 

30 IF(NPCL.EO.NPC) GO TO 40 
I0aNPC*l 

R£AD(5.7) (PC(i,I).PC{2,I).I»l0,NPCL> 

DO 35 1*10, NPCL 
35 Cl=GI*PC( 2 , I)*«2 

WRITE ( 6 , 10 ) (PC(l,n,PC(2,I), Ial0,NPcU 

IF(NZCL.EO.NZC) go TO 50 

I0=NZC*1 
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C 

C 


C 

C 


«5 


READ«5,7> I>.ZC(2»I>o1=I0*N?CL) 

DO ^5 IaI0.M?CL 
GRaG!/rZC(2, I )**2) 

MR I TE ? e 0 1 1 ) ( ZC ( 1 . n . ZC < 2 » n . I S 1 0 . NZCL > 
N“0 

DO 70 4=l«KPl 
PJ,{2)5P0U) 


CAU B0DE(M,l./Pl{?>pNP»Pl,NPC,PC.N2.Z.N7C.ZC.PR,Pl.De»PF(2*K > > 
IP<NS.EO.01 GO TO 71' 

ON FIRST CALI. SOLVE FOR PHI AND ADJUST GAIN SO 
LF( jH0)a-.5*EXP‘(PHI)/(A0/M0<-EXP(PHl ) > . 

CALL 0ODE<h.GM/Pl(?> ,NPL,Pl»NPCL.PC»NZL.7,N7CL,7r.Hfi.HI fDR. AL ) 
PHI0SPH1 (AL. A0> 

GL0a.'i/SORT(A0*«2+2.*A3*COS(PHl0»*l. > 
KF(N)=GL0/SQRT{HR**2*HI*»2) 

K0aCMeKF(N)/Pl(2> 

KlNF«Gl«»<0fiPl(2> 
lF(N.NE.n GO TO 6<5 
GM=GM*KF(i) 

K0=6M 

KINF»GI*K0«PK2J 

KF(1)=1. 

MRITE(6»J2) 

A5 PH10=PH10<»1b0./3. 1415927 
AL=AL*180./3. 1415927 

WRITE (6»13) AL.GLO,PH!C.KF(N),K0.KINF 
70 PF<l,N)=SOflT <PR«®2*P1*«2>*KF{N> 

NSs0 

RETURN 

END 

FUNCTION PNKAL.A0) 

FUNCTIOA SOLVES fO. 6.11 FOR P«I(RAD>. GIVEN 

AL = ARG(L0{JW0)HRAO) and A0SA0/M0 BY SECANT METHOD. 

1 FORMAT{iM0,i0X, ‘FUNCTION PHI FAILED TO CONVERGE Pe.=*.El2.4, 

2* Pl=*>E12.4,» ria«,El2.4./> 

Pi=4.*ATAN(l. ) 

P0=1.5*(AL*Pl > 

Pici,55« ( AL*P I > 

F0sP0-ATAN(SIN(P0>/<Ae>COS(P0))>-AL-Pl 
DO 2 lal.25 

Fl=Pl-ATAN(SINfPl>/< A0*COS(P1) > I-AL-PI 
pHIsPi-Fl*< (Pl-P0>/(Fl-F0)) 

1F(A0S<F1).LE..01.AND.ABS<P1"P0>.LE..01) RFT'IRN 

F0=F1 

P0aPl 

2 Pl=PHI 

WRITE<6»1> P0.P1.F1 

STOP 

END 


c 



